Elimination of cross-talk in diffuse optical tomographic reconstructions through a weighted update procedure: results of simulation study - art. no. 61640S

Reconstructions in optical tomography involve obtaining the images of absorption and reduced scattering coefficients. The integrated intensity data has greater sensitivity to absorption coefficient variations than scattering coefficient. However, the sensitivity of intensity data to scattering coefficient is not zero. We considered an object with two inhomogeneities (one in absorption and the other in scattering coefficient). The standard iterative reconstruction techniques produced results, which were plagued by cross talk, i.e., the absorption coefficient reconstruction has a false positive corresponding to the location of scattering inhomogeneity, and vice-versa. We present a method to remove cross talk in the reconstruction, by generating a weight matrix and weighting the update vector during the iteration. The weight matrix is created by the following method: we first perform a simple backprojection of the difference between the experimental and corresponding homogeneous intensity data. The built up image has greater weightage towards absorption inhomogeneity than the scattering inhomogeneity and its appropriate inverse is weighted towards the scattering inhomogeneity. These two weight matrices are used as multiplication factors in the update vectors, normalized backprojected image of difference intensity for absorption inhomogeneity and the inverse of the above for the scattering inhomogeneity, during the image reconstruction procedure. We demonstrate through numerical simulations, that cross-talk is fully eliminated through this modified reconstruction procedure.

Synthesis and characterisation of co-evaporated tin sulphide thin films

Tin sulphide films were grown at different substrate temperatures by a thermal co-evaporation technique. The crystallinity of the films was evaluated from X-ray diffraction studies. Single-phase SnS films showed a strong (040) orientation with an orthorhombic crystal structure and a grain size of 0.12 mu m. The films showed an electrical resistivity of 6.1 Omega cm with an activation energy of 0.26 eV. These films exhibited an optical band gap of 1.37 eV and had a high optical absorption coefficient (> 10(4) cm(-1)) above the band-gap energy. The results obtained were analysed to evaluate the potentiality of the co-evaporated SnS films as an absorber layer in solar photovoltaic devices.

Photo induced optical changes in Sb/As2S3 multilayered film and (As2S3)(0.93)Sb-0.07 film of equal thickness

The increase in optical band gap (photo bleaching) due to light illumination was studied at room temperature as well as at low (4.2 K) temperature for Sb/As2S3 multilayered film of 640 nm thickness by Fourier Transform Infrared Technique. The interdiffusion of Sb into As2S3 matrix results the formation of Sb-As2S3 ternary solid solutions which is explained by the change in optical band gap (E-g), absorption coefficient (alpha), Tauc parameter (B-1/2), Urbach edge (E-e). At the same time, photo darkening phenomena was observed in (As2S3)(0.93)Sb-0.07 film of same thickness both at low and room temperatures. From our X-ray Photoelectron Spectroscopy measurements,we are able to show that some of the As-As, S-S and Sb-Sb bonds are converted into As-S and S-Sb bonds in case of multilayers. We found that the energetically favoured heteropolar bond formation take place by a phonon-assisted mechanism using the lone pair pi electrons of S-2(0). But in case of (As2S3)(0.93)Sb-0.02 film, the homopolar bonds are playing a major role. (C) 2010 Elsevier B.V. All rights reserved.

Effect of sodium diffusion on the structural and electrical properties of Cu2ZnSnS4 thin films

Cu2ZnSnS4 (CZTS) is a kesterite semiconductor consisting of abundantly available elements. It has a band gap of 1.5 eV and a large absorption coefficient. Hence, thin films made of this material can be used as absorber layers of a solar cell. CZTS films were deposited on soda lime and Na free borosilicate glass substrates through Ultrasonic Spray Pyrolysis. The diffusion of sodium from soda lime glass was found to have a profound effect on characteristics like grain size, crystal texture and conductivity of CZTS thin films. Copper ion concentration also varied during the deposition and it was observed that the carrier concentration was enhanced when there was a deficiency of copper in the films. The effect of sodium diffusion and copper deficiency in enhancing the structural and electrical properties of CZTS films are presented in this paper. (C) 2010 Elsevier B.V. All rights reserved.

Ultrasonic spray pyrolysis of CZTS solar cell absorber layers and characterization studies

