Translator Disclaimer
1 July 2006 Angular measurements of light scattered by turbid chiral media using linear Stokes polarimeter
Author Affiliations +
The effects of turbid chiral media on light polarization are studied in different directions around the scattering samples using a refined linear Stokes polarimeter, which simplifies the signal analysis, and allows for the detailed investigations of scattered light. Because no moving parts are involved in a measurement at a specific detection direction, the determination accuracy of polarization states is increased. The results show that light depolarization increases with both turbidity and detection angle for low and moderately turbid samples; however, the angular dependence decreases with increasing turbidity. When the turbidity is increased to ~100 cm–1, the depolarization becomes higher in the forward than in the backward direction. Polarization sensitive Monte Carlo simulations are used to verify some experimental observations. The results also demonstrate that surviving linear polarization fraction and overall intensity are more sensitive to the increase of glucose concentration in backward than in the forward direction in highly turbid media, indicating that backward geometry may be preferable for potential glucose detection in a biomedical context. Comparison measurements with optically inactive glycerol suggest that the refractive index matching effect, and not the chiral nature of the solute, dominates the observed optical rotation engendered by glucose in highly turbid media.



Diabetes mellitus is a medical condition in which the body cannot adequately produce or effectively use insulin that is needed to maintain a normal blood glucose level. It is reported that diabetes afflicts over 15 million people in North America and 100 million people worldwide.1 To avoid or to delay the complication of diabetes including blindness, heart disease, and kidney failure, diabetics must regulate their insulin injections and caloric intake by monitoring their blood glucose levels. Most current methods for glucose monitoring require drawing blood (typically via a finger prick), which is painful, inconvenient, potentially infectious, and otherwise suboptimal. Significant efforts have centered on the search for noninvasive glucose monitoring methods. Optical methods are potentially well-suited to tackle this currently unsolved problem due to their noninvasive, nonharmful nature and sensitivity to glucose-induced optical changes. 2, 3, 4, 5, 6, 7, 8, 9, 10

Recently, polarized light has been used to extract quantitative information from the optically thick tissues with which it interacts. 3, 10, 11, 12, 13, 14, 15 Because the optics of skin is dominated by scattering effects in the tissue optical window range of 600 to 1200nm , the understanding of polarization properties of light that has escaped from highly scattering media is of great importance. In addition to the forward direction, some research has examined the detection of scattered light in other directions,11, 12, 13 to investigate the detection geometry where the experimental observables such as optical rotation, degree of polarization, and intensity respond most to changes in glucose concentration. These attempts were challenged by the complexity of light-media interactions and experimental difficulties. 12, 13, 14, 15

The approach presented in this paper employs a refined linear Stokes polarimeter to detect the emerging light polarization states in a variety of detection angles around a scattering sample and examines the sensitivity of the derived changes in optical rotation and surviving linear polarization fraction to changes in the glucose concentration. This polarimeter design is similar to the one presented by Hunt and Huffman in the 1970s and subsequently used by many researchers. 16, 17, 18, 19 The main innovation in the present embodiment lies in the measurement of various harmonics of the signal, which are directly related to the linear Stokes vector components of the sample-scattered light. In other words, the current arrangement with the polarization modulation in the analyzer arm after the sample, permits completely assumption-independent determination of sample-induced light polarization parameters. This is a significant improvement over many previous attempts at Mueller matrix or Stokes vector polarimetry, which had to assume specific forms of the sample matrix or of the scattered light components for signal analysis. These assumptions become questionable as one encounters complicated media, such as biological tissues, that may exhibit several optically active effects at once (linear and circular birefringence or dichroism) while also causing light depolarization due to multiple scattering. Also, the electronic synchronous detection at three different modulation frequencies provides sufficient information to determine the Stokes parameters of interest, without the need for mechanical moving parts (rotating linear polarizers and/or wave plates, swapping in and out of polarizing optics). This reduces noise in the polarization measurements. Utilizing the linear Stokes polarimeter, the relative contributions of glucose concentration, sample turbidity, and detection angle geometry to the polarization properties of the emerging light are systematically investigated. The role of optical activity of glucose, the foundation for glucometry in clear and low turbid media, in the observed optical rotation in highly turbid media is also investigated.


Theory and Signal Analysis



