* Address all correspondence to: Nirmalya Ghosh, IISER Kolkata, Mohanpur Campus, Department of Physical Sciences, PO BCKV, Kolkata, West Bengal 741252, India; Tel: 91-9734678247; Fax: 91 33 2587 3020; E-mail: nghosh@iiserkol.ac.in.

# Tissue polarimetry: concepts, challenges, applications, and outlook

**Nirmalya Ghosh**

Indian Institute of Science Education and Research (IISER), Department of Physical Sciences, Kolkata, Mohanpur, West Bengal, India

**I. Alex Vitkin**

Division of Biophysics and Bioimaging, Ontario Cancer Institute, Princess Margaret Hospital, 7th Floor Rm. 415, 610 University Ave, Toronto, Ontario, Canada M5G 2M9

Department of Medical Biophysics, University of Toronto, 610 University Avenue, Toronto, Ontario M5G 2M9 Canada

Department of Radiation Oncology, University of Toronto, 610 University Avenue, Toronto, Ontario M5G 2M9, Canada

*J. Biomed. Opt*. 16(11), 110801 (November 14, 2011). doi:10.1117/1.3652896

#### Open Access

## Abstract

Polarimetry has a long and successful history in various forms of clear media. Driven by their biomedical potential, the use of the polarimetric approaches for biological tissue assessment has also recently received considerable attention. Specifically, polarization can be used as an effective tool to discriminate against multiply scattered light (acting as a gating mechanism) in order to enhance contrast and to improve tissue imaging resolution. Moreover, the intrinsic tissue polarimetry characteristics contain a wealth of morphological and functional information of potential biomedical importance. However, in a complex random medium-like tissue, numerous complexities due to multiple scattering and simultaneous occurrences of many scattering and polarization events present formidable challenges both in terms of accurate measurements and in terms of analysis of the tissue polarimetry signal. In order to realize the potential of the polarimetric approaches for tissue imaging and characterization/diagnosis, a number of researchers are thus pursuing innovative solutions to these challenges. In this review paper, we summarize these and other issues pertinent to the polarized light methodologies in tissues. Specifically, we discuss polarized light basics, Stokes–Muller formalism, methods of polarization measurements, polarized light modeling in turbid media, applications to tissue imaging, inverse analysis for polarimetric results quantification, applications to quantitative tissue assessment, etc.

## Introduction

Polarized light has played important roles in our understanding of the nature of electromagnetic waves,^{1} elucidating the three-dimensional characteristics of chemical bonds,^{2} uncovering the asymmetric (chiral) nature of biological molecules,^{3} determining sugar concentrations in industrial processes,^{4} quantifying protein properties in solutions,^{5} supplying a variety of nondestructive evaluation methods,^{6} developing advanced concepts such as polarization entropy,^{7} contributing to remote sensing in meteorology and astronomy,^{8- 9} and differentiating between normal and precancerous cells in superficial tissue layers,^{10} as well as other biomedical applications.^{11- 12} Traditional polarimetry is well suited for applications in clear media and for studies of surfaces; however, multiple scattering in optically thick turbid media such as biological tissues causes extensive depolarization that confounds the established techniques. Further, even if some residual polarization signal can be measured,^{13} multiple scattering also alters the polarization state; for example, by scattering-induced diattenuation and by scattering-induced changes in the orientation of the linear polarization vector which appears as optical rotation.^{14} The presence of other simultaneously-occurring polarization effects further compromise quantitative polarimetry; in tissues, these include linear birefringence due to anisotropic muscle fibers and structural proteins, and optical rotation due to optically active (chiral) molecules and structures.^{15- 16} Thus, although a wealth of interesting tissue properties can potentially be probed with polarized light, accurate measurements and data analysis leading to unique interpretation of the polarization parameters are difficult, hindering the utility of polarimetric bulk tissue characterization studies.

Some inroads in biomedical polarimetry have been made in the context of optical imaging, specifically using polarization gating to separate out and potentially remove the multiply scattered (depolarized) component of the light beam in order to enhance contrast and to improve tissue imaging resolution.^{12} As discussed below, this has proven moderately successful in selected applications, provided that proper attention is paid to the optimal choice of incident polarization states (e.g., linear versus circular), polarization detection schemes (e.g., Stokes versus Mueller polarimetry), geometry of detection (e.g., transmission versus reflection), etc. For example, some promising results have been reported in skin imaging, whereby a dermatologist uses polarization imaging to selectively concentrate on either surface irregularities or alternatively on deeper epidermal/dermal layers.^{17} Further, polarization effects are extensively used in various forms of biomedical light microscopy, where one is dealing with thin fixed tissue slices, thus obviating the complicating effects of multiple scattering in bulk tissues by physical sectioning. But overall, the full potential of polarization imaging in biomedicine has not been realized, for reasons similar to the polarimetric tissue characterization studies mentioned above. To summarize, these include: 1. extensive loss of polarization signal engendered by tissue multiple scattering, 2. complicated nature of polarization effects in tissue, including simultaneous multiple effects, 3. difficulties in measuring typically small tissue polarization signals, 4. challenges in analysis and quantification of measured signals or images, 5. complexities in understanding and interpreting tissue polarimetry results, and 6. scarcity of data on detailed polarimetric properties of various tissues and their effects on polarized light propagation.

Driven by polarimetry's biomedical potential, a number of researchers are pursuing innovative solutions to these challenges. In this review paper, we summarize these and other issues pertinent to the polarized light methodologies in tissues. Specifically, we discuss polarized light basics, Stokes–Muller formalism, methods of polarization measurements, polarized light modeling in turbid media, applications to tissue imaging, inverse matrix analysis for polarimetric results quantification and applications to quantitative tissue assessment, etc. The intent of the paper is to explain the basics of polarimetry, summarize its current state of research, provide selected illustrative examples, facilitate insight and understanding of the observed trends and findings, indicate the inconsistencies and outstanding issues, and point toward the future of the field. Biomedical polarimetry is still at a relatively early stage of development, with much of its promising potential currently unrealized; it is hoped that this paper will stimulate new ideas and encourage further research into this promising field.

## Polarization Basics and the Different Polarimetry Formalisms