CZTS (Copper Zinc Tin Sulphide) is a wide band gap quartnery chalcopyrite which has a band gap of about 1.45 eV and an absorption coefficient of 10(4) cm(-1); thus making it an ideal material to be used as an absorber layer in solar cells. Ultrasonic Spray Pyrolysis is a deposition technique, where the solution is atomized ultrasonically, thereby giving a fine mist having a narrow size distribution which can be used for uniform coatings on substrates. An Ultrasonic Spray Pyrolysis equipment was developed and CZTS absorber layers were successfully grown with this technique on soda lime glass substrates using aqueous solutions. Substrate temperatures ranging from 523 K to 723 K were used to deposit the CZTS layers and these films were characterized using SEM, EDAX and XRD. It was observed that the film crystallized in the kesterite structure and the best crystallites were obtained at 613 K. It was observed that the grain size progressively increased with temperature. The optical band gap of the material was obtained as 1.54 eV.

Estimation of elasticity map of soft biological tissue mimicking phantom using laser speckle contrast analysis

Scattering of coherent light from scattering particles causes phase shift to the scattered light. The interference of unscattered and scattered light causes the formation of speckles. When the scattering particles, under the influence of an ultrasound (US) pressure wave, vibrate, the phase shift fluctuates, thereby causing fluctuation in speckle intensity. We use the laser speckle contrast analysis (LSCA) to reconstruct a map of the elastic property (Young's modulus) of soft tissue-mimicking phantom. The displacement of the scatters is inversely related to the Young's modulus of the medium. The elastic properties of soft biological tissues vary, many fold with malignancy. The experimental results show that laser speckle contrast (LSC) is very sensitive to the pathological changes in a soft tissue medium. The experiments are carried out on a phantom with two cylindrical inclusions of sizes 6 mm in diameter, separated by 8 mm between them. Three samples are made. One inclusion has Young's modulus E of 40 kPa. The second inclusion has either a Young's modulus E of 20 kPa, or scattering coefficient of mu'(s), = 3.00 mm(-1) or absorption coefficient of mu(a) = 0.03 mm(-1). The optical absorption (mu(a)), reduced scattering (mu'(s)) coefficient, and the Young's modulus of the background are mu(a) = 0.01 mm(-1), mu'(s) = 1.00 mm(-1) and 12kPa, respectively. The experiments are carried out on all three phantoms. On a phantom with two inclusions of Young's modulus of 20 and 40 kPa, the measured relative speckle image contrasts are 36.55% and 63.72%, respectively. Experiments are repeated on phantoms with inclusions of mu(a) = 0.03 mm-1, E = 40 kPa and mu'(s) = 3.00 mm(-1). The results show that it is possible to detect inclusions with contrasts in optical absorption, optical scattering, and Young's modulus. Studies of the variation of laser speckle contrast with ultrasound driving force for various values of mu(a), mu'(s), and Young's modulus of the tissue mimicking medium are also carried out. (C) 2011 American Institute of Physics. doi:10.1063/1.3592352]

Pseudo-time particle filtering for diffuse optical tomography

We recast the reconstruction problem of diffuse optical tomography (DOT) in a pseudo-dynamical framework and develop a method to recover the optical parameters using particle filters, i.e., stochastic filters based on Monte Carlo simulations. In particular, we have implemented two such filters, viz., the bootstrap (BS) filter and the Gaussian-sum (GS) filter and employed them to recover optical absorption coefficient distribution from both numerically simulated and experimentally generated photon fluence data. Using either indicator functions or compactly supported continuous kernels to represent the unknown property distribution within the inhomogeneous inclusions, we have drastically reduced the number of parameters to be recovered and thus brought the overall computation time to within reasonable limits. Even though the GS filter outperformed the BS filter in terms of accuracy of reconstruction, both gave fairly accurate recovery of the height, radius, and location of the inclusions. Since the present filtering algorithms do not use derivatives, we could demonstrate accurate contrast recovery even in the middle of the object where the usual deterministic algorithms perform poorly owing to the poor sensitivity of measurement of the parameters. Consistent with the fact that the DOT recovery, being ill posed, admits multiple solutions, both the filters gave solutions that were verified to be admissible by the closeness of the data computed through them to the data used in the filtering step (either numerically simulated or experimentally generated). (C) 2011 Optical Society of America

