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.

Influence of substrate temperature on the properties of oxygen‐ion‐assisted deposited CeO2 films

Substrate temperature and ion bombardment during deposition have been observed to modify significantly the optical and structural properties of dielectric thin films. Single‐layer films of CeO2 have been deposited by electron beam evaporation with simultaneous oxygen‐ion bombardment using a Kaufman broad beam ion source and maintaining the substrates at elevated temperature. A systematic study has been made on the influence of (a) substrate temperature in the range ambient to 300 °C, (b) ion energy in the range 300–700 eV, and (c) ion current density 100–220 μA/cm2 on optical properties such as refractive index, extinction coefficient, inhomogeneity, packing density, and structural properties. The refractive index increased with in increase in substrate temperature: ion energy up to 600 eV and ion current density. Homogeneous, absorption free and high index (2.48) films have been obtained at 600 eV, 220 μA/cm2 and at substrate temperature of 300 °C. The packing density of the films was observed to be unity for the same deposition conditions. Substrate temperature with simultaneous ion bombardment modified the structure of the films from highly ordered to fine grain structure.

Thermoelectric properties of indium filled and germanium doped Co4Sb12 skutterudites

In-filled and Ge-doped Co4Sb12 skutterudites materials were synthesized by an induction melting process which was followed by annealing at 650 degrees C for 7 days. A structural, compositional, and morphological study was carried out by X-ray diffraction (XRD), electron probe micro analysis (EPMA), and scanning electron microscopy (SEM). The formation of a single skutterudite phase (delta-CoSb3) was confirmed by XRD and the composition of all the samples was verified by EPMA. The homogeneity and morphology of the samples was observed by potential Seebeck microprobe (PSM) and SEM, respectively. The PSM result confirmed the inhomogeneity of the samples. The temperature dependence of the Seebeck coefficient, electrical conductivity, and thermal conductivity were measured in the temperature range of 300-650 K. The samples of In0.16Co4Sb12-xGex (x = 0.05, 0.1, and 0.2) show a negative Seebeck coefficient confirming an n-type conductivity and the In0.16Co4Sb11.7Ge0.3 sample shows a positive Seebeck coefficient confirming a p-type conductivity. There was a change in the Seebeck coefficient from an n-type to a p-type at the doping concentration of x = 0.3 due to the excess Ge which increases in hole carrier concentration. Electrical conductivity decreases with an increase in Ge doping concentrations and with increases in temperature due to the bipolar effect. Thermal conductivity increases with an increase in carrier concentration and decreases when the temperature is increased. The highest ZT = 0.58 was achieved by In0.16Co4 Sb11.95Ge0.05 at 673K and In-filled and Ge-doped Co4Sb12 was not effective in improving the figure of merit. (C) 2012 American Institute of Physics. [doi: 10.1063/1.3677982]

Decoupling of diffusion from viscosity: Difference scenario for translational and rotational motions

Recent experiments have indicated a dramatically different viscosity dependence of the translational and the rotational diffusion coefficients in a supercooled liquid as the glass transition temperature is approached from above. While the translational motion seems to be decoupled from the rising viscosity (eta), the rotational motion seems to remain firmly coupled to eta. In order to understand the microscopic origin of this behavior, we have carried nut detailed theoretical calculations of both the quantities by using a self-consistent mode-coupling theory (MCT). it is found that when the size of the solute is same as that of the solvent molecules, the conventional MCT fails to predict the observed decoupling. The solvent inhomogeneity is found to play a decisive role in determining the decoupling. The difference in the viscosity dependence between rotation and translational diffusion coefficient is discussed.

Studying the resistivity imaging of chicken tissue phantoms with different current patterns in Electrical Impedance Tomography (EIT)

A current injection pattern in Electrical Impedance Tomography (EIT) has its own current distribution profile within the domain under test. Hence, different current patterns have different sensitivity, spatial resolution and distinguishability. Image reconstruction studies with practical phantoms are essential to assess the performance of EIT systems for their validation, calibration and comparison purposes. Impedance imaging of real tissue phantoms with different current injection methods is also essential for better assessment of the biomedical EIT systems. Chicken tissue paste phantoms and chicken tissue block phantoms are developed and the resistivity image reconstruction is studied with different current injection methods. A 16-electrode array is placed inside the phantom tank and the tank is filled with chicken muscle tissue paste or chicken tissue blocks as the background mediums. Chicken fat tissue, chicken bone, air hole and nylon cylinders are used as the inhomogeneity to obtained different phantom configurations. A low magnitude low frequency constant sinusoidal current is injected at the phantom boundary with opposite and neighboring current patterns and the boundary potentials are measured. Resistivity images are reconstructed from the boundary data using EIDORS and the reconstructed images are analyzed with the contrast parameters calculated from their elemental resistivity profiles. Results show that the resistivity profiles of all the phantom domains are successfully reconstructed with a proper background resistivity and high inhomogeneity resistivity for both the current injection methods. Reconstructed images show that, for all the chicken tissue phantoms, the inhomogeneities are suitably reconstructed with both the current injection protocols though the chicken tissue block phantom and opposite method are found more suitable. It is observed that the boundary potentials of the chicken tissue block phantoms are higher than the chicken tissue paste phantom. SNR of the chicken tissue block phantoms are found comparatively more and hence the chicken tissue block phantom is found more suitable for its lower noise performance. The background noise is found less in opposite method for all the phantom configurations which yields the better resistivity images with high PCR and COC and proper IRMean and IRMax neighboring method showed higher noise level for both the chicken tissue paste phantoms and chicken tissue block phantoms with all the inhomogeneities. Opposite method is found more suitable for both the chicken tissue phantoms, and also, chicken tissue block phantoms are found more suitable compared to the chicken tissue paste phantom. (C) 2012 Elsevier Ltd. All rights reserved.

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]