Definitions of polarized light and its properties are described in voluminous literature.^{6- 7,18- 19} Briefly, polarization is a property that arises out of the transverse (and vector) nature of the electromagnetic radiation and it describes the shape and the orientation of the locus of the electric field vector (

^{1,6}The superposition of many such photon states yield the resulting polarization observed for the classical wave. Thus the quantum mechanical view of polarization and the corresponding classical formalisms are mutually consistent.

Mathematical formalisms (in the classical approach) dealing with propagation of polarized light and its interaction with any optical system can be described by two formalisms; the Jones calculus^{1,6,18} which is a field-based representation (assumes coherent addition of the amplitudes and phases of the waves) and the Stokes–Mueller calculus^{1,6,18- 21} that is an intensity-based representation (assumes an incoherent addition of wave intensities). A major drawback of the Jones formalism is that it deals with pure polarization states only and cannot handle partial polarizations and thus depolarizing interactions (which are common in biological tissues). Thus its use in tissue polarimetry is limited. The general case, which does include polarization loss, can be better addressed by the Stokes–Mueller formalism as described below.

In this formalism, the polarization state of the light beam is represented by four measurable quantities (intensities) grouped in a 4 × 1 vector, known as the Stokes vector^{6,18- 21} (introduced by Stokes in 1852). The four Stokes parameters are defined relative to the following six intensity measurements (*I*) performed with ideal polarizers: *I*_{H}, horizontal linear polarizer (0 deg); *I*_{V}, vertical linear polarizer (90 deg); *I*_{P}, 45 deg linear polarizer; *I*_{M}, 135 deg (−45 deg) linear polarizer; *I*_{R}, right circular polarizer, and *I*_{L}, left circular polarizer. The Stokes vector (**S**) is defined as^{18- 21}

*I, Q, U*, and

*V*are Stokes vector elements.

*I*is the total detected light intensity which corresponds to addition of the two orthogonal component intensities,

*Q*is the difference in intensity between horizontal and vertical polarization states,

*U*is the portion of the intensity that corresponds to the difference between intensities of linear +45 deg and −45 deg polarization states, and

*V*is the difference between intensities of right circular and left circular polarization states.

In the Stokes formalism, the following polarization parameters of any light beam are defined:^{6,18- 21}

Net degree of polarization

While the Stokes vectors represent the polarization state of light, a 4 × 4 matrix **M,** known as the Mueller matrix (after its inventor Hans Mueller in the 1940s) describes the transfer function of any medium in its interaction with polarized light:^{6,18- 21}

**S**

_{i}and

**S**

_{o}being the Stokes vectors of the input and output light, respectively. The 4 × 4 real Mueller matrix

**M**possesses at most 16 independent parameters (or 15 if the absolute intensity is excluded), including depolarization information. All the medium polarization properties are encoded in the various elements of the Mueller matrix, which can thus be thought of as the complete “optical polarization fingerprint” of a sample.

The fundamental requirement that real Mueller matrices must meet is that they map physical incident Stokes vectors into physically resultant Stokes vectors [satisfying Eq. 5].^{20- 21} The conditions for physical realizability of **M**, their associated interpretations, the relationships between the two formalisms (Jones and the Stokes–Mueller) have been discussed in the literature.^{20- 22}

Although both the Jones and Stokes–Mueller approaches rely on linear algebra and matrix formalisms, they are different in many aspects. Specifically, the Stokes–Mueller formalism has certain advantages. First of all, it can encompass any polarization state of light, whether it is natural, totally, or partially polarized (can thus deal with both polarizing and depolarizing optical systems). Secondly, the Stokes vectors and Mueller matrices can be measured with relative ease using intensity-measuring conventional (square-law detector) instruments, including most polarimeters, radiometers, and spectrometers. Since biological tissue is a turbid medium where significant depolarization is encountered due to strong multiple scattering effects, the Stokes–Mueller formalism has been used in most tissue polarimetry applications. In contrast, the use of the Jones formalism has been limited as a complementary theoretical approach to the Mueller matrix calculus, or to studies in clear media, specular reflections, and thin films where polarization loss is not an issue. In this paper, we review the use of the Stokes–Mueller approach for noninvasive assessment of biological tissues, discuss inverse analysis methods for extraction/quantification of the intrinsic tissue polarimetry characteristics, and provide selected illustrative application examples of tissue polarimetry. In the following, we define the basic medium polarization properties through the Mueller matrix formalism.

If an incident state is 100% polarized and the exiting state has a degree of polarization less than unity, then the system is said to possess depolarization property. Depolarization is usually encountered due to multiple scattering of photons (although randomly oriented uniaxial birefringent domains can also depolarize light); incoherent addition of amplitudes and phases of the scattered field results in scrambling of the output polarization state. The general form of a pure depolarization Mueller matrix is^{20,23}

*a*| and 1 − |

*b*| are depolarization factors for linear polarization (horizontal/vertical and +45 deg/−45 deg linear polarizations, respectively) and 1 − |

*c*| is the depolarization factor for circular polarization. The net depolarization factor is usually defined as

_{light}∼ 1 − Δ

_{medium}(and similarly for linear and circular states), but this is strictly an equality only in special cases. In general though, a medium with a high value of Δ is significantly depolarizing, and the degree of polarization of the light after interacting with it will be quite low.

The next two polarization effects, retardance and diattenuation, arise from differences in refractive indices for different polarization states, and are often described in terms of ordinary and extraordinary indices and axes. Differences in the real parts of the refractive index lead to linear (circular) birefringence, whereas differences in the imaginary parts of the refractive index cause linear (circular) dichroism (which manifests itself as diattenuation, described below). Specifically, retardance is the phase shift between two orthogonal polarizations of the light. Linear retardance δ (birefringence) arises due to a difference in phase between orthogonal linear polarization states (between vertical and horizontal, or between 45 deg and −45 deg). Circular retardance ψ (optical rotation) arises due to difference in phase between right circularly polarized (RCP) and left circularly polarized (LCP) states. The general form of a Mueller matrix of a pure linear retarder with retardance δ and fast axis oriented at an angle θ with respect to the horizontal is^{20,23}

Similarly, the Mueller matrix for a circular retarder with optical rotation value ψ is^{20,23}

Diattenuation (*d*) of an optical element corresponds to differential attenuation of orthogonal polarizations for both linear and circular polarization states. Accordingly, linear diattenuation is defined as differential attenuation of two orthogonal linear polarization states and circular diattenuation is defined as differential attenuation of RCP and LCP. Mathematically, the Mueller matrix for an ideal diattenuator can be defined using the parameters, *q* and *r,* intensity transmittance (or reflectance) for the two incident orthogonal polarization states (either linear or circular), and the orientation angle of the principal axis (θ). Using this convention, the general Mueller matrix for a linear diattenuator is defined as^{20,23}

Similarly for circular diattenuation, the general form of the Mueller matrix is^{20,23}

*d*= 1 for ideal polarizer), although often with a significant reduction in the overall intensity

*I*. Note that diattenuation is analogous to dichroism, which is defined as the differential absorption of two orthogonal linear polarization states (linear dichroism) or RCP−LCP (circular dichroism), but is more general in a sense that it is defined in terms of differential attenuation (either by absorption or scattering).

## Experimental Tissue Polarimetry Systems

As outlined in Sec. 1, biomedical polarimetry research has two major directions, tissue imaging and tissue characterization. First, polarization can be used as an effective tool to discriminate against multiply scattered light and thus can facilitate higher resolution imaging of tissue and its underlying structure.^{17} Moreover, the intrinsic polarimetry characteristics themselves contain a wealth of morphological, biochemical, and functional information that can be exploited for noninvasive and quantitative tissue diagnosis.^{15- 16,24} For either of these applications, accurate measurement of the polarization retaining signal is extremely important. In this regard, many of the traditional polarimetry systems are not suitable for biological tissue examination (e.g., crossed linear polarizers used in microscopy for examining thin fixed *ex vivo* tissue slices). This follows because multiple scattering in thick tissues leads to depolarization of light, creating a large depolarized source of noise that hinders the detection of the small remaining information-carrying polarization signal. A variety of experimental tools have therefore been developed to maximize measurement sensitivity, so that reliable measurements and analyses of the tissue polarimetry data can be performed. These methods can be employed to perform measurement of both Stokes vector of the light upon interacting with the sample, and/or of the Mueller matrix of the sample itself, as described below.

As mentioned in Sec. 2b, the four Stokes parameters of light can be determined by performing six intensity measurements involving linear and circular polarization states (*I*_{H}, *I*_{V}, *I*_{P}, *I*_{M}, *I*_{R}, and *I*_{L}).^{21} Alternatively, this can be achieved by just four measurements, exploiting the property (*I*_{H} + *I*_{V} = *I*_{P} + *I*_{M} = *I*_{L} + *I*_{R}).^{25} In this approach, a circular polarizer is designed consisting of a linear polarizer whose transmission axis is set at +45 deg with respect to the horizontal direction, followed by a quarter wave-plate with its fast axis parallel to the horizontal direction. Three sets of intensity measurements [*I*_{Cir} (α)] are performed by changing the angle (α) of this circular polarizer to 0, 45, and 90 deg with respect to the horizontal axis. The combined polarizer is then flipped to the other side and the final intensity measurement [*I*_{Lin} (α)] is made by setting α at 0 deg. The four Stokes parameters can be determined from the measured intensities as

^{25}Although this method has been employed in some experimental depolarization studies to measure Stokes parameters of light transmitted (or backscattered) from tissue and tissue-like turbid media,

^{26- 30}a more sensitive detection scheme is desirable, specifically in applications involving accurate quantification of the intrinsic tissue polarimetry characteristics.

^{15- 16}One possible method for improving the sensitivity of the measurement procedure is the use of polarization modulation with synchronous detection. Many sensitive detection schemes are possible with this approach.

^{31- 36}Some of these perform polarization modulation on the light that is incident on the sample; others modulate the sample-emerging light, by placing the polarization modulator between the sample and the detector. The resultant signal can be analyzed to yield sample-specific polarization properties that can then be linked to the quantities of interest.

