Extending the range of turbidity measurement using polarimetry
09823183 ยท 2017-11-21
Assignee
Inventors
Cpc classification
G01N2021/4769
PHYSICS
International classification
Abstract
Turbidity measurements are obtained by directing a polarized optical beam to a scattering sample. Scattered portions of the beam are measured in orthogonal polarization states to determine a scattering minimum and a scattering maximum. These values are used to determine a degree of polarization of the scattered portions of the beam, and concentrations of scattering materials or turbidity can be estimated using the degree of polarization. Typically, linear polarizations are used, and scattering is measured along an axis that orthogonal to the direction of propagation of the polarized optical beam.
Claims
1. An apparatus, comprising: a light source situated to direct a light flux in a first state of polarization (SOP) along first axis; a sample container situated on the first axis so as to receive the light flux; at least one detector situated along a second axis that is at an angle with respect to the first axis, the detector producing an electrical detection signal in response to a portion of the light flux scattered in the sample container; at least one polarizer situated along the second axis between the sample container and the detector, wherein the at least one polarizer is oriented so as to produce at least a first electrical detection signal and a second electrical detection signal simultaneous from the photodetector corresponding to different states of polarization of the portion of the scattered light flux; and a processor that produces an estimate of specimen scattering based on the first and second electrical detection signals.
2. The apparatus of claim 1, wherein the first SOP is a linear SOP and the second axis is perpendicular to a polarization direction associated with the linear SOP.
3. The apparatus of claim 2, wherein the first SOP is a linear SOP and the second axis is parallel to a polarization direction associated with the linear SOP.
4. The apparatus of claim 1, wherein the first SOP is a circular SOP.
5. The apparatus of claim 1, further comprising a rotational stage coupled to rotate the at least one polarizer so as to produce the first and second electrical detection signals.
6. The apparatus of claim 1, wherein the first and second electrical detection signals correspond to a maximum and a minimum, respectively.
7. The apparatus of claim 1, wherein the at least one polarizer includes a first polarizer portion and a second polarizer portion that transmit different SOPs, and the at least one detector comprises a first detector and a second detector situated to receive respective light flux portions from the first and second polarizer portions so as to produce the first and second electrical detection signals.
8. The apparatus of claim 7, wherein the first detector and second detector are detector segments.
9. An apparatus, comprising: a light source situated to direct a light flux in a first state of polarization (SOP) along first axis; a sample container situated on the first axis so as to receive the light flux; at least one detector situated along a second axis that is at an angle with respect to the first axis, the detector producing an electrical detection signal in response to a portion of the light flux scattered in the sample container; at least one polarizer situated along the second axis between the sample container and the detector, wherein the at least one polarizer is oriented so as to produce at least a first electrical detection signal and a second electrical detection signal from the photodetector corresponding to different states of polarization of the portion of the scattered light flux; and a processor that produces an estimate of specimen scattering based on the first and second electrical detection signals, wherein the first SOP is a linear SOP, the at least one polarizer is a linear polarizer, the second axis is perpendicular to the first axis, the at least one polarizer is oriented so that the first electrical detection signal and the second electrical detection signal correspond to a maximum and a minimum, and the processor produces the estimate of specimen scattering based on a degree of linear polarization based on the first and second electrical detection signals.
10. The apparatus of claim 1, wherein the light source includes a polarizer situated so that the light flux has the first state of polarization.
11. The apparatus of claim 10, wherein the polarizer is rotatable so as to select the first state of polarization.
Description
BRIEF DESCRIPTION OF DRAWINGS
(1)
(2)
(3)
(4)
(5)
DETAILED DESCRIPTION
(6) As used in this application and in the claims, the singular forms a, an, and the include the plural forms unless the context clearly dictates otherwise. Additionally, the term includes means comprises. Further, the term coupled does not exclude the presence of intermediate elements between the coupled items.
(7) The systems, apparatus, and methods described herein should not be construed as limiting in any way. Instead, the present disclosure is directed toward all novel and non-obvious features and aspects of the various disclosed embodiments, alone and in various combinations and sub-combinations with one another. The disclosed systems, methods, and apparatus are not limited to any specific aspect or feature or combinations thereof, nor do the disclosed systems, methods, and apparatus require that any one or more specific advantages be present or problems be solved. Any theories of operation are to facilitate explanation, but the disclosed systems, methods, and apparatus are not limited to such theories of operation.
(8) Although the operations of some of the disclosed methods are described in a particular, sequential order for convenient presentation, it should be understood that this manner of description encompasses rearrangement, unless a particular ordering is required by specific language set forth below. For example, operations described sequentially may in some cases be rearranged or performed concurrently. Moreover, for the sake of simplicity, the attached figures may not show the various ways in which the disclosed systems, methods, and apparatus can be used in conjunction with other systems, methods, and apparatus. Additionally, the description sometimes uses terms like produce and provide to describe the disclosed methods. These terms are high-level abstractions of the actual operations that are performed. The actual operations that correspond to these terms will vary depending on the particular implementation and are readily discernible by one of ordinary skill in the art.
(9) In some examples, values, procedures, or apparatus' are referred to as lowest, best, minimum, or the like. It will be appreciated that such descriptions are intended to indicate that a selection among many used functional alternatives can be made, and such selections need not be better, smaller, or otherwise preferable to other selections. In some cases, minimum refers to a minimum value or to values that are within 10%, 20%, or 25% of the minimum value. In some cases, maximum refers to a maximum value or to values that are within 10%, 20%, or 25% of the maximum value.
(10) Examples are described with reference to directions indicated as above, below, upper, lower, and the like. These terms are used for convenient description, but do not imply any particular spatial orientation.
(11) For convenience, beams and light fluxes are described as propagating along one or more axes. Such axes generally are based on one or more line segments so that an axis can include a number of non-collinear segments as the axis is bent or folded or otherwise responsive to mirrors, prisms, lenses, and other optical elements. The term lens is used herein to refer to a single refractive optical element (a singlet) or a compound lens that includes one or more singlets, doublets, or other compound lenses. In some examples, beams are shaped or directed by refractive optical elements, but in other examples, reflective optical elements such as mirrors are used, or combinations of refractive and reflective elements are used. Such optical systems can be referred to as dioptric, catoptric, and catadioptric, respectively. Other types of refractive, reflective, diffractive, holographic and other optical elements can be used as may be convenient.
(12) As used herein, a polarized light flux generally has 90% or more (usually much more) of the associated optical power in the selected state of polarization. Angles are generally intended to be specified with about 10 degrees, 5 degrees, 1, degree, or 0.5 degrees and axes referred to as orthogonal can have similar deviations from 90 degrees. In some cases, a final or other polarizer is referred to as an analyzer, if convenient.
(13) As used herein, optical radiation refers to electromagnetic radiation at wavelengths of between about 100 nm and 10 m, and typically between about 200 nm and 2 m. In some examples, optical radiation is referred to as one or more beams. For convenience, optical radiation is referred to as light or a light flux in some examples, but need not be at visible wavelengths.
(14) As used herein, optical detectors (generally referred to as detectors herein) produce electrical signals corresponding to an incident light flux. Such electrical signals can be represented as voltages, currents, and as either analog or digital signals. Measurement of an optical signal generally refers to detection of at least a portion of an optical beam with a detector.
(15) Various light sources can be used in the disclosed methods and apparatus such as lasers, light emitting diodes, arc lamps and other lamps. In some cases, such sources produce polarized light fluxes and in such cases, additional polarizers may not be necessary to provide a polarized light flux for specimen measurement. Polarizer orientation can be adjusted by rotation. A waveplate can be used to produce an effective rotation, without rotation of the polarizer. Electrically switchable retarders or other devices such as liquid crystal, Faraday rotator, photoelastic modulator (PEM), or Pockels cells can be used, but are generally omitted from the following description.
(16) Referring to
(17) The light source 102 can be a broadband light source such as a lamp, but other light sources such as lasers, laser diodes, and light emitting diodes can be used. Optical power from the light source 102 can be monitored with an additional detector, and some light sources permit direct modulation so that synchronous detection can be used even with the optical chopper 112 omitted. The relative placement of the polarizer 110 and the optical chopper 112 shown in
(18) A right handed coordinate systems 151 is shown having a +Y-axis that extends out of the plane
(19) The detection system 128 can include one or more lock-in amplifiers, transimpedance amplifiers or other circuitry for processing photosignals from some or all detectors.
(20) Data storage, processing, and control of detector bias, rotational stage positions, chopper operation (such as frequency), control of light source output and amplifier gain can also be controlled with a processor using computer-executable instructions stored in a memory such as RAM or a hard drive. Processors can be implemented as stand-alone or dedicated processors, or provided with a general purpose computer, a laptop, a hand held device, or a tablet computer.
(21) Various detectors such as photodiodes, avalanche photodiodes, photomultipliers or others are suitable, but for ease of description, operation based on photodiodes and the associated photocurrents is described. The detector 116 produces a photocurrent I.sub.PD00 associated with detection of a ballistic portion of the input polarized light flux that is transmitted by the specimen 115; the detector 120 produces a photocurrent I.sub.PD+90 associated with detection of a polarized portion of the input polarized light flux as scattered by the specimen 115; the detector 124 produces a photocurrent I.sub.PD90 associated with detection of an unpolarized portion of the input polarized light flux as scattered by the specimen. A photocurrent associated with the detector 126 is noted as I.sub.PDs. The numerous detectors included in the system 100 of
(22) I.sub.PD00, I.sub.PD+90, I.sub.PD90, and I.sub.PDs are used above to refer to photocurrents. However, this notation is used herein to refer to light flux intensities or corresponding electrical signals at or produced by detectors which can be indicated as PD.sub.00, PD.sub.+90, PD.sub.90, and PD.sub.s for convenience.
Analytical Background
(23) In the following, a representative analytical approach to turbidity determination is provided based on a system such as that of
I.sub.0=.Math.I.sub.PD.sub.
wherein .Math.I.sub.PD.sub.
.Math.I.sub.PD.sub.
wherein .sub.t is an optical density corresponding to loss of the original beam intensity due to both absorption and scatter, i.e., .sub.t=.sub.a+.sub.s, respectively. From Eq. (2), the diffuse component therefore consists of all the light scattered out of the original input beam that is not absorbed; this can be written as:
.sub.1,2.Math.I.sub.PD.sub.
When polarization is ignored I.sub.PD.sub.
.sub.1.Math.I.sub.PD.sub.
(24) In terms of Malus' Law, Eq. (4) can be rewritten as:
I.sub.t()=.Math.{[1+P.Math.cos(2)]|.sub.max+[1+P.Math.cos(2)]|.sub.min},(5)
wherein P is a degree of polarization and is an azimuthal angle of analyzer orientation with respect to the plane of incidence (i.e., an input SOP) and min and max reference the associated extrema. Note that for an unpolarized light input, in the absence of polarization effects due to scatter, P=0 and therefore I.sub.t lacks a polarization component and is strictly made up of depolarized light, which does not obey Malus' Law. Likewise in the absence of an analyzer in front of a polarization insensitive detector (i.e., for PD.sub.90),
[1+P.Math.cos(2)]|.sub.max[1+P.Math.cos(2)]|.sub.min=k,(6)
wherein k is a constant conveying the ensuing polarization insensitivity. From Eq. (4), for normalized intensity,
I.sub.t=I.sub.p+I.sub.d=1.(7)
[1+P.Math.cos(2)]|.sub.max and [1+P.Math.cos(2)]|.sub.min occur for |.sub.0 and |.sub.90, respectively, from Eq. (5). Therefore,
I.sub.t=I(|.sub.0)+I(|.sub.90).(8)
Rearranging Eq. (8) to directly map to Eq. (5), yields:
I.sub.t=I(|.sub.0)I(|.sub.90)+2.Math.I(|.sub.90)I.sub.p+I.sub.d,(9)
wherein I.sub.p=I(|.sub.0)I(|.sub.90) and I.sub.d=2.Math.I(|.sub.90). Also, note that
(25)
(26) The fraction of the detected intensity captured by PD.sub.+90 that is polarized is the degree of linear polarization (DOLP) and is given by:
(27)
However, other functions of the difference between signals associated with different SOPs in the scattered beam can be used. From Eq. (11) it is clear therefore that DOLP is a representative metric for input SOP decorrelation.
Example Implementation
(28) Representative measurement results are presented for DOLP based on polarization-analyzed detection with PD.sub.+90 and non-polarized (NP) detection by PD.sub.00 and PD.sub.90 as functions of non-dairy creamer (NDC) concentration, as shown in
(29) For the V-SOP input, note that the ballistic NP-PD.sub.00 trend closely resembles that of the polarized PD.sub.+90 detection. Conversely, the NP-PD.sub.90 detection mode shows an initial increase with increased NDC scatter concentration and then saturates and proceeds to gradually decrease with increasing NDC concentration. This non-monotonic behavior of NP side scatter detection is anticipated by Eqs. (3) and (4), where:
.sub.2.Math.I.sub.PD.sub.
and as sample concentration increases, .sub.t.fwdarw., thus leading to e.sup..sup.
.sub.2.Math.I.sub.PD.sub.
and if .sub.a<<1, then expanding the exponential function yields, to 1.sup.st order approximation,
.sub.2.Math.I.sub.PD.sub.
Clearly, past the initial saturation point, if .sub.a were to increase with NDC concentration due to greater scatterer agglomeration (a typical consequence of increased scatterer density), then the value of I.sub.sat would decrease also as seen in
(30) The maximal verifiable turbidity measurable in this example was that for the 2.067 mg/mL concentration (784 NTUs). The NP side scatter detection mode was able to quantitate up to 7.233 mg/mL with individual errors <10% for a mean quantification error of 3.91%. Polarization-based detection was able to quantitate up to 13.433 mg/mL with 3.38% mean quantification error. The NP ballistic detection mode quantitated up to 11.367 mg/mL with 11.41% mean quantification error, but individual errors varied greatly and were much higher at just under 19%.
(31) The exponential model inverted to predict concentration C is:
y=a.sub.1.Math.e.sup.(a.sup.
wherein y represents either DOLP or NP-intensity The model fit parameters are presented in Table 2 for R.sup.2=0.9984, 0.9934, and 0.9994 (polarized-DOLP and non-polarized, NP-PD.sub.00 and -PD.sub.90 implementations respectively). Additionally, the models were used to predict the 2.067 mg/mL concentration validation sample with the following respective errors: 5.47%, 22.91%, and 8.07%.
(32) Polarization-based quantitation does not require determining DOLP. In other examples, signals for the analyzer oriented at 0 and 90 degrees (i.e., corresponding to max and min respectively) would be measured. A range of values associated with max/min values for one or more reference samples of known turbidity would be used for calibration so that measured max/min values for a sample could be associated with sample turbidity or concentration using a look up table or other calibration table or curve. In some cases (as noted previously) the two photodiode measurements corresponding to analyzer at 0 and 90 degrees are measure simultaneously. In other examples, the detector 124 in
(33) As noted in Table 1, * indicates values of uncertain accuracy as other, slightly larger concentrations (not included in Table 1) also returned the same value. A dagger () indicates that the computed mean did not include the value for the concentration 15.500 mg/mL because the associated error was greater than 20%. The values in the second column are in Nephelometric Turbidity Units (NTUs) and were obtained with a conventional turbidity meter. Other indications in Table 1 are: not available due to turbidimeter measurement error message (NA); not calculated (NC) as either outside of a logarithmic range for NP-PD.sub.00 or that concentration was outside the range that was modeled for NP-PD90; and predicted (Pred).
(34) TABLE-US-00001 TABLE 1 Measurement results. Concentration P (V-SOP) NP-PD.sub.00 NP-PD.sub.-90 (mg/mL) NTUs Pred. % Error Pred. % Error Pred. % Error 1.033 383 1.006 2.63 1.23 18.81 1.01 2.29 2.067 784 3.100 1000* 3.241 4.54 3.16 1.99 3.08 0.69 5.167 NA 4.996 3.30 4.57 11.47 5.59 8.27 7.233 NA 7.375 1.96 6.57 9.16 7.45 3.01 9.300 NA 8.986 3.38 7.59 18.36 7.80 16.08 11.367 NA 11.044 2.84 12.35 8.68 NC NC 13.433 NA 12.757 5.04 NC NC NC NC 15.500 NA 20.525 32.42 NC NC NC NC Mean 3.38 11.41 6.07
(35) TABLE-US-00002 TABLE 2 Model fit parameters. DOLP(V-SOP) PD.sub.00 .sup.PD.sub.90 Coef. Value Value Value a.sub.1 0.3115 0.0082 1.3020 0.054 1.7830 0.0424 a.sub.2 0.2525 0.0165 0.3732 0.042 0.4734 0.0328 a.sub.3 0.0879 0.0045 1.0050 0.030 1.8030 0.0311
(36) In this example, a V-SOP linear polarization input with polarization analysis for detection of the min and max intensity values of side-scatter (i.e., at 90 to the incident propagation direction) was used to measure scattering sample concentration (i.e., turbidity). NDC was used a specimen but other specimens as biological tissues, drinking water, foods, cosmetics, pharmaceuticals and waste water and others can be similarly evaluated.
Additional Examples
(37) Referring to
(38)
(39)
(40) A beam portion scattered along a perpendicular axis 530 is directed to a polarizer assembly 532 that contain polarizer segments 532A, 532B associated with orthogonal SOPs. A detector assembly 536 has detectors 536A, 536B that are situated to receive scattered beam portions from the polarizer segments 532A, 532B, respectively. The detection system 528 can use the electrical signals from the detectors 536A, 536B to produce concentration or turbidity estimates. The detector 536 can be a segmented detector, or two or more individual detectors can be used. The polarizer segments 532A, 532B typically are polarizers having different orientations, and fixed into a common assembly. In some cases, a single polarizer can be used an provide with a wave retarder on an input side, wherein the wave retarder is oriented at 45 degrees with respect to the polarization axis of the single polarizer. In either case, measurements that are effectively at two orthogonal directions are obtained. A lens 540 can be used so that scattered beam portions are directed to each of the detectors 536A, 536B.
(41) As shown in this example, polarization based measurement can use either or both of the axes 508, 530, or an axis or axes at angles between 0 and 180 degrees, or 0 and 90 degrees. The light source can be a broadband source, a source that comprises a plurality spectral lines, and spectral filters can be used to select a spectral region of interest.
(42) In view of the many possible embodiments to which the principles of the disclosed technology may be applied, it should be recognized that the illustrated embodiments are only representative examples and should not be taken as limiting the scope of the disclosure. Alternatives specifically addressed in these sections are merely exemplary and do not constitute all possible alternatives to the embodiments described herein. For instance, various components of systems described herein may be combined in function and use. I claim all that is encompassed by the appended claims.