Numerical study of the characterization of forward scattering Mueller matrix patterns of turbid media with triple forward scattering assumption

We propose a triple forward scattering model to numerically investigate the forward scattering Mueller matrix of a turbid medium. The calculated results show that the Mueller matrix based on the triple scattering model can completely describe the basic symmetries and azimuthal structures of the forward scattering Mueller matrix of a turbid medium. The results show that the forward scattering Mueller matrix is characterized with special symmetric structure compared with backscattering Mueller matrix patterns. The method will extend the investigation to light scattering mechanism from cells and to diagnosis of diabetes and other blood related diseases.     

Influence of single scattering and multiple scattering on backscattered Mueller matrix in turbid media

Monte Carlo algorithm and Stokes-Mueller formalism are used to simulate the propagation behavior of polarized light in turbid media. The influence of single scattering and multiple scattering on backscattered Mueller matrix in turbid media is discussed. Single and double scattering photons form the major part of backscattered polarization patterns, while multiple scattering photons present more likely as background. Further quantitative analyses show that single scattering approximation and double scattering approximation are quite accurate when discussing the polarization patterns near the incident point.     

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.     

随机介质背散射二维Mueller矩阵的数值模拟与分析

推导了随机介质背散射Mueller矩阵的直接计算公式,并运用矢量Monte Carlo方法进行了数值模拟.结果表明随机介质背散射二维Mueller矩阵方位关系随散射系数的减小而增强,而与微粒大小关系不大;Mueller矩阵元素绝对值的空间分布随径向呈近似指数规律衰减,矩阵元素的方位变化具有周期性.对称系统的二维Mueller矩阵的花样图中仅有7幅独立,其余9幅可通过对称、旋转变换得到.     

混浊介质中后向单次漫散射米勒矩阵特征

采用斯托克斯(Stokes)矢量形式,推导出当无限窄的连续光束垂直入射到混浊介质表面时,后向单次漫散射米勒(Mueller)矩阵的解析表达式。基于米氏(Mie)散射模式,详细分析了单次散射米勒矩阵元素的分布模式,以及与介质粒子数密度,粒子尺寸参量之间的关系。研究表明:单次散射米勒矩阵的方位变化随粒子数密度的增加,逐渐消失,而矩阵元素m22,m33,m23,m32随粒子数密度的变化,具有更显著的方位变化特征。矩阵元素m22,m33在方位角=45°时的值随尺寸参量的变化有一定的规律性,当尺寸参量小于某一特征参量时,其值呈下降趋势,反之则呈波动上升趋势。当介质粒子数密度以及粒子尺寸参量改变时,米勒矩阵元素强度的径向分布模式不变,即在任何方位,强度随径向距离都近似成指数规律衰减,方位变化呈周期性。     

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 matrix polarimetry for the characterization of complex random medium like biological tissues

The polarization parameters of light scattered from biological tissues contain wealth of morphological and functional information of potential biomedical importance. But, in optically thick turbid media such as tissues, numerous complexities due to multiple scattering and simultaneous occurrences of many polarization events present formidable challenges, in terms of both accurate measurement and unique interpretation of the individual polarimetry characteristics. We have developed and validated an expanded Mueller matrix decomposition approach to overcome this problem. The approach was validated theoretically with a polarization-sensitive Monte Carlo light propagation model and experimentally by recording Mueller matrices from tissue-like complex random medium. In this paper, we discuss our comprehensive turbid polarimetry platform consisting of the experimental polarimetry system, forward Monte Carlo modelling and inverse polar decomposition analysis. Initial biomedical applications of this novel general method for polarimetry analysis in random media are also presented.     

Measurements of diffuse backscattering Mueller matrices of turbid media

We report on the development of a method that records spatially dependent intensity patterns of polarized light that arise from illuminating a turbid media with a polarized laser beam and being diffusely backscattered. Our technique employs polarized light from a He-Ne laser () which is focused onto the surface of the scattering medium. A surface area of approximately 2×2 cm centered on the light input point is imaged through polarization analysis optics onto a CCD camera. The Mueller matrix is reconstructed by 49 intensity measurements with various orientations of polarizer and analyzer. The measured Mueller matrix of polystyrene spheres was compared with the theory result of incoherent scattering of light by spheres. The experimental and theory results are in excellent agreement. It appears that the azimuthal patterns of the Mueller matrix are determined by the symmetry of the turbid media.     