By way of illustration, a schematic of the experimental polarimetry system employing polarization modulation and synchronous lock-in-amplifier detection is shown in Fig. 1.^{31} Unpolarized light from a laser is used to seed the system. The light first passes through a mechanical chopper operating at a frequency *f*_{c} ∼ 500 Hz; this is used in conjunction with lock-in amplifier detection to accurately establish the overall signal intensity levels. The input optics [a linear polarizer *P*_{1} with/without the quarter wave-plate (QWP_{1})] enables generation of any of the four input polarization states, 0 deg (Stokes vector [1 1 0 0] ^{T}), 45 deg (Stokes vector [1 0 1 0] ^{T}), and 90 deg (Stokes vector [1 −1 0 0]^{T}) linear polarizations, as well as circular polarization (Stokes vector [1 0 0 1] ^{T}) incident on the sample. Following light-tissue interactions, the detection optics begin with a removable quarter wave-plate (QWP_{2}) with its fast axis oriented at −45 deg, when in place allowing for the measurement of Stokes parameters *Q* and *U* (linear polarization descriptors), and when removed allowing for the measurement of Stokes parameter *V* (circular polarization descriptor). The tissue-scattered light then passes through a photoelastic modulator (PEM), which is a linearly birefringent resonant device operating in the kilohertz range (e.g., *f*_{p} = 50 kHz). The fast axis of the PEM is at 0 deg and its retardation is modulated according to the sinusoidal function *δ*_{PEM} (*t*) = *δ*_{o} sin*ωt*, where *ω*_{p} = 2π*f*_{p} and *δ*_{o} is the user-specified amplitude of maximum retardation of PEM. The light finally passes through a linear analyzer orientated at 45 deg, converting the PEM-imparted polarization modulation to an intensity modulation suitable for photodetection. The resulting modulated intensity is collected using a pair of lenses and is relayed to an avalanche photodiode detector. The detected signal is sent to a lock-in amplifier with its reference input toggling between the frequencies of the chopper (500 Hz) and the PEM controller (50 kHz and harmonics) for synchronous detection of their respective signals.

**F1 :**

(a) A schematic of the experimental polarimetry system employing polarization modulation and synchronous lock-in-amplifier detection. C, mechanical chopper; P_{1}, P_{2}, polarizers; QWP_{1}, QWP_{2}, removable quarter wave-plates; L_{1}, L_{2} lenses; PEM, photoelastic modulator; APD, photodetector; *f*_{c}, *f*_{p} modulation frequencies of mechanical chopper and PEM, respectively. The detection optics can be rotated by an angle γ around the sample. (Adopted from Ref. ^{16},) (b) A schematic of the liquid crystal variable retarder polarimeter. P_{1}, P_{2}, linear polarizers; LC_{1}, LC_{2}, LC_{3}, LC_{4}, liquid crystal variable retarders. P_{1}, LC_{1} (with retardance of δ_{1} having orientation angle θ_{1}), LC_{2} (with retardance of δ_{2} having orientation angle θ_{2}) form the (PSG) unit. LC_{3} (with retardance of δ_{2} having orientation angle θ_{2}) LC_{4} (with retardance of δ_{1} having orientation angle θ_{1}), and P_{2} form the (PSA) unit. The schematic is shown for transmission measurements; other detection geometries are also possible using this scheme.

For this particular experimental arrangement, the Stokes vector of light after the analyzing block [*I*_{f}* Q*_{f}* U*_{f}* V*_{f}]^{T} can be related to that of the sample-emerging beam [*I Q U V*]^{T} (with detection quarter wave-plate in place) using Mueller matrix algebra^{31}

*q = Q/I*,

*u = U/I*, and

*v = V/I*, and

*δ*is the time-dependent PEM retardation, δ = δ

_{0}sinω

*t*. By expanding the time varying parts of Eqs. 17- 18 in Fourier series of Bessel functions, and by setting the peak retardance of the PEM

*δ*

_{o}= 2.405 radians (the retardance value at which the zeroth-order Bessel function exhibits its first null), one can relate the normalized Stokes parameters to the synchronously-detected signals at the chopper frequency

*V*

_{1fc}(the dc signal level), and at the first and second harmonics of the PEM frequency

*V*

_{1fp}and

*V*

_{2fp}. When the detection quarter wave-plate in place, this yields

^{31}

and when the detection quarter wave-plate is removed

^{31}

With some additional measurements and analysis, one can go beyond polarimetric light description (Stokes vectors) and determine the polarization transfer function of the sample itself (Mueller matrix). For 4 × 4 Mueller matrix determination, both dc (involving sequential static measurements) and ac modulation-based measurement procedures have been employed. In fact, the polarization modulation approach described in Sec. 3a can also perform sensitive Mueller polarimetry. This can be achieved by sequentially cycling the input polarization between four states (linear polarization at 0 deg, 45 deg, 90 deg, and right circular polarization) and by measuring the output Stokes vector for each respective input states. The elements of the resulting four measured Stokes vectors (16 values) can be algebraically manipulated to solve for the sample Mueller matrix:^{37}

*H*(0 deg),

*P*(45 deg),

*V*(90 deg), and

*R*(right circularly polarized; left circular incidence can be used as well, resulting only in a sign change). The indices

*i,j*= 1, 2, 3, 4 denote rows and columns, respectively.

The described experimental embodiment of the polarization modulation/synchronous detection approach has been used by us for both Stokes vector and Mueller matrix measurements in complex tissue-like turbid media and in actual tissues;^{37- 38} some of the results are presented subsequently in this paper.

Among the various other modulation-based Mueller matrix polarimeters, the dual rotating retarder approach has been widely used in tissue polarimetry investigations.^{39- 41} In this scheme, polarization modulation of the incident state is generated by passing light first through a fixed linear polarizer and then through a rotating linear retarder (retardation δ_{1}) with angular speed ω_{1}. The analyzing optics contains another rotating retarder (retardation δ_{2}, synchronously rotating at angular speed ω_{2}) and a fixed linear polarizer. In the usual configuration, the retardation values of the two retarders are chosen to be the same δ_{1} = δ_{2} = π/2, the axis of the polarizer and the analyzer are kept parallel, and the angular rotation speeds of the retarders are kept as ω_{1} = ω and ω_{2} = 5ω, respectively.^{39,41} The rotation of the retarders at these different rates results in a modulation of the detected intensity signal, as can be understood by sequentially writing the Mueller matrices corresponding to each optical element (polarizers, retarders, and the sample). Note that for this specific scheme, the modulation in the detected intensity arises due to harmonic variation of the orientation angle of the two retarders kept at the polarizing and analyzing end of the polarimeter (θ and 5θ, respectively). It has been shown that the five to one ratio of angular rotation speeds of the two retarders encodes all 16 Mueller matrix elements onto the amplitude and phases of 12 frequencies in the detected intensity signal. The detected signal is Fourier analyzed and the Mueller matrix elements are constructed from the Fourier coefficients. A more general approach based on this scheme may employ arbitrary values for linear retardations (δ_{1} and δ_{2}), rotation speed ratios (ω_{1}/ω_{2}), and axis of the linear polarizers, depending upon which elements (rows/columns) of Mueller matrix are given a priority in terms of higher determination precision and/or SNR.^{39,41}

Snapshot Mueller matrix polarimeter is another important development in modulation-based Mueller matrix polarimetry.^{42} The approach exploits wavelength polarization coding and decoding for high sensitivity, instantaneous measurement of all 16 Mueller matrix elements simultaneously. Briefly, several wavelength-encoded polarization states are generated with a broadband spectrum source, two birefringent retarders, and a linear polarizer. The wavelength decoding is performed using a similar combination of birefringent retarders and a linear polarizer. The thickness of the retarders in the encoding (polarization state generator) and decoding (polarization state analyzer) systems are optimized to generate and analyze sufficient number of polarization states; the resulting spectral signals (over a narrow spectral range Δλ ∼ 10 nm) recorded using a spectrometer are Fourier analyzed to yield all 16 Mueller matrix elements.

Note that the polarization modulation-based polarimeters described above are convenient for point polarimetry measurements; large area imaging is generally not feasible using these approaches (since these typically employ a synchronous detection scheme). Such imaging polarimetry can be accomplished by dc measurements involving sequential measurements with different combinations of source polarizers and detection analyzers, albeit with lower SNR than the synchronous ac detection schemes described above. Because a general 4 × 4 Mueller matrix has 16 independent elements, at least 16 independent measurements are required;^{43} due to the low sensitivity of the dc approach, alternatives have been investigated. For example, polarimetric imaging systems based on liquid crystal variable retarders enable measurement of the Mueller matrix elements with higher sensitivity and precision.^{44- 46} A schematic of such a measurement strategy is shown in Fig. 1. The polarimetry system is comprised of a polarization state generator (PSG) and a polarization state analyzer (PSA) unit coupled to an imaging camera for spatially resolved signal detection. The PSG consists of a linear polarizer (*P*_{1}) and two liquid crystal variable retarders (LC_{1} and LC_{2} having variable retardance of δ_{1} and δ_{2}, respectively). Generally, the birefringence axes of LC_{1} and LC_{2} are kept at angles θ_{1} and θ_{2}, respectively with respect to the axis of the polarizer *P*_{1}.^{45} In a more specific arrangement (which has been widely used), the angles θ_{1} and θ_{2} are chosen to be 45 and 0 deg, respectively. The Stokes vector generated from this arrangement and incident on the sample is^{46}