Microstructure and texture evolution during accumulative roll bonding of aluminium alloys AA2219/AA5086 composite laminates

Accumulative roll bonding of two aluminium alloys, AA2219 and AA5086 was carried out up to 8 passes. During the course of ARB, the deformation inhomogeneity between the two alloy layers results in interfacial instability after the 4th pass, necking of the AA5086 layers after the 6th pass and fracture along the necked regions after the 7th and 8th pass. The EBSD analysis shows deformation bands along the interfaces after 8 passes of ARB. The ARB-processed materials predominantly show characteristic deformation texture components. The weak texture after the 2nd pass results from the combination of a weakly-textured starting AA2219 layer and a strongly-textured starting AA5086 layer. A strong deformation texture forms due to the high imposed strain after a higher number of ARB passes. Subgrain formation and related shear banding induces copper/S components in the case of the small elongated grains, while planar slip leads to the formation of brass component in the large elongated grains.

Evidence of nanoscale structural phase separation in large bandwidth La0.2Sr0.8MnO3

We have investigated the structural evolution of La0.2Sr0.8MnO3 using temperature dependent high resolution synchrotron x-ray diffraction technique. In a wide temperature range, La0.2Sr0.8MnO3 reveals nanoscale structural inhomogeneity consisting of cubic and tetragonal phases. The present results suggest that domains of nanometer size of the tetragonal (low temperature) phase start nucleating in the cubic (high temperature) phase even above the Neel temperature (T-N). The tetragonal phase fraction increases substantially below T-N. Detailed analysis suggests that the twinned phase is tetragonal, orbital ordered, and insulating. At temperatures below 170 K, a small amount of the cubic phase is retained. The present results reveal the significance of the connectivity between the nanoscale structural phase separation with the physical properties.

Common Ground Method of Current Injection in Electrical Impedance Tomography

Surface electrodes are essentially required to be switched for boundary data collection in electrical impedance tomography (Ell). Parallel digital data bits are required to operate the multiplexers used, generally, for electrode switching in ELT. More the electrodes in an EIT system more the digital data bits are needed. For a sixteen electrode system. 16 parallel digital data bits are required to operate the multiplexers in opposite or neighbouring current injection method. In this paper a common ground current injection is proposed for EIT and the resistivity imaging is studied. Common ground method needs only two analog multiplexers each of which need only 4 digital data bits and hence only 8 digital bits are required to switch the 16 surface electrodes. Results show that the USB based data acquisition system sequentially generate digital data required for multiplexers operating in common ground current injection method. The profile of the boundary data collected from practical phantom show that the multiplexers are operating in the required sequence in common ground current injection protocol. The voltage peaks obtained for all the inhomogeneity configurations are found at the accurate positions in the boundary data matrix which proved the sequential operation of multiplexers. Resistivity images reconstructed from the boundary data collected from the practical phantom with different configurations also show that the entire digital data generation module is functioning properly. Reconstructed images and their image parameters proved that the boundary data are successfully acquired by the DAQ system which in turn indicates a sequential and proper operation of multiplexers.

Criteria for thermally induced atomization and catastrophic breakup of acoustically levitated droplet

We report thermally induced instability leading to catastrophic breakup in acoustically levitated vaporizing fuel droplets. Change in surface tension and viscosity with increase in droplet temperature causes wide fluctuations in droplet aspect ratio. If the viscous damping of aspect ratio oscillation is not strong enough, the droplet goes through unbounded stretching. If the droplet exceeds a critical Weber number locally, a bag type and capillary wave induced atomization can occur, which leads to catastrophic breakup. A stability criterion has been established based on the inhomogeneity of Bernoulli (acoustic) pressure and surface tension of the droplet in terms of a local Weber number and Ohnesorge number. This instability is thermally induced in a droplet which does not experience instabilities without heating. (C) 2012 Elsevier Ltd. All rights reserved.

Investigation of the origin of glassiness in La0.5Sr0.5CoO3