Theoretical analysis of backscattered polarization patterns of turbid media containing glucose

Single scattering model and Stokes-Mueller formalism are introduced to investigate the influence of glucose on backscattered polarization patterns in turbid media. Glucose molecules rotate the polarization plane and induce changes in backscattered Mueller matrix patterns. Some Mueller matrix elements present higher optical rotation as the concentration of glucose augments. Using image subtraction and integration,linear relationship between low glucose concentration in the physiological range and optical rotation degree can be derived.     

Application of the polar decomposition to light scattering particle systems

We have applied the polar decomposition theorem (PDT) to the scattering matrix of simple systems such as an isolated sphere, either metallic or dielectric. The discrete dipole approximation (DDA) method has been used as an intermediate tool to calculate these matrices. We show how the parameters introduced by the PDT for describing the scattering by a single sphere can be easily interpreted. It is also shown that the PDT provides us with an alternative frame to conventional Mueller matrix analysis, by characterizing the systems by means of a group of independent parameters representing magnitudes of simple (virtual) elements. This step is necessary in order to translate other complex system matrices into manageable information.     

Magnetic field effects on coherent backscattering of light

We have studied the influence of magneto-optical Faraday rotation on coherent backscattering of light experimentally, theoretically and by computer simulations of Monte-Carlo type. The consistency of these three approaches reveals new aspects of the propagation of vector waves in turbid media with and without Faraday rotation. Experimentally, we have demonstrated that the Faraday rotation may almost completely destroy the reciprocity of light paths. However, as shown by the simulations, the cone of coherent backscattering may not only be destroyed but also shifted off the exact backscattering direction, parallel to the magnetic field, as long as the amount of circular polarization is not completely destroyed by the multiple scattering. The relationship between coherent backscattering, depolarization and Faraday rotation are explained by a simple path model of vector waves. This leads to a new characteristic correlation length required to properly describe the influence of Faraday rotation on multiple light scattering. Received 28 January 2000     

Iterative approach of high-order scattering solution for vector radiative transfer of inhomogeneous random media

By stratifying a random scatter media into multiple thin layers in the vertical z direction, the first-order scattering solution of each thin layer is employed to derive high-order scattering solution of whole random media. Using the Fourier transform and Mueller matrices in discrete ordinates, an iterative approach to solve high-order scattering solution of vector radiative transfer (VRT) equation is newly developed. Numerical results are well compared with the Mueller matrix solutions of the first order for a single layer medium, second order for a half-space, and the results of the discrete ordinate and eigen analysis method. It demonstrates our approach as feasible, effective and especially applicable to high-order solution of VRT for both bistatic scattering and thermal emission of inhomogeneous non-spherical scatter media.     

混浊介质二维后向漫散射穆勒矩阵的测量

研制了一套实验系统用于测量混浊介质的后向漫散射穆勒矩阵.在系统中,激光被调制成所需偏振态,聚焦于混浊介质表面.利用CCD相机配合相应的偏振元件,测量入射点周围区域的后向散射光中特定偏振态的能量空间分布,并由测量值计算出介质的后向漫散射穆勒矩阵.描述了系统的组成、测量原理以及方法,分别测量了消偏振分光棱镜的反射穆勒矩阵和浓度1.5％的脂肪乳溶液的后向漫散射穆勒矩阵.实验证明:该系统能够完全消除样品表面镜面反射光的干扰,精确的测量混浊介质的后向漫散射穆勒矩阵;并可望用于生物组织的研究中.     

Estimation of the scattering coefficients of turbid media using angle-resolved optical frequency-domain imaging

Noninvasive measurements of the scattering coefficients of optically turbid media using angle-resolved optical frequency-domain imaging (OFDI) are demonstrated. It is shown that, by incoherently averaging OFDI reflectance signals acquired at different backscattering angles, speckle noise is reduced, allowing scattering coefficients to be extracted from a single A-line with much higher accuracy than with measurements from conventional OFDI and optical coherence tomography systems. Modeling speckle as a random phasor sum, the relationship between the measurement accuracy and the number of compounded angles is derived. The sensitivity analysis is validated with measurements from a tissue phantom.     