_{1}and δ

_{2}. The PSA unit also consists of similar arrangements of liquid crystal variable retarders (LC

_{3}and LC

_{4}) and linear polarizer, but positioned in reverse order, and followed by a detector (for imaging applications, this is a CCD camera). The polarimetry signal analysis proceeds as follows.

The PSG output can be represented by **W**, a 4 × 4 matrix whose column vectors are the four generated Stokes vectors **S**_{i} incident on the sample. Similarly, after sample interactions, the PSA results can also be described by a 4 × 4 matrix **A**. The Stokes vectors of light to be analyzed are projected on four basis states that are the row vectors of the 4 × 4 analysis matrix **A**. For the construction of the Mueller matrix, a sequence of 16 measurements are performed. This 4 × 4 intensity measurement matrix *M*_{i} can be written as^{20,47}

**M**is the sample Mueller matrix, which when presented as 1 × 16 vector (

**M**

_{vec}), can be related to the intensity measurement (1 × 16) vector

**M**

_{ivec}as

**Q**is a 16 × 16 matrix given as Kronecker product of

**A**with transpose of

**W**

**A**and

**W**).

^{47}Based on this approach, several measurement schemes are possible. In fact, the choice of the values for retardance δ

_{1}and δ

_{2}, the orientation angles of the retarders with respect to the polarizers (analyzers) can be optimized to minimize the noise in the resulting Mueller matrix

**M**. Such optimized measurement strategies have also been explored for performing Mueller matrix measurements in tissues.

^{45}

Recently, a novel stroboscopic illumination technique has been explored to facilitate large area imaging using a polarization modulation scheme.^{48} This approach utilizes a pulsed laser diode to illuminate the object. The short current pulses of this laser diode are precisely controlled by a programmable pulse generator. The temporal reference is triggered by the controller of a PEM operating at 50 kHz. This synchronization procedure facilitates freezing of the intensity variation of the PEM modulated signal at desirable temporal phases. Measurement of the intensity (using a CCD camera) at four specific temporal phases (frozen via stroboscopic illumination) are used to deduce the two-dimensional images of the two well-known ellipsometric parameters.^{48- 49} Note, however, this approach has not yet been explored for complete (16 element) Mueller matrix imaging.

Having described the various experimental strategies for sensitive measurement of Stokes vector and Mueller matrix in turbid medium-like tissue, we now turn to another challenging problem of accurately modeling the polarization signals in turbid media, in the forward sense.

## Modeling Polarized Light Transport in Complex Turbid Media

To aid tissue polarimetry in successful implementation of polarization gated optical imaging and for quantitative determination of intrinsic tissue polarimetry characteristics, accurate forward modeling is enormously useful. This helps in gaining physical insight, designing, and optimizing experiments, and analyzing/interpreting the measured data. The use of electromagnetic theory with Maxwell's equations is the most rigorous and best-suited method for polarimetry analysis, at least in clear media with well-defined optical interfaces. However, unlike optically clear media, tissue is a turbid medium possessing microscopically inhomogeneous complex dielectric structures (macromolecules, cell organelles, organized cell structures, blood and lymphatic networks, extra-cellular matrix, interstitial layers, etc.). Due to the ensuing complexity, the Maxwell's equations approach for polarized light propagation in such a complex turbid medium is impractical and is not presently feasible.^{11,50} Instead, light propagation through such media is often modeled using the radiative transport theory.^{11,50- 51} Although the scalar radiative transport theory and its simplified approximation, the diffusion equation, has been successfully used to model light transport in tissue (specifically light intensity distribution in tissue volume, diffuse reflectance, etc.), both are intensity-based techniques, and hence typically neglect polarization.^{51} Alternatively, the vector radiative transfer equation (VRTE), which includes polarization information by describing transport of the Stokes vectors of light (photon packet) through a random medium,^{11,51} has been explored for tissue polarimetry modeling. However, solving the VRTE in real systems is rather complex. A wide range of analytical and numerical techniques have been developed to solve VRTE, namely, the small angle approximation, the transfer matrix, the singular eigenfunction, the adding-doubling, the discrete ordinates, the successive orders, and the invariant embedding methods.^{11,51} Unfortunately, these are often too slow and insufficiently flexible to incorporate the necessary boundary conditions for arbitrary geometries and arbitrary optical properties as desirable in case of tissue. A more general and robust approach is the Monte Carlo (MC) technique, as described in greater detail below. First, we present a brief overview of some of the simpler analytical approximations developed to deal with depolarization of multiply scattered waves in turbid medium. These are useful for understanding the mechanisms of depolarization of light in turbid media, for performing rough estimates and order-of-magnitude calculations, and for optimizing the polarization gating schemes for tissue imaging.

###### Modeling Depolarization of Multiply Scattered Light in Turbid Medium: Approximate Analytical/Heuristic Approaches

Various methods using photon diffusion formalisms, random walk models, maximum entropy principles, etc., have been explored for modeling of depolarization of multiply scattered waves in random medium.^{52- 59} The aim of these models have been to derive analytical relationships between various quantities of practical interest such as the degree of polarization (either linear or circular) of forward- or backscattered light from a turbid medium, average pathlengths, the optical transport parameters of the medium, and so forth. As in the case for radiative transport theory, in these models also, the turbid medium is considered to have bulk-average scattering and absorption properties, representative of isotropic tissue volumes. For this description, the turbid medium is usually modeled through the optical transport parameters, namely, the absorption coefficient (*μ*_{a}), single scattering coefficient *(μ*_{s}), and single scattering anisotropy (*g*).^{50- 51} As is known from the transport theory, the linear isotropic optical coefficients are defined so that *l*_{a} = μ_{a}^{−1} and *l*_{s} = μ_{s}^{−1} give the absorption and scattering mean free paths, respectively. The anisotropy parameter *g* is defined as the average cosine of scattering angle. The value of *g* ranges from −1 to +1, where *g* = −1 corresponds to fully backward scattering, *g* = 0 corresponds to forward-backward symmetric scattering, and *g* = +1 corresponds to fully forward scattering. In general, the value of *g* depends on the average size of the scatterers in the medium relative to the wavelength of the irradiation. For a medium comprised of scatterers whose size is much smaller than the wavelength (radius *a* ≪ λ), anisotropy parameter *g* is ∼ 0, its value approaching unity (*g* ∼ 1) for medium comprised of larger sized scatterers (*a* ≥ λ). The latter regime applies to most biological tissues in the visible/near-infrared (NIR) spectral range.^{50} The other optical transport parameter frequently used in tissue optics is the reduced scattering coefficient *μ*_{s}^{′} = *μ*_{s} (1−*g*).^{50- 51} This is relevant in the multiple scattering regime and its use assumes that the light intensity metrics (reflectance, transmittance, fluence) of a turbid volume with optical parameters *μ*_{a}, *μ*_{s}, and *g* (≠ 0) are the same as those for an analogous volume with optical parameters *μ*_{a}, *μ*_{s}^{′} and *g* = 0.^{50} The corresponding mean free path, known as the transport mean free path is defined as *l*^{*} = (*μ*_{s}^{′})^{−1}, and is referred to as the typical length scale over which the propagation direction of photons get randomized in a multiply scattering medium. Note that for forward-scattering media such as most biological tissues, *μ*_{s}^{′} < *μ*_{s} by a factor of (1−*g*), thus the transport mean free path *l*^{*} is longer than the mean free path *l* (by a factor of 1/(1−*g*), for example 10× for *g* = 0.9).

Pioneering work in turbid polarimetry modeling was carried out by Bicout et al. ^{53} They related the depolarization of light by multiple scattering to a process of entropy production. Based on the so-called maximum entropy principle, the single path (photon undergoing successive scattering events) degree of polarization decays exponentially with an increasing number of scattering events (*n*).^{52- 53} For a medium comprised of a collection of nonabsorbing, optically inactive, spatially uncorrelated spherical particles whose size is much smaller compared to wavelength (radius *a* ≪ λ, *g* ∼ 0, the so-called Rayleigh regime), expression for single path degree of linear and circular polarization (for incident linearly and circularly polarized light, respectively) can be expressed as^{53}

*l*is the scattering mean free path (note that

*g*∼ 0 here, there is no distinction between

*l*and

*l*

^{*}), and the parameters ξ

_{l}and ξ

_{c}are known as characteristic length scale of depolarization of incident linearly and circularly polarized light, respectively. For a medium comprised of ensemble of Rayleigh scatterers the values for ξ

_{l}and ξ