Figure 1 illustrates the simplified linear Stokes polarimeter setup. The approach is centered on polarization modulation accomplished by the photoelastic modulator (PEM) and synchronous detection via a lock-in amplifier, allowing a sensitive means of detecting low levels of polarization signals in the presence of largely depolarized background. Specifically, the laser beam passes through a linear polarizer P (pass axis horizontal, parallel to the optical table) and impinges on the sample. The scattered light from the sample passes through a quarter wave plate, whose fast-axis orientation is 45° with respect to optical table, and goes through the PEM block orientated horizontally. The PEM is a resonant device operating at 50kHz whose user-defined peak retardance δ0 is selected to optimize polarization signal strength (see below). The light then passes through an analyzer A set at 45° with respect to the optical table and finally reaches the detector. Using standard Stokes vector notation for light polarization and Mueller matrices to represent the optical elements that transform these polarization states, we now analyze the system in Fig. 1. Si represents the Stokes vector of the incident beam

Eq. 1

where horizontal notation is used for convenience, with T standing for the transpose of the vector. Se is the Stokes vector of the emerging beam after sample interactions,

Eq. 2

If the elements of the Stokes vector in Eq. 2 are normalized by intensity I , Eq. 2 becomes

Eq. 3

where q=QI , u=UI , and v=VI . Sf is the Stokes vector of the beam reaching the detector

Eq. 4

Relating the detected signal to the sample-emerging beam via Mueller calculus,

Eq. 5

where MWP is the Mueller matrix of the quarter wave plate, MPEM is the Mueller matrix of the PEM, and MA is the Mueller matrix of the analyzer. For the selected orientations, the Mueller matrices are

Eq. 6


Eq. 7


Eq. 8

Putting Eqs. 6, 7, 8, 3 into Eq. 5, we obtain

Eq. 9

Performing the multiplications, the measurable part of the final Stokes vector reaching the detector is the time-dependent intensity expressed as

Eq. 10

If(t)=12[1qsinδ+ucosδ] =I2[1qsin(δ0sinωt)+ucos(δ0sinωt)],
where the time-dependent PEM retardation follows δ(t)=δ0sin(ωt) . Applying Fourier expansion in terms of Bessel functions

Eq. 11


Eq. 12

to sine and cosine terms in Eq. 10, there results

Eq. 13

If the peak retardance of PEM is set to δ0=2.405 radians so that J0(δ0)=0 , then

Eq. 14

The elements of the Stokes vector of the sample-scattered light, I , q , and u [Eq. 3] can be extracted from the dc component, first harmonic component, and second harmonic component of Eq. 14, respectively. Equation 14 is valid for all scattering directions because no angle-dependent assumptions are used in its derivation. Further, this formulation is based on the general form of sample polarization in Eq. 3, thus avoiding the assumption of specific (and generally unknown) forms of the sample Mueller matrix that may simultaneously exhibit optical rotation, depolarization, and linear anisotropy (birefringence).12, 14 As such, the relationship of the derived polarization parameters to the intrinsic properties of the interrogated sample is further strengthened.

Fig. 1

Simplified linear Stokes polarimeter set-up. P, polarizer; A, analyzer; WP, quarter wave plate; PEM, photoelastic modulator; Si , Stokes vector of the incident beam; Se , Stokes vector of the emerging beam; Sf , Stokes vector of the beam reaching detector. All angles are with respect to the optical table.



Signal Analysis

If the dc component of Eq. 14 is measured by an electronic dc voltmeter, and the first and second harmonic components are measured by a lock-in amplifier, the experimentally observed wave form will be

Eq. 15

where VDC is the dc measurement in volts, V1f is lock-in output in volts (rms) at the PEM’s first harmonic, and V2f is lock-in output in volts (rms) at the PEM’s second harmonic. The factor 2 arises in Eq. 15 because the lock-in amplifier reports rms voltage.

By comparing Eq. 14 with Eq. 15,

Eq. 16


Eq. 17


Eq. 18

where k is an instrumental constant, the same for all three equations assuming constant detector frequency response. Equations 16, 17, 18 yield expressions for q and u in terms of the experimentally measured parameters

Eq. 19


Eq. 20

The optical rotation of the emerging beam α and surviving linear polarization fraction βlinear can be calculated accordingly

Eq. 21


Eq. 22

Strictly speaking, α is the orientation of polarization ellipse of the emerging beam, but can be treated as the optical rotation with reference to the orientation of the polarization ellipse of the incident beam expressed in Eq. 1. In order to increase dc measurement accuracy, for example, to eliminate the effects of light intensity fluctuation, it is useful to employ ac methods based on lock-in amplification. This is convenient in the current set-up that already uses lock-in detection for the PEM-based signals and will ensure the cancellation of the experimentally unknown constant factor k that is needed for the correct derivations of Eqs. 16, 17, 18, 19, 20. If we use a mechanical “square wave” chopper with a duty cycle of 50%, the average dc, 1f and 2f signal intensities will all be halved, and the ratios of VACVDC will remain unchanged. However, the lock-in output of dc measurement Vchopper is not precisely the VDC value. The relationship between the chopper signal measured by the lock-in amplifier Vchopper and VDC is

