Accelerating frequency-domain diffuse optical tomographic image reconstruction using graphics processing units

Diffuse optical tomographic image reconstruction uses advanced numerical models that are computationally costly to be implemented in the real time. The graphics processing units (GPUs) offer desktop massive parallelization that can accelerate these computations. An open-source GPU-accelerated linear algebra library package is used to compute the most intensive matrix-matrix calculations and matrix decompositions that are used in solving the system of linear equations. These open-source functions were integrated into the existing frequency-domain diffuse optical image reconstruction algorithms to evaluate the acceleration capability of the GPUs (NVIDIA Tesla C 1060) with increasing reconstruction problem sizes. These studies indicate that single precision computations are sufficient for diffuse optical tomographic image reconstruction. The acceleration per iteration can be up to 40, using GPUs compared to traditional CPUs in case of three-dimensional reconstruction, where the reconstruction problem is more underdetermined, making the GPUs more attractive in the clinical settings. The current limitation of these GPUs in the available onboard memory (4 GB) that restricts the reconstruction of a large set of optical parameters, more than 13, 377. (C) 2010 Society of Photo-Optical Instrumentation Engineers. DOI: 10.1117/1.3506216]

Accelerated gradient based diffuse optical tomographic image reconstruction

Purpose: Fast reconstruction of interior optical parameter distribution using a new approach called Broyden-based model iterative image reconstruction (BMOBIIR) and adjoint Broyden-based MOBIIR (ABMOBIIR) of a tissue and a tissue mimicking phantom from boundary measurement data in diffuse optical tomography (DOT). Methods: DOT is a nonlinear and ill-posed inverse problem. Newton-based MOBIIR algorithm, which is generally used, requires repeated evaluation of the Jacobian which consumes bulk of the computation time for reconstruction. In this study, we propose a Broyden approach-based accelerated scheme for Jacobian computation and it is combined with conjugate gradient scheme (CGS) for fast reconstruction. The method makes explicit use of secant and adjoint information that can be obtained from forward solution of the diffusion equation. This approach reduces the computational time many fold by approximating the system Jacobian successively through low-rank updates. Results: Simulation studies have been carried out with single as well as multiple inhomogeneities. Algorithms are validated using an experimental study carried out on a pork tissue with fat acting as an inhomogeneity. The results obtained through the proposed BMOBIIR and ABMOBIIR approaches are compared with those of Newton-based MOBIIR algorithm. The mean squared error and execution time are used as metrics for comparing the results of reconstruction. Conclusions: We have shown through experimental and simulation studies that Broyden-based MOBIIR and adjoint Broyden-based methods are capable of reconstructing single as well as multiple inhomogeneities in tissue and a tissue-mimicking phantom. Broyden MOBIIR and adjoint Broyden MOBIIR methods are computationally simple and they result in much faster implementations because they avoid direct evaluation of Jacobian. The image reconstructions have been carried out with different initial values using Newton, Broyden, and adjoint Broyden approaches. These algorithms work well when the initial guess is close to the true solution. However, when initial guess is far away from true solution, Newton-based MOBIIR gives better reconstructed images. The proposed methods are found to be stable with noisy measurement data. (C) 2011 American Association of Physicists in Medicine. DOI: 10.1118/1.3531572]

Contribution of solvent effect to the optical spectra of polyanilines: has it been underestimated

The coefficient of thermochromism of polyaniline solutions has been found to be solvent dependent and the solvent effect is not negligible. Hence, thermochromism of polyaniline solutions cannot be explained solely on the basis of conformational change induced by a change in temperature. Further, comparison of the solvatochromism of polyaniline and polytoluidine shows a higher solvatochromic shift for the former. It implies that the higher energy associated with the exciton peak of polytoluidine is not due to the higher ring torsional angle induced by the higher steric repulsion of the methyl group, as widely accepted, but is due to its less solvatochromic red-shift as compared to polyaniline.

Holographic optical element based single-mode hybrid fiber optic interferometer for realizing zero-order fringe

A holographic optical element (HOE) based single-mode hybrid fiber optic interferometer for realizing the zero-order fringe is described. The HOE proposed and used integrates the actions of a beam combiner and a lens, and endows the interferometer with high tolerance for repositioning errors. The proposed method is simple and offers advantages such as the elimination of in situ processing for the hologram.

Optical and electrical properties of SnS semiconductor crystals grown by physical vapor deposition technique

Tin sulfide (SnS) is a material of interest for use as an absorber in low cost solar cells. Single crystals of SnS were grown by the physical vapor deposition technique. The grown crystals were characterized to evaluate the composition, structure, morphology, electrical and optical properties using appropriate techniques. The composition analysis indicated that the crystals were nearly stoichiometric with Sn-to-S atomic percent ratio of 1.02. Study of their morphology revealed the layered type growth mechanism with low surface roughness. The grown crystals had orthorhombic structure with (0 4 0) orientation. They exhibited an indirect optical band gap of 1.06 eV and direct band gap of 1.21 eV with high absorption coefficient (up to 10(3) cm(-1)) above the fundamental absorption edge. The grown crystals were of p-type with an electrical resistivity of 120 Omega cm and carrier concentration 1.52 x 10(15) cm(-3). Analysis of optical absorption and diffuse reflectance spectra showed the presence of a wide absorption band in the wavelength range 300-1200 nm, which closely matches with a significant part of solar radiation spectrum. The obtained results were discussed to assess the suitability of the SnS crystal for the fabrication of optoelectronic devices. (C) 2011 Elsevier B.V. All rights reserved.