_{c}have the following approximate analytical form

^{53}

_{l}and ξ

_{c}) are known to depend strongly on the size of the scatterers present in the medium (thus on the value of

*g*).

^{53}In fact, empirical relationships between ξ

_{l}, ξ

_{c}, and

*g*has also been obtained based on random walk models,

^{58}radiative transfer theory,

^{55}and extended photon diffusion approximation combined with experimental measurements using diffusing wave spectroscopy (measurements of intensity fluctuations of light scattered from turbid media).

^{57}This dependence can be summarized as follows:

*n*)] as

*n*) obtained from the solution of photon diffusion equation for a given detection geometry.

^{50- 51}The resulting degree of polarization of multiply scattered light transmitted or reflected from a slab of turbid medium can be expressed as follows.

Forward scattering geometry, transmission through a slab of thickness *L* (Ref. ^{53})

*L*→ ∞ (Refs.

^{11}and

^{59})

*l** is the transport mean free path andγ is the correlation decay parameter with its value ranging between from 1.5 to 3.

^{11}

The effect of medium absorption can also be modeled by modifying ρ(*n*) to account for absorption-induced photon loss {multiplying *p(n)* by a factor exp [−μ_{a}*n*]}. Accordingly, the expression for residual degree of polarization of light backscattered from an absorbing turbid medium takes the form^{11}

*a*≪ λ,

*g*∼ 0), depolarization of circularly polarized light is stronger than linearly polarized light (ξ

*l*> ξ

_{c}). The reverse is the case (ξ

_{l}< ξ

_{c}) for media comprised of larger scatterers (

*a*≥ λ,

*g*≥ 0.7, the so-called Mie regime). This is illustrated in Fig. 2, where the computed variations of the length scales of depolarization (ξ

_{l}and ξ

_{c}) are shown as a function of size parameter of scatterer (

*X*) present in the medium. Here,

*X*= 2π

*an*

_{m}/λ is the size parameter for scatterer of radius a (varying between 0.01 to 1.11 μm) having a refractive index

*n*

_{s}(chosen to be 1.59 for the calculation) embedded in a surrounding medium with refractive index

*n*

_{m}(= 1.33).

^{19}The wavelength of light was chosen to be symbol λ = 0.6328 μm. The values for ξ

_{l}and ξ

_{c}are observed to increase with increasing value of

*X*, indicating weaker depolarization with increasing size of the scatterer. Further, as can be seen, at lower value of size parameter (

*X*< 2), ξ

_{l}is larger than ξ

_{c}(depolarization of circular polarization is stronger as compared to depolarization of linear polarization) and reverse is the case for larger size parameter values (

*X*> 2). These theoretical trends (of Fig. 2) have been confirmed by experimental depolarization studies conducted on turbid media comprised of spherical scatterers having varying sizes.

^{26- 30}The reason for the observed size parameter dependence of depolarization of light in turbid medium and the differences in relative rate of depolarization of linear and circular polarization states are discussed subsequently in Sec. 5a in context to the results of the corresponding experimental depolarization studies. Although the predicted depolarization trends have been confirmed by experimental studies for certain detection geometries (e.g., forward scattering geometry), one must be careful in generalizing the applicability of these heuristic models, and the resulting predictions/trends, to arbitrary detection geometry/direction. Specifically, the assumption of detection geometry/direction independent depolarization metrics (ξ

_{l}and ξ

_{c}, see Fig. 2) may not hold for more complex detection geometries than those discussed here. This aspect is also discussed in more detail in Sec. 5a.

**F2 :**

The theoretically computed variation of length scales of depolarization for incident linearly (ξ_{l}) and circularly (ξ_{c}) polarized light as a function of size parameter of scatterer (*X*). The values for ξ_{l} and ξ_{c} increase with increasing value of *X*, indicating weaker depolarization with increasing size of the scatterer. The general decreasing trend beyond a size parameter value of *X* ≥ 10 has been attributed to Mie resonance effects (Ref. ^{53}). See text for details.

The approximate analytical approaches described in Sec. 4a are mainly aimed at understanding the overall depolarization trends, exploring the dependence of depolarization on the scattering properties of the media, and designing general polarization schemes to discriminate against multiply scattered photons for tissue imaging in “simple” geometries. However, these approaches, while useful, are approximate by their very nature and typically neglect other simultaneously occurring complex tissue polarimetry events (such as linear birefringence, optical activity, etc.). A more encompassing, accurate method is clearly needed for further advances in tissue polarimetry. This can be accomplished by the polarization-sensitive Monte Carlo (PSMC) techniques.^{50,60}

Before we describe the polarization-sensitive Monte Carlo models, we briefly discuss the various intrinsic tissue polarimetry characteristics that must be dealt with in accurate forward modeling of polarized light-tissue interactions.

Depolarization caused by multiple scattering is the most prominent polarimetry effect in biological tissues. Multiple scattering is caused by the high density of tissue scattering centers, originating from the random fluctuations of the local refractive index in the tissue microstructure (inside the cell and in the extra-cellular matrix). In fact, the tissue scattering centers vary in size (and shape) from micrometer scale and below (sub-cellular structures such as mitochondria, ribosomes, lysosomes, Golgi apparatus, etc.) to several tens of micrometers (whole cells, collagen fibers, etc.). Typical refractive index fluctuations in these scattering structures vary from n_{s} ∼ 1.4–1.5 (the average background refractive index of cytoplasm and interstitial fluid *n*_{m} ∼ 1.34).^{50} Light scattering from all of these microscopic scattering structures contributes in a complex fashion to the observed depolarization of light in tissue. Note that the underlying mechanism of depolarization due to multiple scattering is the scrambling of photon's reference frame (scattering plane) as a consequence of random sequence of scattering events in a variety of scattering directions.

Linear retardance (or birefringence) is the other important tissue polarimetry characteristic. Although not as pervasive as multiple scattering, the anisotropic organized nature of many tissues stemming from their fibrous structure manifests as anisotropic refractive indices parallel and perpendicular to the fibers. Accordingly, these tissues exhibit linear birefringence, which is defined as the difference in refractive indices of the fibers, Δ*n* (= *n*_{e}−*n*_{o}), where *n*_{e} and *n*_{o} are extraordinary and ordinary refractive indices (the electric field vector or linear polarization of light is perpendicular and parallel to the fiber orientation). This results in phase retardation, also called retardance [δ = (2π/λ)Δ*nL*, *L* is the pathlength] between two orthogonal linear polarization states while propagation through tissue. Extra-cellular matrix proteins (collagen and elastin), actin-myosin fibers, mineralized hydroxyapatite crystals are examples of such birefringent fibers. Various types of tissues, such as muscle, skin, myocardium, bone, teeth, cornea, tendon, cartilage, eye sclera, dura mater, nerve, retina, myelin, etc., possess these birefringent (uniaxial and occasionally biaxial^{12}) fibrous structures. Typical values of linear birefringence of these biological fibers in the visible wavelength range are in the range Δ*n* = 10^{−3} to 10^{−2}.^{12} Interestingly, even though uniform uniaxial birefringence may not be a direct contributor to depolarization *per se*, randomly oriented spatial domains of uniform uniaxial birefringent properties may cause polarization loss.^{11} Similarly, circular birefringence (retardance, also called optical rotation in this context) in tissue arises due to the presence of asymmetric optically active chiral molecules like glucose, proteins, and lipids.^{24} Finally, many biological molecules (such as amino acids, proteins, nucleic acids, etc.) also exhibit dichroism or diattenuation effects. The magnitude of diattenuation effects in tissue is, however, much lower compared to the other polarization phenomena described above.

The MC technique is a general and robust approach for modeling light transport in random medium.^{50,60} In this statistical approach to radiative transfer, the multiple scattering trajectories of individual photons are determined using a random number generator to predict the probability of each scattering event. The superposition of many photon paths approaches the actual photon distribution in time and space. This approach has the advantage of being applicable to arbitrary geometries and arbitrary optical properties, including ability to simulate heterogeneous media. Most Monte Carlo models were developed for intensity calculations only and neglected polarization information; the most commonly used being the code of Wang et al. ^{60} More recently, a number of implementations have incorporated polarization into the Monte Carlo approach.^{56,61- 64}