Eq. 23

which has universal validity for square wave choppers with 50% duty cycle.20 In Eqs. 19, 20, replacing VDC with Eq. 23 and using J1(δ0)=0.51911 and J2(δ0)=0.43169 , we get

Eq. 24


Eq. 25

Equations 24, 25 yield the sample Stokes parameters from data acquired only by synchronous high signal-to-noise ratio lock-in amplifier detection, which in turn yield accurate estimates of sample-induced optical rotations α and linear polarization fraction βlinear via Eqs. 21, 22.


Materials and Experiments


Experiment Setup

Figure 2 shows the diagram of linear Stokes polarimeter, demonstrating the specific embodiment of the general set-up of Fig. 1. The 632.8-nm light from a HeNe laser (Research Electro-Optics, LHRR-120M), chopped by a mechanical chopper at a frequency of 500Hz , passes through a polarizer (pass axis horizontal) and an iris of 1mm diameter, and strikes a round glass cuvette of 1cm inner diameter that contains the sample. For the reported experiments, the effective sampling volume is a function of detection geometry and sample scattering properties. Its approximate range varies from a few cubic millimeters (highly scattering, large detection angles) to tens of cubic millimeters. Further quantification of this important but difficult parameter, and a related measure we call the average photon pathlength, will be reported in a forthcoming publication. The incident power on the sample is 1mW . Part of the emerging light enters an analyzing block consisting of a pair of plano-convex lenses; a quarter wave plate (fast axis at 45° ); a PEM (Hinds Instruments IS-90; horizontal retardation axis, oscillating at 50kHz ), a linear analyzer (pass axis at +45° ), and a polarization-insensitive avalanche photodiode detector (Hamamatsu, C5460), which has a flat frequency response range from dc to 1MHz . The analyzing block is installed on a horizontally rotatable platform forming a detection angle θ with respect to the incident light, with the center of the cylindrical cuvette positioned at the rotation axis of the block. The use of the pair of plano-convex lenses increases measurement sensitivity. The first lens picks up emerging light from an area of the turbid sample that is small enough to avoid significant blurring of the polarization states engendered by spatial averaging. The collected light is then converted into a collimated beam by the lens, an important step as the PEM functions optimally with collimated light. The second lens converges the collimated light, focusing it onto the sensing area of the detector. The signal from the analyzing block is sent to a lock-in amplifier (Standford Research Systems, SR830). The reference input of the lock-in amplifier is toggled between the chopper and PEM controllers. Data acquisition is carried out under computer control, obtaining within 2min 10-measurement averages at 500Hz , 50kHz , and 100kHz as Vchopper , V1f , and V2f for further processing. The detection angle scanning is achieved by rotating the analyzing block platform as it isocentrically pivots around the round cuvette. The angles in the experiment vary from 0 to 155° at 5° intervals. Because of the PEM’s size (15×2.5×5.5cm3) , the measurements at angles larger than θ> 155° are hindered due to the obstruction of incident beam by the analyzing block.

Fig. 2

Schematic diagram of linear Stokes polarimeter. C, optical chopper; P, polarizer; θ , detection angle; L1 and L2 , plano-convex lenses; WP, quarter wave plate; PEM, photoelastic modulator; A, analyzer; D, detector.