A Non-contact measurement technique to measure micro surface stress and obtain deformation profiles of the order of 1nm in Micro-cantilever based structures by single image Optical Diffraction method

A new method based on analysis of a single diffraction pattern is proposed to measure deflections in micro-cantilever (MC) based sensor probes, achieving typical deflection resolutions of 1nm and surface stress changes of 50 mu N/m. The proposed method employs a double MC structure where the deflection of one of the micro-cantilevers relative to the other due to surface stress changes results in a linear shift of intensity maxima of the Fraunhofer diffraction pattern of the transilluminated MC. Measurement of such shifts in the intensity maxima of a particular order along the length of the structure can be done to an accuracy of 0.01mm leading to the proposed sensitivity of deflection measurement in a typical microcantilever. This method can overcome the fundamental measurement sensitivity limit set by diffraction and pointing stability of laser beam in the widely used Optical Beam Deflection method (OBDM).

Deblurring in a noncoherent optical processing system: pupil function synthesis and experimental implementation

The bipolar point spread function (PSF) corresponding to the Wiener filter tor correcting linear-motion-blurred pictures is implemented in a noncoherent optical processor. The following two approaches are taken for this implementation: (1) the PSF is modulated and biased so that the resulting function is non-negative and (2) the PSF is split into its positive and sign-reversed negative parts, and these two parts are dealt with separately. The phase problem associated with arriving at the pupil function from these modified PSFs is solved using both analytical and combined analytical-iterative techniques available in the literature. The designed pupil functions are experimentally implemented, and deblurring in a noncoherent processor is demonstrated. The postprocessing required (i.e., demodulation in the first approach to modulating the PSF and intensity subtraction in the second approach) are carried out either in a coherent processor or with the help of a PC-based vision system. The deblurred outputs are presented.

Diamond-like carbon coating for optical windows

Hard, low stress diamond-like carbon films have been deposited by plasma assisted chemical vapour deposition technique, The various substrates include soft IR components like ZnS and ZnSe windows, Gaseous precursors such as propene, ethyl alcohol and acetone have been used to synthesize the films to study the nature of precursors in determining the film compatibility with the underlying component (substrate), The residual compressive stresses, the Young's modulus and the adhesion energy of the films have been estimated to be 10(10) dynes/cm(2), 10(10) N/m(2) and 1000 ergs/cm(2) respectively. To alleviate film failure, a study on the effects of additive gases such as hydrogen and the use of buffer layers such as ZrO2, has been undertaken, The diamond-like carbon films produced here are hard (5000 kg/mm(2)), specularly smooth in the wavelength region from 2.5 mu m to 20 mu m, with no microstructural features and have excellent adhesion on ZnS and ZnSe windows. The figure of merit of these films for aero-space applications has been evaluated by subjecting the film-buffer layer ZnS or ZnSe composite stack to wind, dust and rain erosion studies and by establishing the integrity of the specular IR transmittance of the stack upto 16 or 20 mu m as the case may be.

Measurement of complex optical flow with use of an augmented generalized gradient scheme

We present a method for measuring the local velocities and first-order variations in velocities in a timevarying image. The scheme is an extension of the generalized gradient model that encompasses the local variation of velocity within a local patch of the image. Motion within a patch is analyzed in parallel by 42 different spatiotemporal filters derived from 6 linearly independent spatiotemporal kernels. No constraints are imposed on the image structure, and there is no need for smoothness constraints on the velocity field. The aperture problem does not arise so long as there is some two-dimensional structure in the patch being analyzed. Among the advantages of the scheme is that there is no requirement to calculate second or higher derivatives of the image function. This makes the scheme robust in the presence of noise. The spatiotemporal kernels are of simple form, involving Gaussian functions, and are biologically plausible receptive fields. The validity of the scheme is demonstrated by application to both synthetic and real video images sequences and by direct comparison with another recently published scheme Biol. Cybern. 63, 185 (1990)] for the measurement of complex optical flow.

Growth and characterization of some novel crystals for nonlinear optical applications

The advent of high intensity lasers coupled with the recent advances in crystal technology has led to rapid progress in the field of nonlinear optics. This article traces the history of materials development that has taken place over the past forty odd years and dwells on the current status in this important area. The materials aspect is discussed under three classes viz. inorganic, organic and semiorganic crystals. In the end, some of the crystal growth work that has been carried out in author's laboratory is presented.

Hydrogen-bond-directed self-assembly of D-(+)-dibenzoyltartaric acid and 4-aminopyridine: optical nonlinearities and stoichiometry-dependent novel structural features