A specific example of a polarization-sensitive Monte Carlo model is shown via a flow chart in Fig. 3.^{15} In addition to modification/incorporation of the position and propagation direction of each photon, polarization information is incorporated by keeping track of the Stokes vectors of propagating photon packets. When the photon encounters a scattering event, a scattering plane and angle are statistically sampled based on the polarization state of the photon and the Mueller matrix of the scatterer. The photon's reference frame is first expressed in the scattering plane and then transformed to the laboratory (experimentally observable) frame through multiplication by appropriate rotation matrices and the Mueller matrix calculated through Mie scattering theory again. For *n* number of successive scattering events, the resulting Stokes vector (**S**_{f}) for a photon packet can be expressed as

*S*

_{O}is the input Stokes vector, and θ and ϕ are the scattering and the azimuthal angles respectively.

*R*(ϕ

_{i}) is the rotation matrix (connecting the two Stokes vectors that describe the same polarization state with respect to the reference plane and the scattering plane) for the

*i*’th scattering event (

*i*= 1,2 ….

*n*) given as

*M*(θ

_{i}) is the Mie-theory calculated scattering Mueller matrix for the

*i*’th scattering event (defined in the scattering plane). For isotropic spherical scatterers, symmetry considerations simplify the general form of

*M*(θ) to yield

^{19}

*M*

_{11}=

*M*

_{22},

*M*

_{12}=

*M*

_{21},

*M*

_{33}=

*M*

_{44}and

*M*

_{43}= −

*M*

_{34}). On the other hand, for scatterers having arbitrary shapes, the form of

*M*(θ) is far more complex, essentially having nonzero values for all the matrix elements.

^{19}

Using the aforementioned approach, the evolution of polarization state of each photon packet is tracked following successive scattering events. Absorption effects are handled as in intensity-based Monte Carlo models.^{50,60} Upon encountering an interface (either an internal one, representing tissue domains of different optical properties, or an external one, representing external tissue boundary), the probability of either reflection or transmission is calculated using Fresnel coefficients. As no coherence effects are considered, the final Stokes vector for light exiting the sample in a particular direction is computed as the sum of all the Stokes vectors of appropriate directional photon sub-populations. As previously described, algebraic manipulation of the resulting Stokes vectors for a variety of different polarization inputs can be performed to yield the Mueller matrix of the simulated turbid medium.

In the PSMC approach that handles Stokes vector evolution due to absorption and scattering in a manner just described, other polarimetry effects such as linear birefringence and optical activity can, in principle, be incorporated by including their corresponding Mueller matrices in Eq. 34. However, this is not an obvious modeling step. Difficulty arises in formulating simultaneous polarization effects. Matrix multiplication of the Mueller matrices for individual polarization effects is not commutative (**M**_{A}** M**_{B} ≠ **M**_{B}** M**_{A,} or the Hermitian is nonzero); thus, different orders in which these effects are applied will have different effects on the polarization. Ordered multiplication of these matrices in fact does not make physical sense, as in biological tissue, these effects (such as optical activity due to chiral molecules and linear birefringence due to anisotropic tissue structures) are exhibited simultaneously and not one after the other as sequential multiplication implies. Fortunately, there exists a method to simulate simultaneous polarization effect in clear media through the so-called *N*-matrix formalism, which combines the effects into a single matrix describing them simultaneously.^{6,65}

The *N*-matrix approach was first developed by Jones,^{65} and a more thorough derivation is provided in Kliger et al. ^{6} Briefly, in this approach, the matrix of the sample is represented as an exponential function of a sum of matrices, where each matrix in the sum corresponds to a single optical polarization effect. The issue of ordering of noncommutative matrices is overcome as matrix addition is always commutative, and applies to differential matrices representing the optical property over an infinitely small optical pathlength. These differential matrices are known as *N*-matrices, and their “parent” nondifferential matrices are known as *M*-matrices. The differential *N*-matrices corresponding to each optical property exhibited by the sample are then summed to express the combined effect. The formalism is expressed in terms of 2 × 2 Jones matrices applicable to nondepolarizing media, rather than the more commonly used 4 × 4 Mueller matrices. For example, the *N* matrix for combined linear birefringence and optical activity effects is given by^{61}

*g*

_{0}= 2π/λΔ

*n*is the phase retardation per unit distance and χ is the optical rotation per unit distance. The “parent” nondifferential

*M*-matrix is calculated from the

*N*matrix to describe the combined effect.

^{61}The resulting Jones

*M*-matrix is then converted to a Mueller matrix.

^{21}Note, however, a Jones matrix description, and thus its conversion to a Mueller matrix, is only valid provided there are no depolarization effects.

^{21}This is indeed applicable in the Monte Carlo model, as depolarization is predominantly caused by the multiple scattering, and no depolarization effects should occur between the scattering events. Once converted to a Mueller matrix, this matrix is then applied to the photons as they propagate between scattering events. This approach thus enables the combination of any number of simultaneously occurring polarizing effects. In tissues, linear birefringence and circular birefringence (optical activity) are the important polarimetry characteristics to be added to the multiple scattering effects.

Similar to conventional Monte Carlo modeling, in PSMC also, the scattering and absorption properties are modeled using the optical transport parameters, scattering coefficient (*μ*_{s}) and absorption coefficient (*μ*_{a}). Mie theory is used to compute the single scattering Mueller matrix for known diameter (*D*) and refractive index of scatterer (*n*_{s}) and refractive index of the surrounding medium (*n*_{m}). Circular and linear birefringence is modeled through the optical activity χ in degrees per centimeter, and through the linear anisotropy in refractive indices Δ*n* (= *n*_{e}−*n*_{o}), respectively. For simplicity, it is generally assumed that the medium is uniaxial and that the direction of the extraordinary axis and the value for Δ*n* is constant throughout the scattering medium^{61} (although our recent research efforts are exploring the effects of multiple uniaxial domains of varying magnitude and orientation of birefringence). Note that as photons propagate between scattering events, the difference in refractive indices [*n*(ϑ)−*n*_{o}] experienced by them depends on their propagation direction with respect to the extraordinary axis (ϑ). The effect is modeled using standard formulae describing the angular variation of refractive index in uniaxial medium,

The validity of this model has been tested on experimental tissue simulating phantoms exhibiting simultaneous scattering and polarization properties, which are known and user-controlled *a priori*.^{15- 16,37} These solid optical phantoms were developed using polyacrylamide as a base medium, with sucrose-induced optical activity, polystyrene microspheres-induced scattering, and mechanical stretching to cause linear birefringence. These phantom system mimics the complexity of biological tissues, in that it exhibits simultaneous linear birefringence, optical activity, and depolarization due to multiple scattering.^{61} Figure 4 shows the change in the normalized Stokes parameter *q* = *Q/I* with increasing birefringence, measured in phantoms and calculated from the PSMC model in the forward direction of a 1 ×1 ×1 cm^{3} sample with input circularly polarized light.^{61} The measurements were performed using the experimental system shown in Fig. 1. Good agreement between the Monte Carlo model and controlled experimental results is seen. As the input light is transferred from circular to linear polarization due to the increasing sample linear birefringence (the sample in effect acting like a turbid wave-plate), optical rotation due to optical activity of dissolved sucrose is seen as an increase in parameter *q*. No such effect is seen in the absence of chirality.

**F4 :**

Experimental measurements (symbols) and Monte Carlo calculations (lines) of the change in the normalized Stokes parameter *q* with and without optical activity (with and without added sucrose) in the forward (γ = 0 deg) detection geometry with input circularly polarized light and a fixed scattering coefficient of μ_{s} = 60 cm^{−1}. Linear retardance was varied from *δ* = 0 to 1.4364 rad (corresponding to birefringence variation *Δn* = 0 to 1.628 × 10^{−5}) and the magnitude of optical activity was *χ* = 1.965 deg cm^{−1}, corresponding to a 1 M sucrose concentration. (Adopted from Refs. ^{16,61}.)

Figure 5 gives an example of the 4 × 4 Mueller matrix experimentally recorded in the forward (transmission) detection geometry from a birefringent (extension = 4 mm for strain applied along the vertical direction, corresponding to a value of linear retardance δ = 1.345 rad for a clear phantom of thickness of 1 cm), chiral (optical activity χ = 1.96 degree cm^{−1}, corresponding to 1 M concentration of sucrose), turbid phantom (*μ*_{s} = 30 cm^{−1} and *g* = 0.95).^{37} The corresponding matrix generated through the PSMC model, using the same controlled input scattering and polarization parameters (linear birefringence Δ*n* = 1.36 ×10^{−5}, corresponding to δ = 1.345 rad, χ = 1.96 degree cm^{−1}, μ_{s} = 30 cm^{−1}, *g* = 0.95) is shown in Fig. 5.^{37} In both the experimental and the MC-generated Mueller matrices, the signature of linear retardance is prominent in the matrix elements *M*_{34} and *M*_{43} [as is expected for a retarder having an orientation angle θ = 90 deg, see Eq. 10]. The effect of optical activity is mainly manifest as a difference in *M*_{23} and *M*_{32} elements [see Eq. 11], whereas depolarization effects are most prominently reflected in the diagonal elements of the Mueller matrix. The excellent agreement between the experimental and the simulated Mueller matrices emphasizes the capability of the PSMC model in simulating complex tissue polarimetry effects, including simultaneous optical activity and birefringence in the presence of multiple scattering in any desired detection geometry. However, the complex nature of the recorded Mueller matrix **M**, with essentially all 16 nonzero elements also underscores a significant problem—how does one extract useful sample metrics from this intertwined array of information? Inverse analysis aimed at quantifying individual polarimetry contributions from “lumped” Mueller matrix are described subsequently (Sec. 6). In Sec. 5, we discuss an interesting application of turbid polarimetry for polarization gated imaging of tissue.