Four sets of samples were used in the experiments. The first set consists of turbid media with different scattering coefficients. These are suspensions of polystyrene microspheres (refractive index ns=1.59 , density ρ=1.05gcm3 , radius r=0.7μm ), yielding a scattering efficiency of Qsca=3.56 at λ=632.8nm in deionized water.21 The weight fraction of microspheres fw ranges from 5.5×104 to 2.7×103 . Because scattering coefficient μs is equal to the product of the scattering cross-section σsca=Qscaπr2 and the scatterer volume density ρscafy[(43)πr3] ( fv is the volume fraction of the scatterers), μs={fv[(43)πr3]}Qscaπr2=3Qscafwρ04rρ [ ρ0 is the water density (1gcm3) ], the corresponding scattering coefficients vary from 20 to 100cm1 . The second set of samples are turbid media containing D -glucose with the microsphere weight fraction kept constant at 8.25×104 ( μs is ∼27.8 to 30cm1 , depending on glucose levels) and glucose concentration ranging from 20 to 900mM . The third set of samples are highly turbid media containing D -glucose with the microsphere weight fraction kept constant at 2.75×103 ( μs is ∼92.75 to 99.9cm1 , depending on glucose levels) and glucose concentration ranging from 30 to 900mM . The fourth set of samples is highly turbid media containing optically inactive glycerol with the microsphere weight fraction kept constant at 2.75×103 (same scattering coefficient range as set 3). The glycerol concentrations are such that they yield the same solution refractive index as glucose in the third set of samples. This is used to make comparison measurements on the contributions to optical rotation due to chirality and due to refractive-index matching effects, in which the increase of glucose concentration reduces the refractive index difference between the scattering particles and the surrounding medium, resulting in the decrease of the scattering coefficient of the turbid medium. 8, 22, 23, 24


Results and Discussion

Figure 3 shows how optical rotation α , surviving linear polarization βlinear , and intensity I of the light emerging from turbid media (no glucose) change as a function of detection angle θ and turbidity of sample, indicated by scattering coefficient μs . The results of about 32 measurements are summarized in the figure. Somewhat surprisingly, Fig. 3a shows large optical rotations even in the absence of glucose, as suggested by previous initial observations.3, 15, 18 As seen, the magnitude of α is a strong and complex function of both sample turbidity and detection angle. For lower turbidity ( μs=20 to 40cm1 ), α increases with θ from about 1° to about 4° and then decreases again for θ> 130° . For turbid media with μs=50 to 80cm1 , α displays a monotonic increase with θ from 1 to 12° . For highly turbid media ( μs=90 to 100cm1 ), α exhibits oscillations at smaller detection angles and then increases monotonically from ∼2 to 12° when θ> 90° . It is noticed that in these highly turbid samples, α is slightly nonzero even in the forward direction, possibly arising from minor detector misalignment about θ=0° . Clearly then, multiple scattering significantly changes the orientation of the linear polarization plane, thus resulting in an apparent optical rotation in most directions. The magnitude of this change in linear polarization orientation depends strongly on both sample turbidity and detector orientation angle. Larger θ generally means larger change of the polarization plane orientation. Further studies are underway to interpret these trends. Regardless of the mechanism, the data unequivocally shows that optical rotation results in a turbid polarimetry setting must be interpreted with caution, in that α is not only caused by chirality of glucose molecules as is the case for clear-media glucometry.

Fig. 3

Measured angular and turbidity dependence of polarization states of emerging light from turbid media in the absence of glucose (a) Optical rotation; (b) surviving linear polarization fraction; (c) normalized intensity. The orientation of the horizontal axes is different for the three plots to enable better visualization.


Figure 3b shows that βlinear decreases significantly with increasing turbidity and detection angle for low and moderately turbid media. In the forward direction (θ=0°) , βlinear drops from 0.9 at μs=20cm1 to 0.35 at μs=60cm1 . For a turbid medium with μs=20cm1 , βlinear drops from 0.9 at θ=0° to 0.2 at θ=155° . With further increases of turbidity of the sample, βlinear still decreases but becomes less dependent on detection angle as its overall levels drop. The data suggests that for low and moderately turbid media, photons exiting at smaller angles may be scattered fewer times than the photons exiting at larger detection angles, yielding the smaller depolarization effects. With the increase of sample turbidity, photons encounter more scattering events on route to reaching the detector, resulting in greater depolarization. With further increase in sample turbidity, photons exiting in any direction encounter almost the same amount of scattering events, as suggested by the decreasing angular dependence of βlinear .

Figure 3c shows that the intensity I is not strongly dependent on the detection angle θ except for low turbidity media (μs25cm1) , where I drops from 0.9 to 0.1 when θ increases from 0 to 155° . For samples with μs> 40cm1 , the intensity is low and largely independent of detection angle, although exhibiting a minor increase for θ> 90° .