Revisit to the symmetry relations in diffusely backscattered polarization patterns of turbid media

As there exists an inconsistency in claiming the symmetrical relations in the 16 Mueller matrix elements used to describe a turbid medium, the author restudies the symmetrical relationships between diffusely backscattered polarization patterns in isotropic turbid media and simulates all two-dimensional elements of diffusely backscattered Mueller matrix in both cases of Rayleigh and Mie scatterings using the doublescattering approximation and the Monte Carlo algorithm, respectively. The previous experimental observatious are compared with the numerically determined matrix elements, showing a good agreement in both double-scattering model and Monte Carlo simulation. The symmetrical relations between the Mueller matrix elements are clarified.     

New technique to simulate angular patterns of polarized signals from multiply scattering media in MFOV and CCD lidars

We present a semi-analytical approach to determine angular patterns of a polarized intensity that form lidar signals from multiply scattering media registered by a monostatic multiple-field-of-view receiver or CCD camera for any polarization state of incident light and any analyzer state at a receiver. Multiply scattering media with a highly forward elongated phase functions are considered. The model of lidar signal formation includes the single near-backscattering and small-angle multiple scattering of light. The developed approach allows computation and analysis of polarized images, including the Mueller matrix images, formed by the laser light backscattered by multiply scattering media, for example, from atmospheric clouds, ocean waters, tissue, etc.     

Polarization characteristics of backscattering of turbid media based on Mueller matrix

Mueller matrix is one approach to characterizing optical polarization of the turbid media. We have simulated the two-dimensional images of Mueller matrix based on single-scattering approximation model and implemented experiments to verify the simulations. By comparing the experimental results to the theoretical simulations, we have obtained some conclusions. When the particle size is smaller than the wavelength, the linearly polarized light propagating through the turbid media of Rayleigh scatterers has better polarization-maintaining ability. Whereas when the particle size is larger than the wavelength, the circularly polarized light propagating through the turbid media of Mie scatterers has better polarization-maintaining ability. Moreover, the radial dependence of the element patterns becomes weak as the transport mean free path decreases. This study can help us understand to the fundamental principle of optical polarization.     

手性介质中后向散射米勒矩阵特性及其在血糖无创检测中的应用初探

利用斯托克斯-米勒矩阵表述分析偏振光在手性介质中的传输规律,利用瑞利近似下的单次散射模型和基于米氏散射理论的多次散射的蒙特卡罗方法,探讨葡萄糖对后向散射偏振光传输特性的影响,计算了不同葡萄糖浓度下的后向散射米勒矩阵. 结果表明在含有葡萄糖的手性介质中后向散射米勒矩阵的若干矩阵元表现出微弱的旋光效应,旋光度随葡萄糖浓度增大而增大. 为了检测低浓度下微小的旋光变化,定义了函数对含-不含葡萄糖的米勒矩阵元图像进行处理,提取由葡萄糖引入的图像差异. 结合无创血糖检测应用要求,分析了葡萄糖生理浓度下的矩阵元图像,获 关键词： 医用光学与生物技术 偏振光 斯托克斯-米勒表述 蒙特卡罗模拟     

Double-scatter Mueller matrices for vector electromagnetic scattering from single grooves and lines on a silicon surface

A recently developed numerical method is used to calculate the single- and double-scattered Mueller matrices for scattering of vector-electromagnetic waves from rough surfaces. Calculations are performed for the case of a single groove and a single rib on a silicon surface and results are compared to the published experimental results. The calculated results show good agreement with the experimental results for the groove case and poorer agreement with the rib case, probably due to errors in the experimental construction of the rib. It is found that, for the cases studied here, variations of the sign with scatter angle in individual Mueller matrix elements are associated with the presence of double- (or multiple-)scattered light, as has been found previously for scalar diffraction calculations in Gaussian randomly rough surfaces.     

Weak localization effects in the second-harmonic light scattered by random systems of particles

We report rigorous numerical simulations that show the presence of coherent backscattering effects in the second-harmonic generation and scattering of light by random systems of two-dimensional particles. Since the medium composing the particles is assumed to be homogeneous and isotropic, the second-harmonic field is generated mainly by surface effects. For the fundamental frequency, the results present a clear enhanced backscattering peak. In contrast, the second-harmonic scattering patterns present an intensity dip in the backscattering direction.