A comprehensive magnetic study has been carried out on the two sets of La0.5Sr0.5CoO3 samples with a view to understand the origin of low temperature glassiness in the ferromagnetic state. The samples prepared by the conventional solid-state synthesis method show a low temperature shoulder in both dc magnetization as well as in the ac susceptibility measurements, which exhibit characteristics of glassiness such as the frequency dependence and memory effect. These observations suggest the existence of a distinct low temperature cluster-glass like phase within dominant ferromagnetic phase. But, once the same sample is properly homogenized by repeated grinding and annealing process, the low temperature glassy phase disappears, and it shows a pure ferromagnetic behavior. Our comparative study clearly reveals that the reentrant spin-glass like nature is not intrinsic to La0.5Sr0.5CoO3 system, in fact this is an outcome of the compositional inhomogeneity.

Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Correlated Conductance Fluctuations Close to the Berezinskii-Kosterlitz-Thouless Transition in Ultrathin NbN Films

We probe the presence of long-range correlations in phase fluctuations by analyzing the higher-order spectrum of resistance fluctuations in ultrathin NbN superconducting films. The non-Gaussian component of resistance fluctuations is found to be sensitive to film thickness close to the transition, which allows us to distinguish between mean field and Berezinskii-Kosterlitz-Thouless (BKT) type superconducting transitions. The extent of non-Gaussianity was found to be bounded by the BKT and mean field transition temperatures and depends strongly on the roughness and structural inhomogeneity of the superconducting films. Our experiment outlines a novel fluctuation-based kinetic probe in detecting the nature of superconductivity in disordered low-dimensional materials.

A MatLAB Based Virtual Phantom for 2D Electrical Impedance Tomography (MatVP2DEIT): Studying the Medical Electrical Impedance Tomography Reconstruction in Computer

Electrical Impedance Tomography (EIT) is a computerized medical imaging technique which reconstructs the electrical impedance images of a domain under test from the boundary voltage-current data measured by an EIT electronic instrumentation using an image reconstruction algorithm. Being a computed tomography technique, EIT injects a constant current to the patient's body through the surface electrodes surrounding the domain to be imaged (Omega) and tries to calculate the spatial distribution of electrical conductivity or resistivity of the closed conducting domain using the potentials developed at the domain boundary (partial derivative Omega). Practical phantoms are essentially required to study, test and calibrate a medical EIT system for certifying the system before applying it on patients for diagnostic imaging. Therefore, the EIT phantoms are essentially required to generate boundary data for studying and assessing the instrumentation and inverse solvers a in EIT. For proper assessment of an inverse solver of a 2D EIT system, a perfect 2D practical phantom is required. As the practical phantoms are the assemblies of the objects with 3D geometries, the developing of a practical 2D-phantom is a great challenge and therefore, the boundary data generated from the practical phantoms with 3D geometry are found inappropriate for assessing a 2D inverse solver. Furthermore, the boundary data errors contributed by the instrumentation are also difficult to separate from the errors developed by the 3D phantoms. Hence, the errorless boundary data are found essential to assess the inverse solver in 2D EIT. In this direction, a MatLAB-based Virtual Phantom for 2D EIT (MatVP2DEIT) is developed to generate accurate boundary data for assessing the 2D-EIT inverse solvers and the image reconstruction accuracy. MatVP2DEIT is a MatLAB-based computer program which simulates a phantom in computer and generates the boundary potential data as the outputs by using the combinations of different phantom parameters as the inputs to the program. Phantom diameter, inhomogeneity geometry (shape, size and position), number of inhomogeneities, applied current magnitude, background resistivity, inhomogeneity resistivity all are set as the phantom variables which are provided as the input parameters to the MatVP2DEIT for simulating different phantom configurations. A constant current injection is simulated at the phantom boundary with different current injection protocols and boundary potential data are calculated. Boundary data sets are generated with different phantom configurations obtained with the different combinations of the phantom variables and the resistivity images are reconstructed using EIDORS. Boundary data of the virtual phantoms, containing inhomogeneities with complex geometries, are also generated for different current injection patterns using MatVP2DEIT and the resistivity imaging is studied. The effect of regularization method on the image reconstruction is also studied with the data generated by MatVP2DEIT. Resistivity images are evaluated by studying the resistivity parameters and contrast parameters estimated from the elemental resistivity profiles of the reconstructed phantom domain. Results show that the MatVP2DEIT generates accurate boundary data for different types of single or multiple objects which are efficient and accurate enough to reconstruct the resistivity images in EIDORS. The spatial resolution studies show that, the resistivity imaging conducted with the boundary data generated by MatVP2DEIT with 2048 elements, can reconstruct two circular inhomogeneities placed with a minimum distance (boundary to boundary) of 2 mm. It is also observed that, in MatVP2DEIT with 2048 elements, the boundary data generated for a phantom with a circular inhomogeneity of a diameter less than 7% of that of the phantom domain can produce resistivity images in EIDORS with a 1968 element mesh. Results also show that the MatVP2DEIT accurately generates the boundary data for neighbouring, opposite reference and trigonometric current patterns which are very suitable for resistivity reconstruction studies. MatVP2DEIT generated data are also found suitable for studying the effect of the different regularization methods on reconstruction process. Comparing the reconstructed image with an original geometry made in MatVP2DEIT, it would be easier to study the resistivity imaging procedures as well as the inverse solver performance. Using the proposed MatVP2DEIT software with modified domains, the cross sectional anatomy of a number of body parts can be simulated in PC and the impedance image reconstruction of human anatomy can be studied.