Selected subsets of the above experimental results were compared with polarization-sensitive Monte Carlo (MC) simulations developed in our group,3, 25 a validated model that expands upon the prior work of Kaplan 26 and Jaillon 27 In these simulations, the photons are propagated between scattering events, as determined by pseudorandom sampling of the scattering mean free path similar to conventional (intensity-only) MC models28 of light propagation in tissue; the polarization information is tracked in the form of individual Stokes vectors, which are summed over a large number of tracked photon histories, to yield the experimentally observable macroscopic properties of interest. Optical activity due to the chiral nature of dissolved glucose is handled by imposing the rotation of the plane of linear polarization between the scattering events, using a Mueller matrix for a standard optical rotator.3, 25 The aim of the comparison is to validate the measurement results and predict the trends of polarization states of scattered light at θ> 155° , where measurements are not currently possible. Typical experiment-simulation comparison curves for glucose-free turbid suspensions are shown in Fig. 4 . The MC simulations of the surviving linear polarization fraction as a function of detection angle [for μs=20cm1 , Fig. 4a] and as a function of scattering coefficient [for θ=0° , Fig. 4b] show similar trends to the experimental measurements, although further refinement may be necessary to improve agreement. The simulation in Fig. 4a predicts a significant increase of βlinear beyond 160° , a region currently unconfirmed with experimental measurements. This suggests the potential benefit of backward detection for highly turbid media, where surviving linear polarization fraction and signal intensity are relatively high. System modifications are underway to enable experimental validations at larger detection angles approaching the exact backscattering direction.

Fig. 4

Experimental data (squares) and MC simulations (line) for glucose-free turbid media in a cylindrical cuvette of 1cm diameter. (a) Surviving linear polarization fraction as a function of detection angle for μs=20cm1 ; (b) surviving linear polarization fraction as a function of scattering coefficient at θ=0° .


Figure 5 shows the effects of increasing glucose (20 to 900mM ) in moderately turbid media (μs=30cm1 in glucose-free suspension) on polarization states of the emerging light detected at five different detection angles of θ=0 , 45, 90, 135, and 155° . Figure 5a shows that optical rotation decreases with glucose concentration at all detection angles except in the forward direction. The observed decrease in α may indicate that the optical rotation due to scattering and the optical rotation due to glucose chirality are in the opposite directions. If former is greater than the latter, addition of glucose can cause a decrease in the net apparent optical rotation. The optical rotation change detected at θ=155° is most sensitive to glucose level change, yielding Δα4.5° for examined concentration change. Selecting a slightly different backward detection angle of θ=135° yields a Δα2.0° for the same glucose range. For comparison, the conventional detection geometry in the forward direction results in a much smaller Δα of +0.8° . The data trends suggest that in the forward direction, glucose chirality is the main contributor to the observed optical activity; in other detection directions, the measured optical activity is dominated by scattering effects (as modulated by the addition of glucose, which (1) reduces scattering (2) is optically active). It is expected that there exists a detection angle at which the measured optical rotation has equal contributions from these reactive index matching and optical activity effects (as borne by the experimental results summarized in Fig. 7 below).

Fig. 5

Effects of glucose concentration (20 to 900mM ) in moderately turbid media ( μs=30cm1 at 0M glucose) on polarization states of emerging light detected at different detection geometries. (a) Optical rotation; (b) surviving linear polarization fraction; (c) normalized intensity. The symbols are data points; the lines are a guide for the eye.


Fig. 7

Glucose effects on optical rotation in turbid media, in order to examine the chirality and refractive index matching contributions. (a) Moderately turbid media (μs30cm1) ; (b) highly turbid media (μs100cm1) . The symbols are data points; the lines are a guide for the eye.


Figure 5b shows the change of corresponding surviving linear polarization fraction due to glucose for the same samples as in Fig. 5a. βlinear increases with glucose concentration at all detection angles with Δβ+0.15 in forward direction and Δβ+0.1 at θ=155° . This change is most likely due to the glucose refractive index matching, which effectively reduces the sample turbidity, thus decreasing the depolarizing effects of multiple scattering. Even through the backward detection displays a slightly reduced Δβlinear and overall lower levels of surviving linear polarization, it is still an attractive detection geometry if both α and βlinear are used as metrics of glucose concentration.

Figure 5c shows the change of normalized intensity I in the presence of varying glucose amounts. As seen, the intensity is higher in the forward hemisphere, and glucose-induced changes are very much directions-dependent, being negligible for θ> 90° and most pronounced in the forward direction. These polarization-independent changes in detected light intensity with varying glucose levels most likely originate from the index-matching effect of glucose, which reduces sample scattering properties as previously outlined.