Activity of TiO2 in the System ZrO2-TiO2 and Gibbs Energy of Formation of ZrTiO4

The activity of Ti02 in single and two··phase regions of ihe system ZrOrTi02 has heen measured lIsing solid state cells based on yttria··doped tho ria (YDT) as the solid electrolyte at 1373 K. The cells used can be represented as: Pt, Tio.07PtO.Y3 + Zrj.,Tix0 2 / YDT / Ti02 + Tio.07Pto.93, Pt Pt, Tio.07Pto.93 + ZrJ.xTix02 + ZrTi04 / YDT / Ti02+ Tio.07PtO.93, Pt In each cell the composition of Pt-Ti alloy was identical at hoth electrodes. The emf of the cell is therefore directly related to the activity of Ti02 in oxide phase or oxide phase mixture: aTiO~ :;: cxp (-4FE/RT). The activity coefficient of Ti02 in th~ zirconia-rich solid solution with monoclinic structure (CUl2 2" XTi02 2" 0) can be expressed as:In the zirconia-rich solid solution with tetragonal structure (0.085 2" X ri02 2" 0.03), the activity coefficient is given by:In YTi02 (± 0.012) = 2.354 (1-XTiO? )2 +0.064 The standard Gibbs energy of formation of ZrTi04 is -5650 (± 200) J/mol at 1373 K .

Analysis of light‐induced processes in bacteriorhodopsin and its application for spatial light modulation

A simplified energy‐level scheme is proposed for the photochemical cycle of the bacteriorhodopsin molecule. Rate equations are solved for the detailed light‐induced processes based on this model and the intensity‐induced population densities in various states of the molecule at steady state are computed which are used to obtain an analytic expression for the absorption coefficient of the modulation beam. Modulation of the probe laser‐beam transmission by the modulation‐laser‐beam intensity‐induced population changes is analyzed. It is predicted that for a probe beam at 412 nm up to 82% modulation can be achieved using a laser beam intensity of 3.2 W/cm2 at 570 nm. For temperatures ∼77 K, the transmission at 610 nm can be switched from zero to 81% for modulating laser intensity of 11 W/cm2. Construction of a spatial light modulator based on bacteriorhodopsin molecules is proposed and some of its features are discussed.

Implementation of a phase array diffuse optical tomographic imager

Diffuse optical tomography (DOT) using near-infrared (NIR) light is a promising tool for noninvasive imaging of deep tissue. This technique is capable of quantitative reconstructions of absorption coefficient inhomogeneities of tissue. The motivation for reconstructing the optical property variation is that it, and, in particular, the absorption coefficient variation, can be used to diagnose different metabolic and disease states of tissue. In DOT, like any other medical imaging modality, the aim is to produce a reconstruction with good spatial resolution and accuracy from noisy measurements. We study the performance of a phase array system for detection of optical inhomogeneities in tissue. The light transport through a tissue is diffusive in nature and can be modeled using diffusion equation if the optical parameters of the inhomogeneity are close to the optical properties of the background. The amplitude cancellation method that uses dual out-of-phase sources (phase array) can detect and locate small objects in turbid medium. The inverse problem is solved using model based iterative image reconstruction. Diffusion equation is solved using finite element method for providing the forward model for photon transport. The solution of the forward problem is used for computing the Jacobian and the simultaneous equation is solved using conjugate gradient search. The simulation studies have been carried out and the results show that a phase array system can resolve inhomogeneities with sizes of 5 mm when the absorption coefficient of the inhomogeneity is twice that of the background tissue. To validate this result, a prototype model for performing a dual-source system has been developed. Experiments are carried out by inserting an inhomogeneity of high optical absorption coefficient in an otherwise homogeneous phantom while keeping the scattering coefficient same. The high frequency (100 MHz) modulated dual out-of-phase laser source light is propagated through the phantom. The interference of these sources creates an amplitude null and a phase shift of 180° along a plane between the two sources with a homogeneous object. A solid resin phantom with inhomogeneities simulating the tumor is used in our experiment. The amplitude and phase changes are found to be disturbed by the presence of the inhomogeneity in the object. The experimental data (amplitude and the phase measured at the detector) are used for reconstruction. The results show that the method is able to detect multiple inhomogeneities with sizes of 4 mm. The localization error for a 5 mm inhomogeneity is found to be approximately 1 mm.

