Mueller matrix polarimetry for improved liver fibrosis diagnosis

An experimental Mueller matrix polarimeter is used to quantify human liver fibrosis by measuring retardance and depolarization of thin biopsies. The former parameter is sensitive to fibrillar collagen, the latter is specifically sensitive to fibrillar collagen around blood vessels, which is not significant for liver fibrosis diagnosis. By using depolarization like a filter, retardance distribution enables distinguishing between disease stages and limits the high degree of observer discrepancy.     

Double-pass axially resolved confocal Mueller matrix imaging polarimetry

We present experimental depth-resolved complete polarization-sensitive measurements of a stack of linear retarders and glass plates by using what is to the best of our knowledge the first combination of a confocal imaging system with a complete Mueller matrix polarimeter. The axially resolved Mueller matrices were compared with a forward simulation, with good agreement.     

Two-dimensional depth-resolved Mueller matrix characterization of biological tissue by optical coherence tomography

We built a polarization-sensitive optical coherence tomographic system and measured the two-dimensional depth-resolved full 4 x 4 Mueller matrix of biological tissue for what is believed to be the first time. The Mueller matrix measurements, which we made by varying the polarization states of the light source and the detector, yielded a complete characterization of the polarization property of the tissue sample. The initial experimental results indicated that this new approach reveals some tissue structures that are not perceptible in standard optical coherence tomography.     

Optical diagnosis of dengue virus infected human blood using Mueller matrix polarimetry

Currently dengue fever diagnosis methods include capture ELISAs, immunofluorescence tests, and hemagglutination assays. In this study optical diagnosis of dengue virus infection in the whole blood is presented utilizing Mueller matrix polarimetry. Mueller matrices of about 50 dengue viral infected and 25 non-dengue healthy blood samples were recorded utilizing light source from 500 to 700 nm with scanning step of 10 nm. Polar decomposition of the Mueller matrices for all the blood samples was performed that yielded polarization properties including depolarization, diattenuation, degree of polarization, retardance and optical activity, out of which, depolarization index clusters up the diseased and healthy in to different separate groups. The average depolarized light in the case of dengue infection in the whole blood at 500 nm is 18%, whereas for the healthy blood samples it is 13.5%. This suggests that depolarization index of polarized light at the wavelengths of 500, 510, 520, 530 and 540 nm, we find that in case of depolarization index values are higher for dengue viral infection as compared to normal samples. This technique can effectively be used for the characterization of the dengue virus infected at an early stage of disease.     

Influence of the order of the constituent basis matrices on the Mueller matrix decomposition-derived polarization parameters in complex turbid media such as biological tissues

The influence of the multiplication order of the constituent basis matrices on the Mueller matrix decomposition-derived polarization parameters in complex tissue-like turbid media exhibiting simultaneous scattering and polarization effects are investigated. A polarization sensitive Monte Carlo (MC) simulation model was used to generate Mueller matrices from turbid media exhibiting simultaneous linear birefringence, optical activity and multiple scattering effects. Mueller matrix decomposition was performed with different selected multiplication orders of the constituent basis matrices, which were further analyzed to derive quantitative individual polarization medium properties. The results show that for turbid medium having weak diattenuation (differential attenuation of two orthogonal polarization states), the decomposition-derived polarization parameters are independent of the multiplication order. Importantly, the values for the extracted polarization parameters were found to be in excellent agreement with the controlled inputs, showing self-consistency in inverse decomposition analysis and successful decoupling of the individual polarization effects. These results were corroborated further by selected experimental results from phantoms having optical (scattering and polarization) properties similar to those used in the MC model. Results from tissue polarimetry confirm that the magnitude of diattenuation is generally lower compared to other polarization effects, so that the demonstrated self-consistency of the decomposition formalism with respect to the potential ambiguity of ordering of the constituent matrices should hold in biological applications.     

Error in determining Mueller matrix elements and its effect on solution of the inverse problem of polarimetry

We present an analysis of the errors in measurement of the Mueller matrix elements in polarimeters that are a combination of different types of Stokes polarimeters in the detection channel and controllable polarization converters in the probing channel. As polarization converters for the probing radiation, we consider a phase plate having four different angular positions in measuring the complete Mueller matrix and a linear polarizer having different angular positions in measuring the structural parts of an incomplete 4 × 3 Mueller matrix. We have shown that the error in determining the Mueller matrix elements is distributed nonuniformly over the matrix. The nature of the error distribution over the elements and its values are different for different combinations of detection and probing channels of the polarimeter, and depend on the anisotropy of the test object. The latter dictates the choice of the optimal layout for a Mueller polarimeter for studying media with the same or different types of anisotropy and the choice of Mueller matrix elements used to solve the inverse problem of determining the anisotropy parameters.     

Differentiating characteristic microstructural features of cancerous tissues using Mueller matrix microscope