Figure 6 is analogous to Fig. 5 but displays the results for the case of higher scattering ( μs=100cm1 in glucose-free suspension). Opposite to Fig. 5a, the optical rotation α in Fig. 6a decreases with glucose concentration at a forward detection direction and increases with glucose concentration at larger detection angles (135 and 155° ). The monotonically decreasing α in the forward direction indicates an overall dominance of the glucose refractive index matching effect that reduces scattering and the apparent optical rotation it causes in this high-scattering regime. The increase in α at larger angles implies that more complicated mechanisms are contributing to optical rotation. For example, the reduction in scattering due to glucose-induced refractive index matching effect can cause the observed increase in α at large scattering angles, by lengthening the effective pathlength over which the detected photons can accrue additional optical rotation. Figure 6b shows that βlinear levels and their sensitivity to glucose changes are different in their dependence on θ as compared with optical rotation. The magnitude and increase of βlinear caused by increasing glucose levels are both largest in backward directions. The intensity dependence of Fig. 6c shows that more photons escape the highly turbid medium in the backward hemisphere and the glucose-induced intensity reduction is the largest at θ=155° .

Fig. 6

Effects of glucose concentration (30 to 900mM ) in highly turbid media ( μs=100cm1 at 0M glucose) on polarization states of emerging light detected at different detection geometries. (a) Optical rotation; (b) surviving linear polarization fraction; (c) normalized intensity. The symbols are data points; the lines are a guide for the eye.


Figure 7a displays the relative importance of scattering-induced and optical activity-induced glucose contributions to the observed optical rotation in moderately turbid samples. As θ is increased from 5 to 30° in samples with μs30cm1 , two distinct regimes of optical rotation contribution are evident. Near the forward direction ( θ=5 and 10° ), increasing the glucose concentration leads to nearly monotonic increase in measured optical rotation. Conversely, for detection angles of 25 and 30° , a decrease in optical rotation is seen. We interpret the former behavior as being dominated by optical activity, whereas the latter is indicative of scattering effect dominance as engendered by glucose refractive index matching. The θ=15° detection geometry represents an intermediate case of roughly equal contributions, at least for this sample geometry and this scattering coefficient (μs30cm1) . An analogous graph for higher scattering (μs100cm1) seen in Fig. 7b shows that the transition region occurs at much higher detection angles and also indicates overall higher optical rotation levels. It appears that α detected at θ=125° has roughly equal contributions from chirality and from index matching for the examined range of glucose concentration. Given the larger magnitude of the drop, we conclude that the refractive index matching effect is the stronger one in highly scattering media.

Glycerol comparison measurement results shown in Fig. 8 reveal that the same solution refractive index change (from 1.337 to 1.355) due to presence of optically active glucose (243 to 900mM glucose concentration) and optically inactive glycerol (2.49 to 17.5%) yields similar amounts of optical rotation detected at θ=0 , 45, and 135° in media with μs100cm1 . This implies that the refractive index matching effect dominates over optical activity in its contribution to optical rotation in highly turbid media, as is also suggested by Fig. 7b. A similar conclusion was reached by Hadley for media with lower scattering.29

Fig. 8

Comparison measurements in highly turbid suspensions (μs100cm1) containing glycerol or glucose. For the same amount of refractive index change caused by both optically active glucose and optically inactive glycerol, similar amounts of optical rotation in both forward and backward detection geometries are seen. For the range of index matching values shown, the scattering coefficient changes from 98.3cm1 (nmedium=1.337) to 92.7cm1 (nmedium=1.355) . The symbols are data points; the lines are a guide for the eye.




We have demonstrated that the polarization parameters α, βlinear and overall scattered intensity I , simultaneously obtained from polarimetric measurements in different detection directions around a scattering sample, depend in a varied and complex way on a number of factors, including sample glucose concentration. For example, in scattering media with μs in the 30 to 100cm1 range, large changes in α and βlinear were detected in the backward direction as glucose levels increased to 1M . The three parameters ( α , βlinear , and I ) form a potentially attractive three-channel probe for glucose level monitoring in turbid media, although their dependence on sample properties and experimental geometry is complicated. The weak dependence of the derived optical rotation on the chiral nature of glucose in multiply scattering samples, while surprising and disappointing initially, may actually remove some of the difficulties of polarimetry caused by other chiral confounders in biological tissues. The range of validity of the observed trends is currently being explored in other chiral scattering systems. Further research will focus on the optimal detection geometry for maximizing sensitivity to glucose concentration levels closer to the physiological level. Other alternative detection geometries, such as out-of-plane detection, are beginning to be explored.30 Expanding upon the presented methodology, the effects of media absorption, linear birefringence, and the possibility of multiwavelength (spectroscopic) polarimetry coupled with chemometric analysis are also under investigation.



