# Evolution of Statistic Moments of 2D-Distributions of Biological Liquid Crystal Net Mueller Matrix Elements in the Process of Their Birefringent Structure Changes

Authors: A. G. Ushenko; I. Z. Misevich; V. Istratiy; I. Bachyns'ka; A. P. Peresunko; Omar Kamal Numan; T. G. Moiysuk
Journal: Advances in Optical Technologies
 (2010)

---

## Abstract

This research is aimed to investigate the reliability of Mueller-matrix differentiation of birefringence change of optically thick layers of biological liquid crystals at the early stages of the change in their physiological state. This is performed by measuring the set of skewness and kurtosis values of Mueller matrix image of the phase elementM44 in various points of the object under investigation.


---

## Body

## 1. Introduction

Among many methods of optical diagnostics of organic phase-inhomogeneous object a new technique—laser polarimetry [1]—has been formulated within recent 10 years. It enables to obtain information about optical anisotropy [2–5] of biological tissues (BT) in the form of coordinate distributions of BT Mueller matrix elements, azimuths, and ellipticities of their object field polarization.To analyze this polarimetric information the following model approach was elaborated [1, 2, 6–13]:(i)
all the variety of human BT can be represented as four main types—connective, muscle, epithelial, and nerve tissues;(ii)
morphological structure of any BT type is regarded as a 2-component amorphous-crystalline structure;(iii)
the crystalline component or extracellular matrix is an architectonic net consisting of coaxial cylindrical protein (collagen, myosin, elastin, etc.) fibrils;(iv)
optically, the protein fibrils possess the properties of uniaxial birefringent crystals;(v)
interaction of laser radiation with the BT layer is considered in the single scattering approximation, when the attenuation factor corresponds toτ≤0,1.Specifically, the above mentioned model was used for finding and substantiating the interconnections between the ensemble of statistic moments of the 1st–4th orders that characterize the orientation-phase structure (distribution of optical axes and phase shifts of protein fibrils networks directions) of birefringent BT architectonics and that of 2D distributions of the elements of the corresponding Mueller matrix [6, 10]. It was determined [8, 11] that the 3rd and the 4th statistic moments of coordinate distributions of “phase” matrix elements (z24, z34, z44) are the most sensitive to the change (dystrophic and oncological processes) of optical anisotropy of protein  crystals. These statistic moments characterize the BT extracellular matrix birefringence. On this basis the criteria of early diagnostics of muscle dystrophy, precancer states of connective tissue, collagenosis, and so forth were determined.However, such techniques do not take into account the coordinate heterogeneity of orientation-phase structure of protein crystals nets of the BT layer as well as the order of scattering in its depth. Thus it is important to investigate the distribution of statistic moments of the 1st–4th orders that characterize the 2D elements of Mueller matrix not only in the section of the probing laser beam but also within the whole BT layer of various optical thickness and physiological state.For this we shall consider the potentiality of matrix modeling of polarization properties of optically thick BT layer.


## 2. Mueller Matrix Modeling of the Properties of Optically Thick Biological Tissue