Diffuse optical tomographic imager using a single light source

Near-infrared diffuse optical tomography (DOT) technique has the capability of providing good quantitative reconstruction of tissue absorption and scattering properties with additional inputs such as input and output modulation depths and correction for the photon leakage. We have calculated the two-dimensional (2D) input modulation depth from three-dimensional (3D) diffusion to model the 2D diffusion of photons. The photon leakage when light traverses from phantom to the fiber tip is estimated using a solid angle model. The experiments are carried for single (5 and 6 mm) as well as multiple inhomogeneities (6 and 8 mm) with higher absorption coefficient in a homogeneous phantom. Diffusion equation for photon transport is solved using finite element method and Jacobian is modeled for reconstructing the optical parameters. We study the development and performance of DOT system using modulated single light source and multiple detectors. The dual source methods are reported to have better reconstruction capabilities to resolve and localize single as well as multiple inhomogeneities because of its superior noise rejection capability. However, an experimental setup with dual sources is much more difficult to implement because of adjustment of two out of phase identical light probes symmetrically on either side of the detector during scanning time. Our work shows that with a relatively simpler system with a single source, the results are better in terms of resolution and localization. The experiments are carried out with 5 and 6 mm inhomogeneities separately and 6 and 8 mm inhomogeneities both together with absorption coefficient almost three times as that of the background. The results show that our experimental single source system with additional inputs such as 2D input/output modulation depth and air fiber interface correction is capable of detecting 5 and 6 mm inhomogeneities separately and can identify the size difference of multiple inhomogeneities such as 6 and 8 mm. The localization error is zero. The recovered absorption coefficient is 93% of inhomogeneity that we have embedded in experimental phantom.

A pseudo-time EnKF incorporating shape based reconstruction for diffuse optical tomography