**F5 :**

The 4 × 4 Mueller matrices (a) experimentally recorded in the forward (transmission) detection geometry from a birefringent (extension = 4 mm, corresponding to a value of linear retardance δ = 1.345 rad for a clear phantom), chiral concentration of sucrose = 1 M, χ = 1.96 deg/cm) turbid (μ_{s} = 30 cm^{−1}, *g* = 0.95) phantom of thickness of 1 cm; (b) generated through the PSMC model, using the same controlled input scattering and polarization parameters (linear birefringence Δ*n* = 1.36 ×10^{−5}, corresponding to δ = 1.345 rad, χ = 1.96 deg cm^{−1}, μ_{s} = 30 cm^{−1}, *g* = 0.95) as that of the experimental phantom. The results of the decomposition analysis on the experimental Mueller matrix are presented in Fig. 7 and Table 1.

## Polarized Light Scattering as a Gating Mechanism for Tissue Imaging and Spectroscopy

Investigation of optical techniques for noninvasive, high resolution imaging of tissue and its underlying structure is an actively pursued area.^{66- 67} However, a major difficulty in realizing the potential of the optical techniques for medical imaging is the fact that multiple scattering in tissue leads to loss of directionality that causes severe image blurring. Polarization gating, along with other approaches such as time gating, spatial filtering, and coherence gating have shown some promise in isolating the image bearing photons [the unscattered (ballistic) or weakly scattered (snake) photons] from the image blurring (diffusive or multiply scattered) photons for Ref. ^{66}. Among these attempts, polarization gating has received particular attention because of the relative simplicity of instrumentation, making this a potential tool amenable for clinical use.^{68} The underlying principle of this approach is that multiply scattered light gets depolarized, and therefore by detecting primarily the polarization-retaining component of light scattered from tissue, one can filter out the multiply scattered photons and thus improve imaging resolution. However, for practical implementation of this approach, it is essential to understand the mechanism of depolarization and its dependence on the different morphological parameters of tissue (e.g., the density, size and its distribution, shape and refractive index of the tissue scatterers). A number of experimental studies have therefore investigated these issues.^{26- 30,68- 71} We shall briefly summarize these studies (with illustrative examples) and discuss their implications for polarimetric tissue imaging.

###### Depolarization of Light in Tissues and Tissue Phantoms: Implications for Polarized Light Tissue Imaging

Since it is difficult to control and quantify relevant properties in a complex system like tissue, the usual approach followed in most experimental depolarization investigations has been to study the influence of the individual morphological parameters (user controlled and known *a priori*) in well-characterized tissue simulating phantoms, and relate these to depolarization behavior of light in actual tissue. Suspensions of Intralipid and aqueous suspension of polystyrene (or silica) microspheres have commonly been used to prepare tissue phantoms for these depolarization studies.^{26- 30} Usually the size of the scatterers is chosen such that the value of anisotropy parameter (*g*) is comparable to that of tissue. The concentration of scatterers is then adjusted to yield the value of the scattering coefficient μ_{s} or the reduced scattering coefficient μ_{s}^{′} desired to mimic the tissue of interest. Molecular dyes of know molar extinction are used to simulate tissue absorption.

For scattering media comprised of Rayleigh scatterers (radius of scatterer *a* ≪ λ, size parameter *X* < 2 and anisotropy parameter *g* ≤ 0.2), in the forward scattering geometry, depolarization of incident circular polarization is stronger compared to the depolarization of incident linear polarization (as predicted by heuristic models described in Sec. 4a). The reverse is the case for scattering media comprised of large sized Mie scatterers (*a* ≥ λ, *X* >2, *g* ≥ 0.7). Moreover, the overall strength of depolarization is considerably weaker in the latter case. Typical results for such experimental depolarization measurements are shown in Fig. 6, where the degree of polarization for transmitted light is shown as a function of optical thickness τ (= μ_{s} × *L*, *μ*_{s} is the scattering coefficient and *L* is the physical thickness) from samples prepared using aqueous (refractive index *n*_{m} = 1.33) suspensions of polystyrene microspheres (refractive index *n*_{s} = 1.59) having mean diameter of 0.11 μm (Rayleigh scatterers, *g* = 0.09, *X* = 0.72 at 632.8 nm) and 1.08 μm (Mie scatterers, *g* = 0.92 and *X* = 7.13).^{28,30} Again, these results are in qualitative agreement with the predictions of the corresponding theoretical models for polarized light transport in turbid medium (described in Sec. 4a). The weaker depolarization for media comprised of large sized scatterers can be explained by noting that with increasing value of size parameter of scatterer *X*, forward scattering becomes more predominant (value of *g* increases), resulting in weaker randomization of photon's propagation direction when forward (transmission) detection geometry is employed.^{53- 56,68} The difference in relative rate of depolarization of linearly and circularly polarized light can be attributed to the different mechanism of depolarization of the two. While randomization of the incident field vector's direction as a consequence of multiple scattering events is solely responsible for the depolarization of incident linear polarization state, the depolarization of circular polarization state is additionally influenced by the randomization of the helicity (handedness).^{54} Note that scattering at large angles flips the helicity of the circular polarization state. For Rayleigh scatterers, the critical scattering angle, above which the helicity of circularly polarized photons is reversed, is 90 deg. However, this helicity flip angle increases (>90 deg) with increasing size of the scatterers.^{28,68} In a turbid medium, light travels along many possible zig-zag paths, having contributions from scattering at various angles. For Rayleigh scatterers (where forward and back scattering is ∼ equally likely), the contribution of the large angle scattered photons is greater as compared to the larger sized Mie scatterers (where forward scattering predominates), ensuing stronger randomization of helicity and thus resulting in stronger depolarization of circularly polarized light in Rayleigh media (*X* < 2, *g* ≤ 0.2). As the scatterer size increases, the additional cause of depolarization of circularly polarized light, the flipping of helicity also gets considerably reduced (due to a relatively lesser contribution of large angle scattered photons coupled with the fact that the helicity flip angle is also larger), resulting in weaker depolarization of circularly polarized light for anisotropic (*g* ≥ 0.7) scattering media comprised of Mie scatterers (*X* >2). It may also be pertinent to note that the corresponding theoretical variation of the length scales of depolarization (ξ_{l} and ξ_{c}) showed a general decreasing trend beyond a size parameter value of *X* ≥ 10 (see Fig. 2, Sec. 4a). This anomalous depolarization behavior has been attributed to the Mie resonance effects present in the large sized scatterers.^{53}

**F6 :**

(a) Measured degrees of linear (DOLP, triangles) and circular (DOCP, circles) polarization as a function of optical thickness τ (= μ_{s} × *L*, μ_{s} is the scattering coefficient and *L* is the thickness) from samples prepared using aqueous (refractive index *n*_{m} = 1.33) suspension of polystyrene microspheres (refractive index *n*_{s} = 1.59) having mean diameter of 0.11 μm (Rayleigh scatterers, *g* = 0.09, *X* = 0.72 at 632.8 nm) and 1.08 μm (Mie scatterers, *g* = 0.92 and *X* = 7.13) as scatterers. The measurements were performed in transmission geometry. The thickness of the sample was 1 cm and the optical thickness τ was varied by changing the value of scattering coefficient *μ*_{s} through dilution. The symbols are experimental data and the lines are guide for the eye. (Adopted from Ref. ^{28}.) (b) Measured degrees of linear (triangles) and circular polarization (circles) of optical thickness τ from samples prepared using aqueous suspension of 0.65-μm diameter polystyrene (*n*_{s} = 1.59, *m* = 1.03, *g* = 0.86, *X* = 4.29) and silica microspheres (*n*_{s} = 1.37, *m* = 1.2, g = 0.89, *X* = 4.29). The measurements were performed in transmission geometry. The thickness of the sample was 1 cm and the optical thickness τ was varied by changing the value of scattering coefficient *μ*_{s} through dilution. The symbols are experimental data and the lines are guide for the eye. (Adopted from Ref. ^{29}.)