Polarized light imaging can provide rich microstructural information of samples, and has been applied to the detections of various abnormal tissues. In this paper, we report a polarized light microscope based on Mueller matrix imaging by adding the polarization state generator and analyzer (PSG and PSA) to a commercial transmission optical microscope. The maximum errors for the absolute values of Mueller matrix elements are reduced to 0.01 after calibration. This Mueller matrix microscope has been used to examine human cervical and liver cancerous tissues with fibrosis. Images of the transformed Mueller matrix parameters provide quantitative assessment on the characteristic features of the pathological tissues. Contrast mechanism of the experimental results are backed up by Monte Carlo simulations based on the sphere–cylinder birefringence model, which reveal the relationship between the pathological features in the cancerous tissues at the cellular level and the polarization parameters. Both the experimental and simulated data indicate that the microscopic transformed Mueller matrix parameters can distinguish the breaking down of birefringent normal tissues for cervical cancer, or the formation of birefringent surrounding structures accompanying the inflammatory reaction for liver cancer. With its simple structure, fast measurement and high precision, polarized light microscope based on Mueller matrix shows a good diagnosis application prospect.     

Measurement and calculation of the two-dimensional backscattering Mueller matrix of a turbid medium

We present both experimental and Monte Carlo-based simulation results for the diffusely backscattered intensity patterns that arise from illumination of a turbid medium with a polarized laser beam. A numerical method that allows the calculation of all 16 elements of the two-dimensional Muller matrix is used; moreover, it is shown that only seven matrix elements are independent. To validate our method, we compared our simulations with experimental measurements, using a turbid medium consisting of 2.02-microm -diameter polystyrene spheres suspended in deionized water. By varying the incident polarization and the analyzer optics for the experimental measurements, we obtained the diffuse backscattering Mueller matrix elements. The experimental and the numerical results are in good agreement.     

随机非球形粒子全极化散射的时间相关Mueller矩阵解

从与时间相关的矢量辐射传输方程推导一阶Mueller矩阵解,用来模拟Gauss型平面脉冲波入射下,一层随机、非均匀取向非球形粒子的全极化双站散射.数值计算了同极化和去极化脉冲响应,与入射脉冲进行了比较,说明了随机介质的物理参数,如粒子的取向和占空比、入射角、极化以及层厚等对脉冲响应的影响 关键词： 平面脉冲波 非球形粒子 Mueller矩阵     

The temporal Mueller matrix solution for polarimetric scattering from a layer of random non-spherical scatterers

The Mueller matrix solution of the vector radiative transfer equation with time dependence is derived in this paper. It is applied to simulation of polarimetric bistatic scattering from a layer of non-uniformly oriented, random non-spherical scatterers when a Gaussian plane pulse is incident upon. Co-polarized and cross-polarized bistatic scattering are numerically calculated. The pulse echoes are compared with the incidence, and demonstrate its functional dependence on the physical parameters of random medium, such as spatial orientation and fractional volume of scatterers, incidence angle and polarization, the layer depth and others.     

Optical polarimetry of random fields

The state of polarization of an optical field provides detailed information concerning both the radiation emission processes and the intricate interaction between light and matter. We report here a novel approach for characterizing the polarization properties of electromagnetic fields for which the electric field vector at a point may fluctuate in three dimensions. Using probes which couple all three components of the field, we were able to extract the polarized and unpolarized components of such fields. Our results constitute the proof of concept for what could be called three-dimensional optical polarimetry.     

Least-squares analysis of the Mueller matrix

In a single-mode fiber excited by light with a fixed polarization state, the output polarizations obtained at two different optical frequencies are related by a Mueller matrix. We examine least-squares procedures for estimating this matrix from repeated measurements of the output Stokes vector for a random set of input polarization states. We then apply these methods to the determination of polarization mode dispersion and polarization-dependent loss in an optical fiber. We find that a relatively simple formalism leads to results that are comparable with those of far more involved techniques.     

Measurement and calculation of the two-dimensional backscattering Mueller matrix of a turbid medium: errata

In our recent Letter,(1) several typographical errors were present. On p. 487, in Fig. 2, the equations for the following Mueller matrix elements should read as S(14) = (RO - LO), S(22) = (HH + VV) - (HV + VH), S(23) = (PH + MV) - (PV + MH), S(24) = (RH + LV) - (RV + LH), S(32) = (HP + VM) - (HM + VP), S(33) = (PP + MM) - (PM + MP), S(34) = (RP + LM) - (RM + LP), S(41) = (OR + OL), S(42) = (HR + VL) - (HL + VR), S(43) = (PR + ML) - (PL + MR), and S(44) = (RR + LL) - (RL + LR). Also on p. 487, in the left-hand column, line 10 from the top should read as follows: mfp? = 1/[mua + mus(1 - g)], was 0.957 cm.     

Mueller matrix for characterization of one-dimensional rough perfectly reflecting surfaces in a conical configuration

Theoretical results of the use of a Mueller matrix to characterize a one-dimensional rough perfectly reflecting, single-scattering surface in a conical configuration are presented. The conical Mueller matrix (CMM) is derived from the known Mueller matrix of this kind of surface in the plane of incidence [the plane Mueller matrix (PMM)]. The key argument is that, as the PMM is considered to be a Mueller-Jones matrix, an appropriate rotation of the complex amplitude matrix provides the conic Mueller matrix.     