We have demonstrated the feasibility of studying the polarization states of scattered light from turbid media in different directions using a refined linear Stokes polarimeter. An extensive experimental study was performed, concentrating on sample (turbidity, chirality) and system (detection angle) variables. Selected validation with polarization-sensitive MC simulation model was performed. It is found that in turbid media, the scattering process itself can cause a large optical rotation, even in the absence of glucose. In general, the magnitude of the optical rotation is strongly dependent on sample turbidity, detection angle, and glucose levels. For example, in highly turbid media (μs=100cm1) , optical rotation can reach up to 20° levels and can exhibit large changes ( Δα up to 10 to 15° ) over the 1M glucose concentration range. The surviving linear polarization fraction decreases with increasing turbidity and also with the detection angle for low and moderately turbid media. Its angular dependence becomes less pronounced with increasing turbidity and, in fact, reverses for highly turbid samples. It is demonstrated that the three derived polarization parameters ( α,βlinear , and I ) depend on glucose levels, albeit in a complicated way. The sensitivity to glucose is strongly detection-angle–dependent, suggesting that the backward detection angles may be preferable for glucose sensing in highly turbid media, such as biological tissue, although the complicated angular and turbidity dependence of optical rotation α may make this metric less glucose-specific. Comparison measurements in turbid media containing (optically inactive) glycerol indicate that the optical rotation dependence on (optically active) glucose concentration is not due to its chiral nature, but is mainly caused by the refractive index matching effect.


The authors thank Dr. Theodore Oakberg (Hinds Instruments, Hillsboro, Oregon) for helpful discussions and the Natural Sciences and Engineering Research Council of Canada (NSERC) for financial support of this research.



J. A. Tamada, M. Lesho, and M. J. Tierney, “Keeping watch on glucose: new monitors help fight the long-term complication of diabetes,” IEEE Spectrum, 39 52 –57 (2002). 0018-9235 Google Scholar


A. J. Berger, “Minimally invasive technologies for optical glucose monitoring,” Opt. Photonics News, 12 S8 –S11 (2001). 1047-6938 Google Scholar


D. Côté and I. A. Vitkin, “Balanced detection for low-noise precision polarimetric measurements of optically-active, multiply-scattering tissue phantoms,” J. Biomed. Opt., 9 213 –220 (2004). 1083-3668 Google Scholar


O. S. Khalil, “Spectroscopic and clinical aspects of noninvasive glucose measurements,” Clin. Chem., 45 165 –177 (1999). 0009-9147 Google Scholar


K. V. Larin, M. S. Eledrisi, M. Motamedi, and R. O. Esenaliev, “Noninvasive blood glucose monitoring with optical coherence tomography,” Diabetes Care, 25 2263 –2267 (2002). 0149-5992 Google Scholar


R. R. Ansari, S. Bockle, and L. Rovati, “New optical scheme for a polarimetric-based glucose sensor,” J. Biomed. Opt., 9 103 –115 (2004). 1083-3668 Google Scholar


M. Kohl and M. Cope, “Influence of glucose concentration on light scattering in tissue-simulating phantoms,” Opt. Lett., 19 2170 –2172 (1994). 0146-9592 Google Scholar


J. S. Maier, S. A. Walker, S. Fantini, M. A. Franceschini, and E. Gratton, “Possible correlation between blood glucose concentration and the reduced scattering coefficient of tissues in the near infrared,” Opt. Lett., 19 2062 –2064 (1994). 0146-9592 Google Scholar


J. T. Bruulsema, J. E. Hayward, T. J. Farrell, M. S. Patterson, L. Heinemann, M. Berger, T. Koschinsky, J. Sandahl-Christiansen, H. Orskov, M. Essenpreis, G. Schmelzeisen-Redeker, and D. Böcker, “Correlation between blood glucose concentration in diabetics and noninvasively measured tissue optical scattering coefficient,” Opt. Lett., 22 190 –192 (1997). 0146-9592 Google Scholar


R. J. McNichols and G. L. Coté, “Optical glucose sensing in biological fluids: an overview,” J. Biomed. Opt., 5 5 –16 (2000). 1083-3668 Google Scholar


M. Mehrübeoğlu, N. Kehtarnavaz, S. Rastegar, and L. V. Wang, “Effect of molecular concentrations in tissue-simulating phantoms on images obtained using diffuse reflectance polarimetry,” Opt. Express, 3 286 –297 (1998). 1094-4087 Google Scholar