The depolarization is also influenced by the refractive index contrast present in the turbid medium.^{28,55} Depolarization characteristics of Mie-like anisotropic scattering media (*g* ≥ 0.7) comprised of low refractive index scatterers [e.g., silica in water, low value of the relative refractive index *m* (= *n*_{s}/*n*_{m}) ∼ 1.03 to 1.05] is distinctly different from those comprised of high refractive index scatterers (*m* ≫ 1.05, e.g., polystyrene in water). In fact, the depolarization behavior in low-index-contrast Mie medium surprisingly resembles the behavior for Rayleigh scatterers.^{28} This is illustrated in Fig. 6, where the measured degrees of linear and circular polarization from samples prepared using aqueous suspension of 0.65-μm diameter polystyrene (*n*_{s} = 1.59, *m* = 1.2, *g* = 0.86, *X* = 4.29) and silica microspheres (*n*_{s} = 1.37, *m* = 1.03, *g* = 0.89, *X* = 4.29) are displayed.^{29} Despite having a slightly larger value for *g*, depolarization of both linearly and circularly polarized light is stronger for the samples having silica microspheres as compared to those having polystyrene microspheres as scatterers. Moreover, for the silica microspheres samples, circular polarization actually decays more rapidly as compared to linear polarization, opposite of the Mie trends previously described and in fact more associated with Rayleigh scatterers exhibiting small *g* (≤ 0.2). This interesting behavior originates from the fact that the anisotropic scatterers (*g* ≥ 0.7, *X* ≫ 2) having a lower value of relative refractive index *m* [*m* ≤ 1.05, (*m*−1) X ≪ 1] can be treated to be in the weak scattering regime, where each volume element inside the scatterer can be assumed to be giving electrical dipole scattering in an independent manner.^{29,56} The corresponding scattering matrix becomes similar in nature to Rayleigh scatterers (ideal dipole scattering). The retention of this dipolar behavior thus leads to depolarization characteristics of such weakly fluctuating anisotropic scattering medium much similar to that of Rayleigh scatterers. Since biological tissue can be considered as a random continuum of inhomogenities having weak fluctuations in refractive indices (refractive indices of scattering inhomogenity *n*_{s} ∼ 1.37 to 1.41, of background medium *n*_{m} ∼ 1.33 to 1.35, thus *m* ≤ 1.05), its depolarization behavior is thus expected to be similar to that of Rayleigh scatterers.^{56} This indeed has been observed in experimental depolarization studies on actual tissues,^{26- 27} despite the potentially mistaken paradigm that “circular is better in Mie-scattering.”

Once again, caution must be exercised in generalizing the depolarization behavior in a turbid medium based on the experimental results described above. As noted previously, the quantitative rate of depolarization and the relative depolarization behavior of linearly and circularly polarized light may also be additionally influenced by detection direction/geometry. In order to illustrate this point, let us consider the following: incident circularly polarized photons getting multiply scattered in an anisotropic (*g* ≥ 0.7) scattering medium, with a photodetector positioned for measurements in the backward direction. At detection points sufficiently close to the point of illumination, a major component of the detected photons are the helicity flipped backscattered photons. In contrast, as one moves away from the point of illumination, there would be significant contribution of the sub-population of helicity preserving photons which have undergone a series of forward scattering events to emerge through the backward direction. The observed depolarization at a chosen detection position in the backward direction will thus be the resultant of a complex combination of these many possible zig-zag photon paths, which may differ from the experimental trends presented above.

Finally, the efficiency of the polarization gating scheme for filtering out the multiply scattered light is determined by the degree to which incident polarized light becomes randomized or depolarized. It is thus desirable that the multiply scattered photons (whose direction information has been randomized) should also be depolarized to a greater extent. These photons can then be efficiently filtered out by imposing a polarization gate that will preferentially detect the polarization-retaining component of light emerging from the medium. Based on the results of the depolarization studies discussed above, the following inferences can be made: i. For scattering media comprised of Rayleigh scatterers, polarization (specifically circular polarization) can be used as an effective scheme to filter out multiply scattered photons. ii. For anisotropic scattering media comprised of high refractive index, large-sized Mie scatterers, the polarization filtering schemes will be less efficient. iii. As in the case of Rayleigh scattering, for anisotropic scattering media comprised of low refractive index, large-sized scatterers (or weakly fluctuating random media), polarization gating should be an efficient method for discriminating ballistic and snake from diffusive photons. Since in terms of the scattering properties, biological tissue falls in category iii., polarization-gated methods may be advantageous compared to corresponding polarization-blind methods for tissue imaging. Several research groups have therefore actively pursued various polarization gating schemes and strategies for imaging of tissue and its underlying structures, as summarized below.

The simplest form of polarized light imaging method employs illumination of tissue by either linearly or circularly polarized light, and recording of images with either the same orientation of the polarizer and the analyzer (co-polarized, *I*_{co}) or with the opposite orientation (cross-polarized, *I*_{cross}). A similar approach for linear polarization illumination has already become common practice in dermatology.^{72} Viewing the skin with co- (or crossed-) polarization channels allows the dermatologist to either observe the glare from the air/skin surface and surface details (or to suppress it for a better view the deeper subsurface tissue structures). Jacques et al. ^{17,73- 74} modified this scheme further for high resolution imaging of the texture of the superficial (sub-surface) skin structures. In this approach, the glare from the air/glass/skin interfaces was removed through optical coupling (using index matching liquid) of a glass plate in contact with the skin and by off-normal illumination. The multiply scattered photons were filtered out by extracting the polarization preserving component of the scattered light employing either the polarization difference scheme (*I*_{co} – *I*_{cross}) or the degree of polarization scheme [DOP = (*I*_{co} − *I*_{cross})/(*I*_{co} + *I*_{cross})]. In fact, the DOP imaging has subtle advantages over the polarization difference imaging in that the DOP imaging is relatively insensitive to spatial variation of illuminating light and variations in pigmentation. This follows because the DOP image is formed by the ratio of the numerator (primarily comprised of superficial subsurface reflectance) and the denominator (representative of total subsurface reflectance). Based on this simple principle, a polarized light camera system has been developed and used in clinical studies for finding skin cancer margins.^{17}

Morgan et al. ^{75} developed another elegant polarization scheme for eliminating surface reflection in backscattering imaging of superficial tissue layers, obviating the use of index optically flat plates (and index matching fluid) and angled illumination (thus enabling co-axial detection). It was shown that a suitable combination of co- and cross-polarized images acquired using linear and circular polarization illumination can be used to simultaneously filter out both the multiply scattered light and the surface reflection. The relative efficacies of the different polarization gating schemes for tissue imaging have also been investigated in detail by this group.^{69}

Demos et al. ^{76} have developed a methodology based on both spectral and polarization discrimination of backscattered photons for deep subsurface imaging of tissue. This approach, known as spectral polarization difference imaging, exploits the difference in light penetration of different wavelengths in combination with polarization filtering for selectively imaging different tissue depths.

Orthogonal polarization spectral imaging is another important development in polarimetric tissue imaging.^{77} In this approach, the tissue is illuminated with linearly polarized light and the backscattered light is detected in cross-polarization channel. Since the detected depolarized scattered light (passing through orthogonal polarizer), penetrates deeper in tissue, it effectively back-illuminates the absorbing material (e.g., blood present in tissue) in the foreground. Appropriate choice of illuminating wavelengths (typically in the range of 550 nm, optimized based on the absorption of the hemoglobins present in blood) enables one to image blood vessels in tissue. This approach has therefore been exploited for quantitative assessment of the blood vessels in microcirculation for a number of diagnostic applications.^{78}

Apart from the few aforementioned examples, several other experimental approaches such as rotating polarizer imaging,^{79} complete 4 × 4 (16 element) Mueller matrix imaging,^{80} and others have also been explored for polarimetric tissue imaging.

Investigation of the optical spectroscopic techniques (such as fluorescence, Raman, and elastic scattering spectroscopy) for biomedical diagnosis is an area of considerable current interest, as optical methods can facilitate noninvasive quantitative diagnosis and can provide (functional/biochemical/compositional/metabolic) information beyond tissue structure. Indeed, very encouraging results have already been obtained in using these approaches for early detection