Dynamic light scattering method and apparatus

Abstract: A dynamic light scattering apparatus comprises at least two lasers (11) of different wavelengths which illuminate the same volume (SV) of a fluid containing scattering or fluorescent centres. The scattered or fluoresced radiation is detected at each discrete wavelength and properties, such as the concentration of different particle sizes, of the volume (SV) are derived. The detection at multiple wavelengths allows the effects of multiple scattering to be overcome. Multiple volumes (SV) may be illuminated and detected simultaneously so as to measure the spatial variation of particle size concentration. The rotational diffusion constants of the particles may be measured from the depolarized component of scattered light. The apparatus or method may be used to measure the velocity of particles in the volume (SV). The apparatus or method may also be applied to remote sensing of the atmosphere by scattering of laser pulses, the intensity variation functions of the received scattered pulses being concatenated prior to auto- or cross-correlation to detect particle size concentrations.

Description: The present invention relates to a dynamic light scattering (DLS) or photon correlation spectroscopy (PCS) method and apparatus, particularly but not exclusively lidar applications. PCS/DLS involve the illumination of a sample with coherent light. The intensity of photons scattered from the sample fluctuates as a result of the Brownian motion of scattering particles in the sample. The number of photons scattered in a narrow range of angles, defined by specific apertures and light propagating and focusing elements, is detected as a function of time and the resultant function is autocorrelated, preferably in real time. This allows the distribution of relaxation time constants and therefore the distribution of sizes of the scattering particles to be determined. DLS-Lidar involves the active remote sensing of a portion of the atmosphere by illuminating the portion with coherent light and detecting the intensity of photons scattered from that portion as a function of time. The resultant function is autocorrelated in real time, in order to derive the particle size distribution within the sampled portion of the atmosphere. A problem associated with the above techniques is that, as the number per unit volume of scattering particles increases, there is an increased probability of multiple scattering in which a detected photon has been scattered by more than one particle, which means that the desired properties can no longer be measured directly. Furthermore, the range in velocity of Brownian motion becomes restricted because of collisions between scattering particles. In order to compensate for these effects, the normal practice is to detect the scattered photon intensity at a range of different scattering angles, typically from 10° to 180° (backscattering). As there is a quadratic dependence of time constants on scattering angle for the ideal case, any deviation from this quadratic function, due to multiple scattering or restricted particle collision, is quantifiable and can be compensated for. However, such procedures are time-consuming because of the number of readings that need to be taken, and require expensive goniometers. In lidar applications, it is not practically feasible to perform angle-dependent scans of a portion of the atmosphere, so that the particle size distribution and degree of multiple scattering cannot readily be measured. US 5,155,549 describes a dynamic light scattering apparatus comprising a laser source for directing a laser beam into a probe and a photodector for converting the scattered light detected by the probe into electrical signals. The probe includes at least one optical fiber connected to the laser source, and a second optical fiber connected to the photodetector. Each of the fibers may adjoin a gradient index microlens which is capable of directing a collimated laser beam into a scattering medium. Thereby, the use of a fiber directional coupler can be avoided by assigning one or more optical fibers for transmitting an optical field to a scattering medium and using several optical fibers for receiving the signal at different scattering angles. US 5,208,651 describes an apparatus and a method for measuring emission and lifetime properties of fluorescence radiation emitted from a sample. A source of a multiple-wavelength electromagnetic radiation is passed through a first interferometer modulated at a frist frequency, the output thereof being directed into a sample to be investigated. The light emitted from the sample as a result of the interaction thereof with the excitation radiation is directed into a second interferometer modulated at a second frequency, and the output is detected and analyzed. An object of the present invention is to provide a simplified dynamic light scattering apparatus and method as well as a simplified fluorescence correlation spectroscopy method and apparatus, respectively. These objects are accomplished by a dynamic light scattering method according to claim 1 and a dynamic light scattering apparatus according to claim 15. According to one aspect of the present invention, there is provided a dynamic light scattering apparatus and method in which a sample is illuminated by laser light of different wavelengths and the fluctuations in light scattered or fluoresced by the sample at each wavelength are detected. The time constants of particles in the sample are derived from the detected light fluctuations by auto-correlation or cross-correlation of the detected light intensity with respect to time. Measurement of light scattering or fluorescence at different wavelengths replaces measurement at different scattering angles, thus simplifying the process of measurement. According to another aspect of the present invention, there is provided a method and apparatus for measuring the velocity of particles in a sample by means of dynamic light scattering, in which the frequency of oscillations in correlation functions of the detected intensity fluctuations as a function of time are measured and used to determine components of velocity in the direction of detection. In an embodiment of the invention, a dynamic light scattering apparatus comprises discrete numbers of transmitting lasers with various wavelengths ranging from the visible to the near infrared. Each discrete laser transmitter forms a transmitter channel, which operates in conjunction with one or more receiver channels. The laser light propagates to the light scattering sample via apertures, focusing, telescope and/or fibre optical elements. By these elements, beam waists are produced in the sample. Corresponding to the channels of transmitting lasers are a second set of functionally identical elements, for collecting and receiving each wavelength of scattered light separately and simultaneously. The apertures, light collecting and propagating elements for one channel are adjusted in combination to produce an intercept of the dynamic part of the autocorrelation function close to unity. The automated functions of the apparatus from one back scattering position replace the function of a goniometer for scattering angle scans. This is achieved by (quasi-) simultaneous, n-fold channel operation. The apparatus enables classification and quantification of multiple scattering sample systems, the normalisation of particle cross-sections, sizes and shapes to the wavelengths used and the quantification of repulsive forces exerted to particles or restricted particle collision as in a gel-like system with cage functions. If multiple scattering is encountered the apparatus allows adjustments to be selected to derive the true particle size. The simultaneous three general fields: particle sizing in the atmosphere by active remote sensing, particle sizing and visco-elastic property determination in fluids and suspensions as for example in crystal growth of bio molecules from solutions and property determination of light emitting samples by fluorescence recording as for example in immune assay developments. The apparatus utilises cross-correlation, depolarised and polychromatic dynamic light scattering in pulsed and continuous wave operation and filtering methods to derive appropriate autocorrelation functions. Specific embodiments of the present invention will now be described with reference to the accompanying drawings, in which: Figure 1 is a schematic diagram of the optical front end of one channel of an apparatus used to measure backscattering from suspended particles; Figures 1a and 1b show alternative examples of the optical front end in Figure 1; Figure 1c depicts different arrangements of optical fibres of different types; Figure 2a shows an example of recorded intensity fluctuation as measured by the apparatus of Figure 1; Figure 2b shows the autocorrelation function of the fluctuation of Figure 2a; Figure 2c shows the particle radius distribution derived from Figure 2b; Figure 3 shows the variation of time constants with scattering angle for different particle sizes in a sample of aqueous solution or fluid suspension with a given viscosity and refractive index; Figure 4 shows the variation of time constants with scattering angle for different particle sizes in a sample of air, gas or atmosphere with a given viscosity and refractive index; Figure 5 shows different autocorrelation functions and bulk velocity measurements from particles in samples of fluid at different velocities; Figure 6 shows autocorrelation functions and corresponding radius distributions from samples with contained particles having a relative velocity to each other, showing oscillations of autocorrelation functions prior to nucleation of these particles to form a crystal; Figure 7 shows the measured angle-dependent deviation of decay time constants from a theoretical curve as shown in Figures 3 and 4 with the given viscosity and refractive index; Figure 8 shows the variation of scattering volume length with atmospheric height at different separations between transmitter and receiver; Figure 9 shows the equivalence of time constants and wavelengths with scattering angle as calculated and normalised for two different wavelengths; Figure 10 shows a generalised embodiment of the present invention for the measurement of a sample; Figure 10a shows beam tracings in an implementation of the embodiment of Figure 10; Figure 10b shows an implementation of the embodiment of Figure 10 for measuring the depolarised cross-correlation DLS; Figure 11 shows the measured values of time constants for two different wavelengths; Figures 12a and 12b shows two alternative embodiments for active remote sensing of the atmosphere; and Figures 13a to 13f show the signal processing steps for measuring particle radius distribution in the embodiments of Figures 12a and 12b. Figure 1 shows in schematic form the optical front end of one channel (i.e. one transmitter and one receiver) of the apparatus consisting of a transmitting optical fibre or fibres 4 connected to a laser source and a receiving optical fibre or fibres 2 connected to a detector, such as a photomultiplier tube or avalanche photo-diode, mounted together so that the axes of their ends are parallel (as in Figure 1a) or at a defined angle (as in Figure 1b) to each other. Laser light emitted by the transmitting optical fibre(s) 4 diverges over a beam B while the receiving optical fibre(s) 2 collects light emitted from a reception volume RV. In a scattering volume SV which is the intersection of the illuminating beam B and the reception volume RV, particles scatter the laser light in different directions and the receiving optical fibre(s) 2 collects light scattered within a narrow range of directions. As shown in Figures 1a and 1b, the defined angles between the transmitting fibre(s) 4 and the receiving fibre(s) 2 are adjusted by movement of tight-fitting coaxial micro-positioning stages (MPS) in the reciprocating direction shown by double arrows, so as to locate the scattering volume SV at distances L1 and L2 in Figures 1a and 1b respectively, and to compensate for the effect of Snell's law on different wavelengths when co-locating the beam waists of different wavelengths. Additional beam forming is produced by suitable optical elements (OE) such as lenses, prisms, wavelength filters and/or telescopes, according to the specific application. The receiving and transmitting optical fibres 2 and 4 are selected as either mono-mode or few-mode (i.e. more than one but preferably less than about two modes) fibres and may either receive and transmit through the same optical elements or have different optical elements. Different possible combinations will now be discussed with reference to Fig 1c. In arrangement I, the transmitting fibre 4 comprises a single mono-mode fibre and the receiving fibre 2 comprises a single few-mode fibre (shown as option Ia), or the transmitting fibre 4 comprises a single few-mode fibre and the receiving fibre 2 comprises a single mono-mode fibre (shown as option Ib); in either case, the transmitting and receiving fibres share the same optical elements. In arrangement II, the transmitting fibre 4 and the receiving fibre 2 have separate sets of optical elements. In option a (indicated by the reference IIa in Figure 1c), the receiving fibre 2 is a few-mode fibre and the transmitting fibre 4 is a mono-mode fibre, as in arrangement Ia, while in option b (indicated by the reference IIb in Figure 1c), the transmitting fibre 4 is a few-mode fibre while the receiving fibre 2 is a mono-mode fibre, as in arrangement Ib. In arrangement III, the transmitting fibre(s) 4 and the receiving fibre(s) 2 have separate sets of optical elements. The transmitting fibre 4 is a single mono-mode fibre, while the receiving fibres 2 comprise at least two mono-mode fibres in parallel, each having a separate detector, discriminator and pulse former. Preferably, the receiving fibres 2 are arranged in pairs of fibres arranged to receive the same wavelength, with the outputs of each pair then being cross-correlated. This cross-correlating arrangement eliminates detector dead-times, after-pulsing and dark counts and will be described in more detail with reference to Figure 10b. The detector outputs a current proportional to the rate of photons received by a receiving fibre 2, which current varies with time as shown for example in Figure 2a. The current is converted to a digital value by an A/D converter connected to a computer, which stores a set of digital values as a function of time. Alternatively, individual received photons are counted by means for example of a photomultiplier tube, a discriminator for setting a threshold for the signal of a detected photon, a pulse former for generating a pulse when the photon is detected, a counter for counting the generated pulses within a repeating period and a storage device for storing the count values for each period. The computer performs an autocorrelation function on the stored values as follows: where I(t) is the fluctuating received light intensity as a function of time t, such that l(t) equals the average scattered intensity over the integration time, and τ is the sampling interval, which forms quasi-logarithmic channels in the structures of the correlator. A sample autocorrelation function is shown in Figure 2b. The ACF function is related to the particle sizes by the following equations: ACF (τ) = g (2) = 1+β· e (-2Γ·Δτ) where g (2) is the second order laser light intensity fluctuation, "l" is an arbitrary baseline, β is a quality factor having a theoretical maximum of l, Δτ is the entire sampling time structure of the correlator, and τ C = l/Γ = 1/( q 2 × D ) where τ C is the time constant, D is the diffusion constant and q is the scattering vector, where q=(4π/λ)sin(Θ/2) with Θ being the scattering angle and λ the wavelength of the impinging light. The diffusion constant D relates, via the Stokes-Einstein equation, a geometrical factor, the viscosity η and the hydrodynamic radius r h to the thermal energy, which is Boltzmann's constant k multiplied by the absolute temperature T as follows: D = kT 6πη r h Under ideal conditions, the time constants τ C are directly related to the radius r h and to the viscosity of the suspension of scattering particles. The term β (also called the intercept) is inversely proportional to the number of observed coherence areas (number of speckles) N, which fluctuate due to Brownian motion. The measured ACF is formed from many different individual functions of the form of equation (2), for each of the different particle sizes. The conversion of the composite ACF shown for example in Fig. 2b to the relative amplitude intensities for different particle sizes shown in Figure 2c may be performed by a suitable computer program such as is described in 'CONTIN: a general purpose constrained regularization program for inverting noisy linear algebraic and integral equations', Provencher, S. W., Computerphysics Communication 27, pp. 229-242, North Holland Publishing Company, 1982. As a result of the above relationships, DLS can be thought of as a point-like measurement technique with virtually no spatial dimension. Three dimensional (voxel) resolution requires an embodiment using a plurality of channels. The spatial mode of the received light defines N together with its inverse proportionality to the number of scattering centres in the scattering volume. The configuration of the apparatus used in this embodiment includes a minimum of two mono-mode receiver units, in order to obtain the maximum intercept value and to enable signal processing with cross-correlation for the various polychromatic embodiments. The cross-correlation function is given by: where I(t) is the intensity measured by one receiver channel and J(t) is that measured by another receiver channel. Equation (6) is based on equation (1), but with I(t + τ) replaced by J(t + τ). Thus, the distribution of particle radius sizes in a sample may be determined from the autocorrelation function (ACF), as shown in Figure 2c, or from the cross-correlation function (CCF) as in equation (6). The relation between the decay time constant, which depends on the particle size, and scattering angle is shown in Figure 3 for a fluid with a viscosity of 1.98 ×10 -3 Nsm -2 (1.98 cP) and refractive index n=1.33, and in Figure 4 for particles in air with a viscosity of 1.71 ×10 -5 Nsm -2 (0.0171 cP) and a refractive index n=1.0. The variation of time constants as a function of scattering angle shows that backscattering measurements at high scattering angles are preferable to scattering measurements at low scattering angles for the following reasons: 1) A 1° error in scattering angle gives a measurement error of approximately 17.2 % at a scattering angle of 10°, but a measurement error of only approximately 0.2 % at 170°. 2) Assuming that a measurement time of 1000 time constants τ C is required to give a reasonably small standard deviation of about 0.3 %, for a 1 µm particle a measurement time of about 5.8 minutes is required at a scattering angle of 170° , but 11.6 hours at 10°, about 120 times as long. 3) The light intensity scattered in a forward direction increases with particles above about 200 nm in size to about 10 times the backscattered intensity. The probability of getting unacceptably high count rates is therefore greater for forward scattering. Dynamic Light Scattering may also be used to measure the velocities of particles in a fluid. In one example, the laser light is incident on a capillary tube through which fluid flows. Light is reflected from the capillary wall onto the receiver and is also scattered by the fluid. The reflected light then combines with the scattered light at the receiver to produce interference. Flow velocities are derived from the time constant τ C according to the following equation: Δ f = 1/(2πτ C ) The velocity component ν in the direction of observation is given by: ν = c (1-cos() n Δ f f where c = velocity of light n = refractive index  = observation angle relative to the velocity vector, and f = frequency of the laser light. Figure 5 shows the ACF's and velocity profiles measured from a sample with flow rates of 0 to 25 mm/s. For more detailed velocity measurements, the velocity profile and the extent of the measured sample need to be taken into account as indicated above in the discussion of the point-like measurement character of DLS and the advantage of scaling up systems to a multi-channel (n-fold) approach enabling the resolution of phenomena in three dimensions. The interference between the reflected and scattered light (heterodyning) can be used to resolve velocity gradients, but extracting the velocity of each of the n-fold channels simultaneously provides a complete profile acquisition on the basis of homodyning. Pronounced oscillations of ACF's are detected prior to nucleation as a result of currents or velocity gradients within a suspension. Figure 6 shows the oscillating ACF (ACF-33) close to the baseline. The time constant of oscillation is 185 ms giving a frequency of approximately 5 Hz which is used in the equation (8) as Δf, so that the velocity gradient in the volume can be determined. This oscillation is no longer present in a measurement of an ACF taken approximately 10 minutes later (ACF-40). Under ideal DLS, with no particle interaction, light is collected from one coherence area (1 × N) of the scattering volume and is analysed by autocorrelation. Using two detectors receiving scattered light from the same area reduces the effect of noise and imperfections of the detector by cross-correlating the signal from both detectors, and multiple scattering is suppressed if one coherence area is observed. The effect of multiple scattering and deviation from free diffusivity is determined conventionally by an angle-dependent test on a sample having a known distribution of particle sizes, using a goniometer. For example, Figure 7 shows the measured deviation from the theoretical curve (dotted) of a sample with a particle radius of 392 nm, with viscosity of 9.0 × 10 -4 Nsm -2 (0.90 cP). If a measured sample includes fluorescent particles, the sample may be illuminated simultaneously with light of a number of discrete wavelengths and light corresponding to the wavelength of fluorescence is detected as described in the third mode-selecting arrangement described above. The detector includes a polychroic element to separate out the desired wavelength. In lidar remote-sensing applications, a portion of the atmosphere is illuminated and backscattering from that portion is detected. The scattering volume is defined by the intersection between the field of view of the transmitter and the receiver. The scattering volume length (svl) varies as a function of atmospheric height and separation a between the transmitter and the receiver. An example of this variation with a 0.1 mrad full angle receiver aperture is shown in Figure 8. As it is very difficult to observe the same portion of the atmosphere at different viewing angles, lidars are usually restricted to a single backscattering angle, which restricts the number of properties which can be measured. With pulsed illumination, backscattered light can be resolved according to time of flight. The temperature profiles in the atmosphere can be measured from elastic backscattering signals and from Raman effect return signals from nitrogen. The air viscosity can then be determined from the air temperature at the relevant height. In embodiments of the present invention, in order to overcome the problems associated with angle-dependent measurements, measurements are made instead at multiple different wavelengths while the transmitter and the receiver remain in the same position. The deviation of the measured from the theoretical time constants τ C is measured at each of the wavelengths, so as to determine the effect of multiple scattering or the degree of particle interaction which leads to a decrease in free diffusivity. Figure 9 shows the equivalence of time constants τ C and wavelengths as calculated and normalised for two different wavelengths of 266 nm and 690 nm. The numbers shown next to the points of the graph correspond to the equivalent wavelengths, in µm, that would produce the same time constant τ C as the scattering angle shown on the abscissa and the actual wavelength (266 nm or 690 nm), if the equivalent wavelength were to illuminate the sample at an angle of 180° (backscattering). The scattering particle size is 2 nm. Different embodiments of the present invention will now be described for different applications of the present invention. Particle Size and Multiple Scattering Measurement Figure 10 shows an apparatus according to an embodiment of the invention, for the measurement of particle size, particle collision and multiple scattering effects in a sample. There are a plurality n of transmitters 10 and corresponding receivers 12. By means of appropriate optical elements (OE), as shown in Figures 1a and 1b, the scattering volumes of each transmitter/receiver channel may be superimposed so that the measurements of each channel are made on the same scattering volume, or the scattering volumes of the channels may be arranged in an array, so as to provide two or three dimensional resolution of the measured properties. Each transmitter 10 comprises one or more lasers 11 each coupled to a respective mono-mode fibre 4, for transmitting a single wavelength λ 1 , λ 2 , λ 3 ... Each receiver 12 comprises a plurality of mono-mode fibres, each having a corresponding detector and discriminator, as in arrangement III described with reference to Figure 1c. The mono-mode fibres may be selected to transmit respective different wavelengths or some of the fibres may transmit the same wavelength. The detectors for each receiver 12 are shown collectively as reference 14 in Figure 10, while the discriminators are shown collectively by reference 16 and the pulse formers by reference 18. For each channel 1 to n, there is a correlator 20a to 20n which auto-correlates the signal from one detector, or cross-correlates the signals from two or more different detectors of light of the same wavelength. The outputs of the cross-correlators, in the form of cross-correlation functions CCF as derived in equation (6) above are input to a processor 22. In a typical implementation, the discriminators 16, pulse formers 18 and correlators 20 are provided by a dedicated card installed in a general purpose computer which is represented by the processor 30. The apparatus shown in Figure 10 allows measurements of the time constant τ C at multiple wavelengths simultaneously from the same sample. In the example shown in Figure 11, measurements of the time constant τ C for particles of radius 15 nm are shown for wavelengths of 632.8 nm and 689.5 nm, with their equivalent scattering angles normalised and calculated for a beam of 618 nm impinging at 180°. Four measurements are made for each wavelength and the maximum, minimum and average measurements being shown. The large systematic error for measurements at 632.8 nm are due to temperature control instabilities in the laser diode used for that wavelength, leading to limited coherence lengths or spectral purities. The measurements are shown in Figure 11 on a graph of time constant τ C against scattering angle to show that the equivalencies exist and the sample has the ideal properties within the standard deviation for these tests. The equivalence is calculated by combining equations (3) and (4) and normalising this to an assumed wavelength impinging at 180°. The array arrangement allows for particle sizing as a function of extended length in a cuvette or a micro-vial plate. To ensure test reproducibility and the homogeneity of the sample, control of sample properties is provided by the quasi-simultaneous sampling of one scattering volume. Figure 10a shows an example of an implementation in which multiple sampling volumes are defined in the lens LE of an eye. The front of the lens LE is positioned in upper part of Figure 10a, and the scattering apparatus is positioned in front of the lens, so as to measure scattering particle size at different parts of the lens in situ. A first approximately parallel beam B1 and a second approximately parallel beam B2 generated by first and second transmitters 10a, 10b (not shown in this Figure) are incident along the lens axis and approximately parallel to the lens surface respectively. Each beam may contain one or more discrete wavelengths generated by different lasers. For each beam B1, B2 there are five receivers 12a to 12e which receive scattered light from respective reception volumes RV11-RV15, RV21-RV25 which intersect the beams B1, B2 to define scattering volumes SV11-SV15, SV21-SV25 respectively for measuring particle size at different parts of the lens LE. A similar arrangement of intersecting beams B and detected volumes V can be used to measure particle size across a micro-vial plate. Depolarised Cross-correlation DLS While the translational diffusion velocity is proportional to the particle radius r , the rotational diffusion velocity is proportional to r 3 and is therefore more sensitive to small changes in size. Moreover, by comparison of rotational and translational diffusion velocities, information on particle shape can be obtained. Rotational velocity can be derived from the depolarised component of scattered light. Figure 10b shows apparatus for depolarised cross-correlation DLS measurement. A laser 11 generates a beam having vertical polarisation which illuminates a sampling volume SV and the scattered light from the sampling volume SV is received by a receiver 12 having at least two mono-mode fibres 2 through respective parallel horizontal polarizing filters Fa, Fb. As in the arrangement III described with reference to Figures 1a and 1b, each fibre 2 is coupled to a detector 14, discriminator 16 and pulse former 18, and the different outputs are correlated by a cross-correlator 20, as described with reference to Figure 10. The complementary adjustment of the polarization direction of the laser 11 compared to the polarizing filters F removes system noise by cross-correlating the two signals. Thus, very subtle changes in signal strength and in the count rate of the depolarized component can be resolved. The depolarized component is used to calculate the time constant of rotational diffusion of the particles, which is proportional to r 3 whereas the transitional diffusion time constant is proportional to r 2 , as shown by equations (3) and (5). Active Remote Sensing Figures 12a and 12b show two alternative embodiments for polychromatic active remote sensing of the atmosphere. These embodiments represent opposite extremes of design, in which the distance between transmitter and receiver is respectively either very large and the scattering volume very short as in Figure 12a, or the distance is not very large and the scattering volume very long as in Figure 12b. For both embodiments, the signal processing apparatus is as shown in Figures 13a to f. In both embodiments, pulsed lasers 11 of respective different wavelengths are used to generate the illuminating beam. The lasers 11 are either triggered by a pulse generator 22, or synchronised with the sampling timing of the correlator 20. In the embodiment shown in Figure 12b, a time-of-flight register 24 measures the time between the pulsing of the laser 11 and the received scattered pulse and thereby calculates the time of flight. The received pulse intensity from each of the detectors is input to a multi-channel integrator 27 which combines the intensity as received from each detector 14 and passes this to a store 28 which stores the intensity as a function of time. This function is passed to the processor 30 which performs evaluation of the lidar data in the conventional manner, to calculate scattering intensity as a function of height. Both the embodiments have a receiver comprising a plurality of mono-mode fibres 4, as in the arrangement III shown in Figures 1a and 1b, and detectors 14. Figure 13a shows the waveform of one of the pulsed lasers 11, while Figure 13b shows the received intensity signal as measured by one of the detectors 14. The received signal is processed by removing the intervals between the received signal pulses to give a concatenated signal as shown in Figure 13c. The concatenated signal is then digitally high-pass filtered by filter 26 to give the waveform shown in Figure 13d. The filtering removes the effect of fluctuation of the number of scattering particles in the scattering volume. The resultant signal is then cross-correlated to give the CCF shown in Figure 13e, from which the radius distribution as shown in Figure 13f is derived by the method described above. In the embodiment shown in Figure 12b, the auto-correlation function is calculated as a function of atmospheric height. In both Figures 12a and 12b, the optical equipment OE comprises a telescope defining the transmitted beam B and a separate telescope defining the receiving volume RV. In one example, one of the lasers is a continuous wave 1 W NdYAG laser emitting light at 532 nm is used to illuminate a scattering volume as given in Figure 8 at a height of 10000 m and a transmitter-receiver distance of 1000 m. Applications The above embodiments, particularly the configurations of the electronics and the receiver fibres, can be used for many different applications, as described below. The apparatus shown in Figure 8 may be used both for laboratory-based studies, such as aerosol or combustion study, and for reduced gravity environment studies, such as may be carried out on spacecraft. As the apparatus can be miniaturized, it is particularly suitable for a spacecraft payload. Possible microgravity applications include spectroscopy on complex chemical systems, fibre protein aggregations, kinetic studies and nucleation onset together with the study of early crystal growth. For a given payload size, the apparatus allows a higher density of devices; this is particularly important for a Protein Crystallization Diagnostics Facility. The apparatus could also be applied to most types of ground-based lidar for the determination of multiple scattering or sizing of particles, or performance of PCS on droplets or scattering centres, or determination of visco-elastic properties. The apparatus is particularly useful for applications in ophthalmology, where it is difficult to obtain angle-dependent measurements under in vivo conditions. For example, the mechanisms leading to opacification of the eye lens and the changes in the vitreous body as a function of age may be measured. The lidars shown in Figures 12a and 12b could be installed on two satellites, with the transmitter on one, and the receiver on the other or the transmitter could be located on a satellite with the receiver on the ground. In the atmospheric sensing application, the effect of multiple scattering, the limitation on diffusion by particle collision and the velocity components of particles in clouds can be differentiated. This could be used for example to determine the horizontal displacement of a cloud of pollutants in atmospheric boundary layers, by determining the vertical profile of velocity vectors. The embodiment can be used to measure properties of aerosols in the atmosphere, such as particle or droplet size, and tracking of pollution especially in regions of high aerosol loading. The determination of particle sizes in the troposphere is particularly important, as ozone depletion is thought to be linked to the nature of particles in tropospheric clouds ('Unexpected Low Ozone Concentration in Midlatitude Tropospheric Ice Clouds: A Case Study', J. Reichardt et. al., Geophysical Research Letters, Vol. 23, No. 15, pp. 1929-1932, 15 July 1996). The active remote sensing embodiment of the present invention may be used to determine in situ the pressure and temperature dependency of condensation of atmospheric water vapour. The measurement of smallest particle sizes as a function of atmospheric height may be used to estimate water vapour cross-sections, for which the published data are in considerable disagreement. The application to atmospheric measurements allows the measurement of the wavelength dependency of scattering by ice-crystals or droplets in optically thick layers. Dutch patent no. 1001369 describes how DLS or PCS test results may be applied to data derived from low concentrations of particles suspended in air as aerosols. Embodiments of the present invention allow a correction factor to be determined between multiple scattering, the scattering coefficient and the diameter adjustment of the field of view. The correction factor is determined by calculating the separation between the theoretical curve and the measured curve for the sample under observation, as shown for example in Figure 7. For example, if multiple scattering occurs, the derived ACF will not be plotted as a function of intercept and sampling time structure of the correlator, but as a function of the square root of time. This is called diffusive wave spectroscopy, as described in Chapter 16, 'Diffusing Wave Spectroscopy' by David A. Weitz and David J. Pine of the book 'Dynamic Light Scattering - The Method and Some Applications' edited by Wyn Brown, Clarendon Press 1993, and could be applied to the above embodiments in real time.

Claims: A dynamic light scattering method, comprising the following steps: illuminating a fluid scattering volume (SV) with temporally and spatially coherent radiation (B) at a plurality of discrete wavelengths, detecting an intensity fluctuation of radiation scattered by scattering centers contained in said fluid scattering volume (SV) for each of said plurality of discrete wavelengths, and determining a property of said fluid scattering volume from said detected intensity fluctuation of radiation wherein it further comprises the following step : detecting said intensity fluctuation of radiation at a predetermined constant detection angle for each of said wavelengths, and : characterized in that said step of determining a property comprises the step of determining a plurality of equivalent scattering angles related to said scattering centers from said detected intensity fluctuation of radiation of said plurality of discrete wavelengths. Method according to claim 1, further comprising the steps of detecting an electronic signal resulting from said intensity fluctuation of scattered radiation and correlating the detected signal for each of said wavelengths as a function of time so as to derive said property. Method according to claim 1 or 2, wherein the fluid scattering volume (SV) is simultaneously illuminated with each of said plurality of discrete wavelengths. Method according to any of the preceding claim, wherein each of said wavelengths is detected at approximately the same scattering angle. Method according to any of the preceding claims, wherein the scattered radiation is received by an optical arrangement including a mode-selective and mono-mode-transmission medium (2, 12a, 12n) prior to detection. Method according to any preceding claims, wherein an impinging beam is transmitted by an optical arrangement including a mode-selective and mono-mode- transmission medium (4, 10a, 10b, 10n) prior to transmission. Method according to any preceding claims, wherein the scattered radiation is received by a plurality of mono-mode optical fibres (2) each coupled to a respective receiver (14) for detecting the radiation received by said fibre (2). Method according to any preceding claims, wherein the intensity fluctuation of the scattered radiation of one or more of said wavelengths is received by at least one pair of said receivers (14) corresponding to each of said one or more wavelengths and the outputs of each said pair are cross-correlated. Method according to any preceding claims, wherein at least one further fluid scattering volume (SV) is illuminated with said radiation (B), and the intensity fluctuation of said radiation scattered from said at least one further fluid scattering volume (SV) is detected. Method according to any of the preceding claims further comprising the following steps: detecting the velocity of scattering centers in a fluid including the steps ofilluminating the fluid scattering volume (SV) with temporally coherent radiation (B) of a discrete wavelength, detecting the fluctuation of radiation scattered by said fluid volume (SV). correlating the detected intensity as a function of time to generate acorrelation function (ACF-33); determining one or more oscillation frequencies of components of said correlation function (ACF-33); and deriving velocity components of scattering particles in said fluid volume (SV) from said one or more oscillation frequencies. Method according to claim 10, wherein the velocity components are derived by means of the relation ν = c Δ f / nf where ν  is the velocity component is the direction of detection, Δf is the oscillation frequency, f is the frequency of the radiation, n is the refractive index of the fluid volume, and c is the speed of light. Method according to any of the preceding claims, further comprising the steps of: illuminating a fluid scattering volume (SV) constituted by the earth's atmosphere with pulses of temporally and spatially coherent radiation, detecting the fluctuation of said radiation pulses scattered by said fluid scattering volume (SV), and correlating said detected fluctuation as a function of time so as to derive a property of said fluid scattering volume (SV) of the earth's atmosphere. Method according to claim 12, further comprising the step of: concatenating the detected pulses to generate a fluctuation function; wherein said correlating step comprises correlating said fluctuation function. Method according to claim 12 or 13, further comprising the step of: high-pass filtering said fluctuation function prior to correlation. Dynamic light scattering apparatus, comprising: a radiation source (11) for illuminating a fluid scattering volume (SV) with temporally and spatially coherent radiation (B) at a plurality of discrete wavelengths, a detector (14) for detecting the fluctuation of said radiation scattered by said fluid volume (SV)for each of said plurality of discrete wavelengths, and means for determining a property of said fluid scattering volume from said detected fluctuation of radiation wherein said detector (14) is arranged to detect the intensity fluctuation of radiation for each of said wavelengths at a predetermined constant detection angle, and characterized in that said means for determining a property of said fluid scattering volume comprises means for determining a plurality of equivalent scattering angles from said fluctuation of radiation of said plurality of discrete wavelengths. Apparatus according to claim 15, further comprising a plurality of lasers and a plurality of detectors. Apparatus according to claim 15 or 16 wherein said detector (14) comprises a plurality of mono-mode fibres (2) each arranged to detect said radiation at approximately the same scattering angle and each coupled to a respective radiation detecting element. Apparatus according to any of claims 15 to 17, wherein said detector comprises mono mode and/or few mode fibres. Apparatus according to any of claims 15 to 18, wherein each of said one or more lasers (11) is coupled to a respective one or more mono-mode fibre (4), respectively. Apparatus according to any of claims 15 to 19, wherein the radiation source (11) is arranged so that it transmits said radiation through, and the detector (14) is arranged so that it receives said scattered radiation through, respective different optical elements (OE). Apparatus according to any of claims 15 to 20, further including a cross-correlator (20) for cross-correlating as a function of time the outputs of at least two of said radiation detecting elements (14) receiving said scattered radiation at the same wavelength. Apparatus according to any one of claims 15 to 21, wherein said illuminating radiation (B) is polarized, and the detector (14) includes polarizing means (Fa, Fb) so as to detect a depolarized component of said scattered radiation. Apparatus according to any one of claims 15 to 22, comprising : a radiation source (11) for illuminating a scattering volume (SV) of the atmosphere with pulses of temporally and spatially coherent radiation, a detector (14) for detecting the fluctuation of said radiation pulses scattered by said fluid volume (SV), and a correlator (20) for correlating said detected fluctuation as a function of time so as to derive said property of said scattering volume (SV).

Inventors: DIECKMANN MATTHIAS WILHELM MAR, DIECKMANN, MATTHIAS WILHELM MARIA
Applicants: EUROP AGENCE SPATIALE, EUROPEAN SPACE AGENCY
Application year: 1999
Main classification: G01J3
IPC codes: G01J3