R. C. N. Studinski and I. A. Vitkin, “Methodology for examining polarized light interactions with tissues and tissuelike media in the exact backscattering direction,” J. Biomed. Opt., 5 330 –337 (2000). 1083-3668 Google Scholar


M. P. Silverman, W. Strange, J. Badoz, and I. A. Vitkin, “Enhanced optical rotation and diminished depolarization in diffusive scattering from a chiral liquid,” Opt. Commun., 132 410 –416 (1996). 0030-4018 Google Scholar


I. A. Vitkin, R. D. Lazslo, and C. L. Whyman, “Effects of molecular asymmetry of optically active molecules on the polarization properties of multiply scattered light,” Opt. Express, 10 222 –229 (2002). 1094-4087 Google Scholar


I. A. Vitkin and E. Hoskinson, “Polarization studies in multiply scattering chiral media,” Opt. Eng., 39 353 –362 (2000). 0091-3286 Google Scholar


A. J. Hunt and D. R. Huffman, “A new polarization-modulated light scattering instrument,” Rev. Sci. Instrum., 44 1753 –1762 (1973). 0034-6748 Google Scholar


P. S. Hauge, “Recent developments in instrumentation in ellipsometry,” Surf. Sci., 96 108 –140 (1980). 0039-6028 Google Scholar


J. C. Kemp, G. D. Henson, C. T. Steiner, and E. R. Powell, “The optical polarization of the sun measured at a sensitivity of parts in 10 million,” Nature (London), 326 270 –273 (1987). 0028-0836 Google Scholar


R. Anderson, “Measurement of Mueller Matrices,” Appl. Opt., 31 11 –13 (1992). 0003-6935 Google Scholar


M. Diem, G. M. Roberts, O. Lee, and A. Barlow, “Design and performance of an optimized dispersive infrared dichrograph,” Appl. Spectrosc., 42 20 –27 (1988). 0003-7028 Google Scholar


C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles, John Wiley, New York (1983). Google Scholar


M. Tarumi, M. Shimada, T. Murakami, M. Tamura, M. Shimada, H. Arimoto, and Y. Yamada, “Simulation study of in vitro glucose measurement by NIR spectroscopy and a method of error reduction,” Phys. Med. Biol., 48 2373 –2390 (2003). 0031-9155 Google Scholar


H. Liu, B. Beauvoit, M. Kimura, and B. Chance, “Dependence of tissue optical properties on solute-induced changes in refractive index and osmolarity,” J. Biomed. Opt., 1 200 –211 (1996). 1083-3668 Google Scholar


J. Y. Qu and B. C. Wilson, “Monte Carlo modeling studies of the effect of physiological factors and other analytes on the determination of glucose concentration in vivo by near infrared optical absorption and scattering measurements,” J. Biomed. Opt., 2 319 –325 (1997). 1083-3668 Google Scholar


D. Côté and I. A. Vitkin, “Robust concentration determination of optically active molecules in turbid media with validated three-dimensional polarization sensitive Monte Carlo calculations,” Opt. Express, 13 148 –163 (2005). 1094-4087 Google Scholar


B. Kaplan, G. Ledanois, and B. Drévillon, “Muller Matrix of dense polystyrene latex sphere suspensions: measurements and Monte Carlo simulation,” Appl. Opt., 40 2769 –2777 (2001). 0003-6935 Google Scholar


F. Jaillon and H. Saint-Jalmes, “Description and time reduction of a Monte Carlo code to simulate propagation of polarized light through scattering media,” Appl. Opt., 42 3290 –3296 (2003). 0003-6935 Google Scholar


L. Wang, S. L. Jacques, and L. Zheng, “MCML-Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed., 47 131 –146 (1995). 0169-2607 Google Scholar


K. C. Hadley and I. A. Vitkin, “Optical rotation and linear and circular depolarization rates in diffusively scattered light from chiral, racemic and achiral turbid media,” J. Biomed. Opt., 7 291 –299 (2002). 1083-3668 Google Scholar


X. Guo, M. F. G. Wood, I. A. Vitkin, “Effects of detection geometry on polarimetric measurements of scattered light from turbid media containing optically active glucose molecules,” Proc. SPIE, 5969 59691K (2005). 0277-786X Google Scholar
© (2006) Society of Photo-Optical Instrumentation Engineers (SPIE)
Xinxin Guo, Michael Wood, and Alex I. Vitkin "Angular measurements of light scattered by turbid chiral media using linear Stokes polarimeter," Journal of Biomedical Optics 11(4), 041105 (1 July 2006).
Published: 1 July 2006

Back to Top