Attempts to prepare hydrogen-bond-directed nonlinear optical materials from a 1:1 molar mixture Of D-(+)-dibenzoyltartaric acid (DBT, I) and 4-aminopyridine (4-AP, II) resulted in two salts of different stoichiometry. One of them crystallizes in an unusual 1.5:1 (acid:base) monohydrate salt form III while the other one crystallizes as 1:1 (acid:base) salt IV. Crystal structures of both of the salts were determined from single-crystal X-ray diffraction data. The salt III crystallizes in a monoclinic space group C2 with a = 30.339(l), b = 7.881(2), c = 14.355(1) angstrom, beta = 97.48(1)degrees, V = 3403.1(9) angstrom3, Z = 4, R(w) = 0.058, R(w)= 0.058. The salt IV also crystallizes in a monoclinic space group P2(1) with a = 7.500(1), b = 14.968(2), c = 10.370(1) angstrom, beta = 102.67(1)degrees, V = 1135.9(2) angstrom3, Z = 2, R = 0.043, R(w) = 0.043. Interestingly, two DBT molecules with distinctly different conformation are present in the same crystal lattice of salt III. Extensive hydrogen-bonding interactions are found in both of the salts, and both of them show SHG intensity 1.4-1.6 times that of urea.

Optical absorption and photoluminescence studies of Nd3 doped alkali boro germanate glasses

Optical absorption and photoluminescence studies have been carried out at room temperature in 25 R2O-25 GeO2-49.5 B2O3-0.5 Nd2O3 glass systems, (Composition in mol%, R= Li, Na, K and Rb). Judd Ofelt Intensity parameters and other parameters like Racah (E-1, E-2 and E-3), Slater-Condon-Shortley (F-2, F-4 and F-6) Spin-Orbit Coupling (xi(4f)) and Configuration Interaction (alpha,beta and gamma) for Nd3+ ion in the glass system are calculated. The variation of the 02 parameters are interpreted in terms of the covalency of the RE ion in the glass matrix. Further the hypersensitive transition I-4(9/2) -> (4)G(5/2), (2)G(7/2) is analyzed with respect to the intensity ratio I-L/I-S and is found to be dependent on the type of alkali in the glass matrix. The Photoluminescence studies do not show any appreciable shift in the peak emission wavelength of the F-4(3/2) to I-4(11/2) transition with the change in alkali type. (C) 2010 Elsevier B.V. All rights reserved.

Microscopic Mechanism of 1/f Noise in Graphene: Role of Energy Band Dispersion

A distinctive feature of single-layer graphene is the linearly dispersive energy bands, which in the case of multilayer graphene become parabolic. A simple electrical transport-based probe to differentiate between these two band structures will be immensely valuable, particularly when quantum Hall measurements are difficult, such as in chemically synthesized graphene nanoribbons. Here we show that the flicker noise, or the 1/f noise, in electrical resistance is a sensitive and robust probe to the band structure of graphene. At low temperatures, the dependence of noise magnitude on the carrier density was found to be opposite for the linear and parabolic bands. We explain our data with a comprehensive theoretical model that clarifies several puzzling issues concerning the microscopic origin of flicker noise in graphene field-effect transistors (GraFET).

Linear electro‐optic dispersion in (‐)‐2‐(α‐methylbenzylamino)‐5‐nitropyridine single crystals

The variation of the linear electro-optic effect in (-)-2-(alpha-methylbenzylamino)-5-nitropyridine with the wavelength of the incident light at room temperature has been measured. The reduced half-wave voltages have been found to have the values 2.1, 2.8, and 6.0 kV at 488, 514.5, and 632.8 nm respectively and the corresponding values of the linear electro-optic coefficient have been evaluated.;The interpretation of the results in terms of the structures of the molecule and the crystal is discussed. The thermal variation of the birefringence has also been investigated and the coefficient for the temperature variation of the refractive index difference is found to have the value (d Delta n/dT)=9.3X10(-5) K-1.

Measurement of complex optical flow with use of an augmented generalized gradient scheme

We present a method for measuring the local velocities and first-order variations in velocities in a time-varying image. The scheme is an extension of the generalized gradient model that encompasses the local variation of velocity within a local patch of the image. Motion within a patch is analyzed in parallel by 42 different spatiotemporal filters derived from 6 linearly independent spatiotemporal kernels. No constraints are imposed on the image structure, and there is no need for smoothness constraints on the velocity field. The aperture problem does not arise so long as there is some two-dimensional structure in the patch being analyzed. Among the advantages of the scheme is that there is no requirement to calculate second or higher derivatives of the image function. This makes the scheme robust in the presence of noise. The spatiotemporal kernels are of simple form, involving Gaussian functions, and are biologically plausible receptive fields. The validity of the scheme is demonstrated by application to both synthetic and real video images sequences and by direct comparison with another recently published scheme [Biol. Cybern. 63, 185 (1990)] for the measurement of complex optical flow.