Purpose: The authors aim at developing a pseudo-time, sub-optimal stochastic filtering approach based on a derivative free variant of the ensemble Kalman filter (EnKF) for solving the inverse problem of diffuse optical tomography (DOT) while making use of a shape based reconstruction strategy that enables representing a cross section of an inhomogeneous tumor boundary by a general closed curve. Methods: The optical parameter fields to be recovered are approximated via an expansion based on the circular harmonics (CH) (Fourier basis functions) and the EnKF is used to recover the coefficients in the expansion with both simulated and experimentally obtained photon fluence data on phantoms with inhomogeneous inclusions. The process and measurement equations in the pseudo-dynamic EnKF (PD-EnKF) presently yield a parsimonious representation of the filter variables, which consist of only the Fourier coefficients and the constant scalar parameter value within the inclusion. Using fictitious, low-intensity Wiener noise processes in suitably constructed ``measurement'' equations, the filter variables are treated as pseudo-stochastic processes so that their recovery within a stochastic filtering framework is made possible. Results: In our numerical simulations, we have considered both elliptical inclusions (two inhomogeneities) and those with more complex shapes (such as an annular ring and a dumbbell) in 2-D objects which are cross-sections of a cylinder with background absorption and (reduced) scattering coefficient chosen as mu(b)(a)=0.01mm(-1) and mu('b)(s)=1.0mm(-1), respectively. We also assume mu(a) = 0.02 mm(-1) within the inhomogeneity (for the single inhomogeneity case) and mu(a) = 0.02 and 0.03 mm(-1) (for the two inhomogeneities case). The reconstruction results by the PD-EnKF are shown to be consistently superior to those through a deterministic and explicitly regularized Gauss-Newton algorithm. We have also estimated the unknown mu(a) from experimentally gathered fluence data and verified the reconstruction by matching the experimental data with the computed one. Conclusions: The PD-EnKF, which exhibits little sensitivity against variations in the fictitiously introduced noise processes, is also proven to be accurate and robust in recovering a spatial map of the absorption coefficient from DOT data. With the help of shape based representation of the inhomogeneities and an appropriate scaling of the CH expansion coefficients representing the boundary, we have been able to recover inhomogeneities representative of the shape of malignancies in medical diagnostic imaging. (C) 2012 American Association of Physicists in Medicine. [DOI: 10.1118/1.3679855]

Structural Transformations, Composition Anomalies and a Dramatic Collapse of Linear Polymer Chains in Dilute Ethanol-Water Mixtures

Water-ethanol mixtures exhibit many interesting anomalies, such as negative excess partial molar volume of ethanol, excess sound absorption coefficient at low concentrations, and positive deviation from Raoult's law for vapor pressure, to mention a few. These anomalies have been attributed to different, often contradictory origins, but a quantitative understanding is still lacking. We show by computer simulation and theoretical analyses that these anomalies arise from the sudden emergence of a bicontinuous phase that occurs at a relatively low ethanol concentration of x(eth) approximate to 0.06-0.10 (that amounts to a volume fraction of 0.17-0.26, which is a significant range!). The bicontinuous phase is formed by aggregation of ethanol molecules, resulting in a weak phase transition whose nature is elucidated. We find that the microheterogeneous structure of the mixture gives rise to a pronounced nonmonotonic composition dependence of local compressibility and nonmonotonic dependence in the peak value of the radial distribution function of ethyl groups. A multidimensional free energy surface of pair association is shown to provide a molecular explanation of the known negative excess partial volume of ethanol in terms of parallel orientation and hence better packing of the ethyl groups in the mixture due to hydrophobic interactions. The energy distribution of the ethanol molecules indicates additional energy decay channels that explain the excess sound attenuation coefficient in aqueous alcohol mixtures. We studied the dependence of the solvation of a linear polymer chain on the composition of the water-ethanol solvent. We find that there is a sudden collapse of the polymer at x(eth) approximate to 0.05-a phenomenon which we attribute to the formation of the microheterogeneous structures in the binary mixture at low ethanol concentrations. Together with recent single molecule pulling experiments, these results provide new insight into the behavior of polymer chain and foreign solutes, such as enzymes, in aqueous binary mixtures.

Practical fully three-dimensional reconstruction algorithms for diffuse optical tomography

We have developed an efficient fully three-dimensional (3D) reconstruction algorithm for diffuse optical tomography (DOT). The 3D DOT, a severely ill-posed problem, is tackled through a pseudodynamic (PD) approach wherein an ordinary differential equation representing the evolution of the solution on pseudotime is integrated that bypasses an explicit inversion of the associated, ill-conditioned system matrix. One of the most computationally expensive parts of the iterative DOT algorithm, the reevaluation of the Jacobian in each of the iterations, is avoided by using the adjoint-Broyden update formula to provide low rank updates to the Jacobian. In addition, wherever feasible, we have also made the algorithm efficient by integrating along the quadratic path provided by the perturbation equation containing the Hessian. These algorithms are then proven by reconstruction, using simulated and experimental data and verifying the PD results with those from the popular Gauss-Newton scheme. The major findings of this work are as follows: (i) the PD reconstructions are comparatively artifact free, providing superior absorption coefficient maps in terms of quantitative accuracy and contrast recovery; (ii) the scaling of computation time with the dimension of the measurement set is much less steep with the Jacobian update formula in place than without it; and (iii) an increase in the data dimension, even though it renders the reconstruction problem less ill conditioned and thus provides relatively artifact-free reconstructions, does not necessarily provide better contrast property recovery. For the latter, one should also take care to uniformly distribute the measurement points, avoiding regions close to the source so that the relative strength of the derivatives for measurements away from the source does not become insignificant. (c) 2012 Optical Society of America

Tailoring the optical properties of amorphous heavily Er3+-doped Ge-Ga-S thin films

This study deals with the influence of Er-doping level and thermal annealing on the optical properties of amorphous Ge-Ga-S thin films. Nominal compositions of (GeS2)(75)(Ga2S3)(25) doped with high concentrations of 2.1 and 2.4 mol% Er2S3 (corresponding to 1.2 and 1.4 at% Er, respectively) have been chosen for this work. The results have been related to those obtained for the un-doped samples. The values of the refractive index, the absorption coefficient and optical band gap have been determined from the transmittance data. It has been found that the optical band gap of un-doped and 2.1 mol% Er2S3-doped films slightly increases with annealing temperature, whereas at 2.4 mol% Er2S3-doping level it is decreased. The dependences of the optical parameters on the erbium concentration and effect of annealing in the temperature range of 100-200 degrees C have been evaluated and discussed in relation to possible structural changes.