Polarimetry group theory analysis in biological tissue phantoms by Mueller coherency matrix

The characterization of biological tissues by optical techniques provides several advantages over other techniques. Optical techniques enable to perform high resolution and contrast imaging, in a non-invasive way and with no-contact. Biological tissues are turbid media that strongly scatter light. The ultrastructure of some tissues makes them present a certain degree of anisotropy. Both scattering and anisotropy affect light polarization. Some pathologies alter these characteristics of the tissue. As a consequence polarized light can be used to extract additional information and achieve a better diagnosis.In this work, Group Theory is applied to analyse the polarization behavior of several samples. Firstly, the Mueller matrix for each sample is measured. Then, the Mueller Coherency matrix is obtained by means of the SU(4)-O + (6) homomorphism. Finally, the target decomposition theorem is applied by analyzing the eigenvalues and eigenvectors, and subsequently the different polarimetric effects are separated. In this way, the contrast of tissue imaging can be increased. This analysis is applied to biological tissue phantoms, which consisted on glucose suspensions of polystyrene spheres with different scatterer concentrations. Their behaviour can be modeled by means of single or multiple scattering depending on the concentration, either in the Rayleigh or Mie regimes. The same procedure could be used in a wide range of applications, like the study of cancerous cells that grow without control in cell cultures, or erythrocytes monitoring in anemia. The technique also has a great potential to be applied in Polarization Sensitive Optical Coherence Tomography (PS-OCT).     

Combined Mueller and Jones matrix method for the evaluation of the complex modulation in a liquid-crystal-on-silicon display

We apply the polar decomposition of the Mueller matrix describing a liquid-crystal-on-silicon display to identify the diattenuator, depolarizer, and retarder contributions as a function of the gray level. The retarder contribution is expressed in terms of the equivalent Jones matrix to apply previous techniques to evaluate the phase modulation. This allows searching for optimized polarization configurations for phase- or amplitude-only modulation responses. We present results for lambda=633 nm showing a phase-only modulation up to 2 pi rad and flat intensity modulation.     

Nonlinear integro-differential models for intense waves in media like biological tissues and geostructures with complex internal relaxation-type dynamics

The paper discusses a universal scheme for constructing nonlinear integro-differential models to describe intense waves in media with a complex internal relaxation-type dynamics. Examples of such media are presented. Various forms of kernels are described. Situations are shown in which the models can be simplified by reducing them to differential or differential-difference equations with partial derivatives. Integral relations for the linear momentum and energy transferred by the wave are obtained. Exact solutions are found. The mapping method is used to obtain approximate solutions and analyze them in the form of difference relations.     

Characteristic polynomials of complex random matrix models

We calculate the expectation value of an arbitrary product of characteristic polynomials of complex random matrices and their Hermitian conjugates. Using the technique of orthogonal polynomials in the complex plane our result can be written in terms of a determinant containing these polynomials and their kernel. It generalizes the known expression for Hermitian matrices and it also provides a generalization of the Christoffel formula to the complex plane. The derivation we present holds for complex matrix models with a general weight function at finite-N, where N is the size of the matrix. We give some explicit examples at finite-N for specific weight functions. The characteristic polynomials in the large-N limit at weak and strong non-hermiticity follow easily and they are universal in the weak limit. We also comment on the issue of the BMN large-N limit.     

Characteristic properties of the conical Mueller matrix

We present results of mathematical relations existing between the Mueller matrix obtained for an in-plane of incidence scattering geometry (plane Mueller matrix, PMM) and the Mueller matrix obtained for an out-the-plane of incidence scattering geometry (conical Mueller matrix, CMM), for light scattered from a rough surface. We obtain a similarity relation between the CMM and the PMM for one- (1-D) and for two-dimensional (2-D) surfaces. This similarity relation implies that the PMM and the CMM have the same determinant, trace and eigenvalues for 1-D and 2-D surfaces, respectively. We can say that measurements made in the conical geometry are “Similarity Equivalent” to those in the in-plane geometry for both kind of surfaces, respectively.     

Stokes-Mueller matrix polarimetry system for glucose sensing

A Stokes-Mueller matrix polarimetry system consisting of a polarization scanning generator (PSG) and a high-accuracy Stokes polarimeter is used to sense the glucose concentration in aqueous solutions with and without scattering effects, respectively. In the proposed system, an electro-optic (EO) modulator driven by a saw-tooth waveform voltage is used to perform full state of polarization (linear/circular) scanning, while a self-built Stokes polarimeter is used to obtain dynamic measurements of the output polarized light intensity. It is shown that the measured output Stokes vectors have an accuracy of 10⁻⁴, i.e., one order higher than that of existing commercial Stokes polarimeters. The experimental results show that the optical rotation angle varies linearly with the glucose concentration over the range of 0-0.5 g/dl. Moreover, glucose sensing is successfully achieved at concentrations as low as 0.02 g/dl with a resolution of 10⁻⁶ deg/mm and an average deviation of 10⁻⁴ deg. In general, the polarimetry system proposed in this study provides a fast and reliable method for measuring the Stokes vectors, and thus has significant potential for biological sensing applications.