Let us represent the layer of such a biological object as the set of successively located optical thin partial layers (Figure1).Figure 1
On the analysis of modeling polarization properties of anisotropic component of optically thick biological tissue.Polarization properties of birefringent nets of every other BT layer are described by the Mueller matrix{Z}j representing a superposition of matrix operators of optically separate coaxial protein fibrils {z}j [1]:(1){z}j=∥10000z22z23z240z32z33z340z42z43z44∥,
where(2)zik(ρ,δ)={z22=cos22ρ+sin22ρcosδ,z23,32=cos2ρsin2ρ(1-cosδ),z33=sin22ρ+cos22ρcosδ,z34,43=±cos2ρsinδ,z24,42=±sin2ρsinδ,z44=cosδ.
Here ρ is the direction of optical axis determined by the direction of packing of the birefringent fibril, δ=(2π/λ)Δnd is the phase shift introduced between the orthogonal components of the amplitude of laser wave with the length λ passing through the fibril with linear size of its geometrical section d and birefringence index Δn.Mueller matrix(Zik)j elements of the net of protein fibrils of partial BT layer (j) are determined by the following algorithms: (i)
for finite number (N) of fibrils
(3)Zik=∑u=1N∑u=1Nzik(ρn,δn),(ii)
for “infinite” (N→∞) number of fibrils
(4)Zik=∫0π∫02πQρWδzik(ρ,δ)dρdδ.HereQρ and Wδ are the distribution function of orientation (ρ) and phase (δ) parameters of biological crystals network.Mueller matrix of optically thick BT is determined by multiplication of partial matrix operators:(5){Z}={Z}n{Z}n-1⋯{Z}2{Z}1.To make it simpler (without decreasing the analysis depth) further we shall consider only a 2-layered BT:(6){Z}={Z}2{Z}1≡{Y}{X}.In the expanded form the matrix operator elements (6) are written as follows:(7)Zik(ρx,ρy,δx,δy)={Z22=y22x22+y23x32+y24x42,Z23=y22x23+y23x33+y24x43,Z32=y32x22+y33x32+y34x42,Z33=y32x23+y33x33+y34x43,Z34=y32x24+y33x34+y34x44,Z43=y42x23+y43x33+y44x43,Z24=y22x24+y23x34+y24x44,Z42=y42x22+y43x32+y44x42,Z44=y42x24+y43x34+y44x44.The analysis of relations (7) shows that optical properties of anisotropic component of 2-layered BT (Zik) are described by complex superposition of orientation and phase parameters of its partial layers' biological crystals (xik(ρx,δx),yik(ρy,δy)). The concrete form of such dependences for “phase” z42,43,44 matrix elements [10] is illustrated by the following relations:(8)z24=sinδy[cos2ρxsin2(ρx-ρy)+cosδxsin2ρxcos2(ρx+ρy)]+sinδxcosδysin2ρx;z34=sinδx[sin2ρysin2(ρx-ρy)+cosδycos2ρycos2(ρx+ρy)]+cosδxsinδycos2ρy;z44=cos(2ρx-2ρy)sinδxsinδy+cosδxcosδy.The analysis of relations (8) shows that a strict solution of the inverse problem {ρx=q(zik,ρy,δy);δx=g(zik,ρy,δy)}, revealing the changes in the structure of biological crystals network of one of the layers on the basis of data about matrix elements zik and yik, is both mathematically incorrect and physically ambiguous.Thus, it is important to elaborate approximated statistical methods of experimental solution of such diagnostic task.Generally speaking, the increase of optical thickness of biological layer results in multiple scattering regime in its volume. Such process appeared in formation of integral depolarization of laser field and results in averaging within angular aperture of photodetector of azimuth and ellipticity random values of fully polarized isolated speckles [8]. In this case the direct interconnection between the scattered field polarization parameters and the orientation and phase structure of biological crystals will be lost. From the other hand, for such integrally depolarized fields the Mueller matrix approach remains adequate [1]. Therefore, the search of interconnections “object-laser field” on the level of indirect, statistic parameters, which characterizes their structure is topical.


## 3. Scheme of Experimental Changes of Coordinate Distributions of Biological Tissue Mueller Matrix Elements

The histological sections of different geometrical (optical) thickness were used as the objects of investigation. The technique of obtaining of such objects is convenient: biological tissue is freezing to nitrogen temperature with the following obtaining, by means of medical microtome, the histological sections (from 10μm to 100 μm). Such a technique allows preserving the optical properties of biological tissues over a period of 24 hours.Conventional optical scheme of polarimeter for measuring 2D distributions of the BT Mueller matrix elements is presented in Figure2 [1].Figure 2
Optical scheme of polarimeter, where 1: He-Ne laser; 2: collimator; 3: stationary quarter-wave plate; 5, 8: mechanically movable quarter-wave plates; 4, 9: polarizer and analyzer correspondingly; 6: object of investigation; 7: microobjective; 10: CCD camera; 11: PC.Histological sections of BT were illuminated by a parallel beam of He-Ne laser (λ = 0.6328 μm, W = 5.0 μW) with the radius r=1mm. Polarization illuminator consists of quarter-wave plates 3, 5 and polarizer 4, providing the formation of laser beam with random azimuth 0∘≤α0≤180∘ or ellipticity 0∘≤β0≤90∘ of polarization. Polarization images of BT by means of microobjective 7 were projected into the plane of sensitized plate (m×n=800×600 pixels) if with CCD-camera 10. The analysis of BT images was carried out by means of polarizer 9 and quarter-wave plate 8. As a result, the Stokes vector parameters for every pixel of the BT image {Si=1,2,3,4} were determined and the set of elements of Mueller matrix was calculated according to algorithm [10]:(9)Zi1=0,5[Si(1)+Si(1)],Zi2=0,5[Si(1)-Si(1)],Zi3=Sl(3)-Zi1,Zi4=Sl(4)-Zi1,i=1,2,3,4.
Indices 1–4 correspond to the following polarization states of the beam probing the BT layer: 1–0°; 2–90°; 3–+45°; 4–⊗ (right-hand circulation).The technique of measuring the ensemble (j—the number of probing areas) of 2D matrix elements Zik(j)(m×n) of BT consisted in the following sequence of actions: (i)
the plane of BT histological sections (≈10 mm × 20 mm) BT was scanned line by line by the laser beam with the radius r=1mm with the step of linear displacement Δr=2mm;(ii)
within each of 50 areas (j = 1,2,…,50) of illumination of the BT layer plane according to the algorithm (9) the local array (m×n=800×600) of values of Mueller matrix elements Zik(j)(m×n) was determined:
(10)Zik(j)(m×n)=(Zik11⋯Zik1m⋯⋯⋯Zikn1⋯Ziknm)j=1,2,...,50;(iii)
for every local arrayZik(j)(m×n) statistic moments of the 3rd-4th M(s=3,4) orders were calculated;(iv)
the histogramsN(Mj=1,2,...,50(s=3;4)) of the values of statistical moments of higher orders within the whole plane of the investigated BT histological sections were determined.


## 4. Investigation of Skewness and Kurtosis Distributions of 2D Phase Element of Mueller Matrix of the Samples of Skeletal Muscle Tissue

The following histological sections of different types of rat’s tissues were investigated—connection tissue, muscular tissue, epithelial tissue, and nervous tissue. All these tissues have common optical peculiarity—the presence of birefringent network of protein fibrils [1]. Therefore, in this work, without decrease of analysis completeness, we have concentrated on investigation of optical anisotropy of rat's skeletal muscle tissue birefringent fibrils.DistributionsN(M(s=3,4)) of the values of skewness M(s=3) and kurtosis M(s=4) of 2D phase (Z44) element of Mueller matrix of rat's skeletal muscle tissue were investigated. Histological sections of various optical thickness (τ1=0,08 and τ2=1,47) and physiological state were the object of investigation. Such choice of optical thickness of histological sections (d1=15μm;d2=50μm) allows us to investigate the influence of a single (τ1=0,08) and multiple (τ2=1,47) scattering on coordinate structure of Mueller matrix elements.The choice of matrix elementZ44 as an analytical parameter is explained by the fact that it is the most sensitive to the changes of birefringence of protein fibrils nets, connected with their pathological changes [1, 11]. The choice of a skeletal muscle tissue of a rat as the object of investigation is connected with the possibility of direct experimental formation of septic inflammation and monitoring the control of their optical manifestations under the condition of scattering of laser beam of various multiplication factor.The series of Figures3 and 4 presents the distributions of the element Z44(m,n) of optically thin histological sections of healthy (Figure 3) and inflamed (Figure 4) muscle tissue.Coordinate and 3D distribution of the values of phase matrix elementZ44(m,n) of histological section of optically thin (τ1=0,08) layer of the healthy rat's muscle tissue.
(a)(b)Coordinate and 3D distribution of the values of phase matrix elementZ44(m,n) of histological section of optically thin (τ1=0,08) layer of the inflamed rat's muscle tissue.
(a)(b)The following can be seen from the data obtained:(i)
2D distributions of matrix elementsZ44(m,n) of both types of tissues are coordinately heterogeneous;(ii)
for the sample of healthy muscle tissue (Figure3) the lesser range of values change (from −0.7 to 0.7) of the matrix element Z44(m,n) is typical in comparison with the changes (from −0.85 to 0.85) of the analogous element of Mueller matrix of the inflamed tissue (Figure 4).The structure of 2D distributions of elementsZ44(m,n) can be connected with the following peculiarities of birefringent architectonics of the skeletal tissue samples. Firstly, while a laser wave (with the wave-length λ) propagates through the network of optically anisotropic (birefringence index Δn=1,5×10-3 [7]) myosin fibrils, a wide range of values of phase shifts δk, proportional to their geometrical sizes dk, is formed(11)δk~(2πΔnλ)×{d1⋮dk.Thus, phase elementsZ44(m,n) distributions are coordinately heterogeneous and dependent on the peculiarities of morphological structure of extracellular matrix of muscle tissue samples.Secondly, the inflammation process is accompanied by the increase of birefringenceΔn due to formation of myosin fibrils edema [7]. Optically, it causes the increase of the range of phase shifts δk changes and the growth of fluctuations of the element Z44(m,n) connected with it (Figure 4).The histograms of distributionN(M(s=3,4)) of statistic moments of higher-orders M(s=3,4) of coordinate distributions Z(j=1÷50)44(m,n) of investigated samples of optically thin histological sections are presented in Figure 5.Histograms of the set of skewnessM(s=3) values (a, b) and kurtosis M(s=4) (c, d) of coordinate distribution of the element Z(j=1-50)44(m,n) of histological section of optically thin (τ1=0,08) layer of the healthy (a, c) and inflamed (b, d) rat's skeletal tissue.
(a)(b)(c)(d)The analysis of experimentally measured histogramsN(M(s=3,4)) proved the following:(i)
the change ranges of the values of skewnessM(s=3) and kurtosis M(s=4) of coordinate distributions Z(j=1-50)44(m,n), measured for the healthy (a,c) and inflamed (b,d) muscle tissues, do not actually coincide;(ii)
extreme values of skewness in the distributionN(M(s=3)) for phase element Z44 of the healthy tissue are by 4–6 times higher than analogous values of the given statistic moment M(s=3) for the inflamed muscle tissue;(iii)
extreme values of kurtosis in the distributionN(M(s=4)) for the inflamed tissue are by 2-3 times higher than the values of the statistic moment M(s=4) for the healthy skeletal tissue.Thus, it can be stated that for single differentiation of physiological state of optically thin layers of muscle tissue it is enough to measure the 2D phase matrix elementZ44(m,n) in one domain (πr2) of irradiation by the laser beam and to calculate the skewness and kurtosis of its values distribution.The series of Figures6 and 7 presents the results of experimental measurements of phase elements Z44(m,n) of Mueller matrix of optically thick (τ2=1,47) of the rat's muscle tissue layers.Coordinate and 3D distribution of the values of phase matrix elementZ44(m,n) of histological section of optically thick (τ1=1,47) layer of the healthy rat's muscle tissue.
(a)(b)Coordinate and 3D distribution of the values of phase matrix elementZ44(m,n) of histological section of optically thick (τ1=1,47) layer of the inflamed rat's muscle tissue.
(a)(b)It can be seen from the data obtained the following:(i)
for the samples of the healthy and inflamed muscle tissue the same range of values change (from −0,55 to 0,55) is typical in the distribution of matrix elementsZ44(m,n);(ii)
comparative visual analysis of coordinate distributions of matrix elementsZ44(m,n) of histological sections of skeletal tissue of both types did not show any sufficient difference between them.Similarity of the distribution structureZ44(m,n) of both types of rat's skeletal tissue samples can be explained by multiple light-scattering. As a result of every local (ith) act of interaction between laser radiation and separate fibril random value of phase shift δ(i) is formed, which is multiplied (δ*=∑i=1Rδi) in the process of propagation in the BT depth reaching equiprobable random values from 0 to 2kπ, k=1,2,3….Figure8 shows histograms N(M(s=3,4)) of the set of values of statistic moments of skewness M(s=3) and kurtosis M(s=4) of distributions of the matrix element Z(j=1-50)44(m,n) of optically thick samples of muscle tissue.Histograms of the set of skewnessM(s=3) values (a, b) and kurtosis M(s=4) (c, d) of coordinate distribution of the element Z(j=1-50)44(m,n) of histological section of optically thick (τ2=1,17) layer of the healthy (a, c) and inflamed (b, d) rat's skeletal tissue.
(a)(b)(c)(d)The analysis of experimental data showed the following:(i)
the skewnessM(s=3) values distribution of 2D phase Z(j=1-50)44(m,n), and changes within the whole plane of samples of the layer of the healthy (a) and inflamed (b) muscle tissue practically coincide and cannot be used as an objective criterion for their optical properties differentiation;(ii)
histogramsN(M(s=4)) of the values of the statistic moment M(s=4) of elements Z(j=1-50)44(m,n) distribution possess individual structure that depends on physiological state of muscle tissue;(iii)
extreme values of kurtosisM(s=4) in the histogram N(M(s=4)) of the healthy tissue (c) are by 4-5 times less than analogous values of the given statistic moment of the inflamed rat's skeletal tissue (d).


## 5. Conclusion

Thus for realization of Mueller-matrix differentiation of optically thick BT layers at early stages of changing their physiological state the measurement of distributions of the kurtosisN(M(s=4)) of the ensemble of phase elements Z(j=1÷50)44(m,n) is objective.The Mueller-matrix approach in diagnostics and differentiation of the biological crystals networks appears to be effective for optically thick histological sections layers of human tissues (skin derma, tissues of women reproductive sphere, parenchymatous tissues) also.



---

*Source: 423145-2010-06-07.xml*