Singular value decomposition based computationally efficient algorithm for rapid dynamic near-infrared diffuse optical tomography

Purpose: A computationally efficient algorithm (linear iterative type) based on singular value decomposition (SVD) of the Jacobian has been developed that can be used in rapid dynamic near-infrared (NIR) diffuse optical tomography. Methods: Numerical and experimental studies have been conducted to prove the computational efficacy of this SVD-based algorithm over conventional optical image reconstruction algorithms. Results: These studies indicate that the performance of linear iterative algorithms in terms of contrast recovery (quantitation of optical images) is better compared to nonlinear iterative (conventional) algorithms, provided the initial guess is close to the actual solution. The nonlinear algorithms can provide better quality images compared to the linear iterative type algorithms. Moreover, the analytical and numerical equivalence of the SVD-based algorithm to linear iterative algorithms was also established as a part of this work. It is also demonstrated that the SVD-based image reconstruction typically requires O(NN2) operations per iteration, as contrasted with linear and nonlinear iterative methods that, respectively, requir O(NN3) and O(NN6) operations, with ``NN'' being the number of unknown parameters in the optical image reconstruction procedure. Conclusions: This SVD-based computationally efficient algorithm can make the integration of image reconstruction procedure with the data acquisition feasible, in turn making the rapid dynamic NIR tomography viable in the clinic to continuously monitor hemodynamic changes in the tissue pathophysiology.

Transient response of nonlinear magneto-optic rotation in a paraffin-coated Rb vapor cell

We study resonant nonlinear magneto-optic rotation (NMOR) in a paraffin-coated Rb vapor cell as the magnetic field is swept. At low sweep rates, the nonlinear rotation appears as a narrow resonance signal with a linewidth of about ``300 mu G''(2 pi x 420 Hz). At high sweep rates, the signal shows transient response with an oscillatory decay. The decay time constant is of order 100 ms. The behavior is different for transitions starting from the lower or the upper hyperfine level of the ground state because of optical pumping effects.

Synthesis of Diketopyrrolopyrrole Containing Copolymers: A Study of Their Optical and Photovoltaic Properties

The diketopyrrolopyrrole-based copolymers PDPP-BBT and TDPP-BBT were synthesized and used as donor for bulk heterojunction photovoltaic devices. The photophysical properties of these polymers showed absorption in the range 500-600 nm with a maximum peak around 563 nm, while TDPP-BBT showed broadband absorption in the range 620 - 800 nm with a peak around 656 nm. The power conversion efficiencies (PCE) of the polymer solar cells based on these copolymers and [6,6]-phenyl C61 butyric acid methyl ester (PCBM) were 0.68% (as cast PDPP-BBT:PCBM), 1.51% (annealed PDPP-BBT:PCBM), 1.57% (as cast TDPPBBT: PCBM), and 2.78% (annealed TDPP-BBT:PCBM), under illumination of AM 1.5 (100 mW/cm2). The higher PCE for TDPP-BBT-based polymer solar cells has been attributed to the low band gap of this copolymer as compared to PDPP-BBT, which increases the numbers of photogenerated excitons and corresponding photocurrent of the device. These results indicate that PDPP-BBT and TDPP-BBT act as excellent electron donors for bulk heterojunction devices.

Photophysics of conjugated polymers: interplay between Forster energy migration and defect concentration in shaping a photochemical funnel in PPV

Recent single molecule experiments have suggested the existence of a photochemical funnel in the photophysics of conjugated polymers, like poly[2-methoxy-5-(2'-ethylhexyl)oxy-1,4-phenylenevinylene] (MEH-PPV). The funnel is believed to be a consequence of the presence of conformational or chemical defects along the polymer chain and efficient non-radiative energy transfer among different chromophore segments. Here we address the effect of the excitation energy dynamics on the photophysics of PPV. The PPV chain is modeled as a polymer with the length distribution of chromophores given either by a Gaussian or by a Poisson distribution. We observe that the Poisson distribution of the segment lengths explains the photophysics of PPV better than the Gaussian distribution. A recently proposed version of an extended particle-in-a-box' model is used to calculate the exciton energies and the transition dipole moments of the chromophores, and a master equation to describe the excitation energy transfer among different chromophores. The rate of energy transfer is assumed to be given here, as a first approximation, by the well-known Forster expression. The observed excitation population dynamics confirms the photochemical funneling of excitation energy from shorter to longer chromophores of the polymer chain. The time scale of spectral shift and energy transfer for our model polymer, with realistic values of optical parameters, is in the range of 200-300 ps. We find that the excitation energy may not always migrate towards the longest chromophore segments in the polymer chain as the efficiency of energy transfer between chromophores depends on the separation distance between the two and their relative orientation.

Optical nonlinearity in lead sulfide microtowers

Lead sulfide (PbS) microtowers on silicon substrates, having the physical properties of bulk PbS, have been synthesized. Optical nonlinearity studies using the open aperture z-scan technique employing 5 ns and 100 fs laser pulses reveal effective two-photon type absorption. For nanosecond excitation the nonlinear absorption coefficients (beta(eff)) are in the order of 10(-11) m W-1, two orders of magnitude less than the values reported for quantum confined PbS nanocrystals. For femtosecond excitation beta(eff) is of the order of 10(-14) m W-1. These results obtained in bulk PbS experimentally confirm the importance of quantum confinement in the enhancement of optical nonlinearities in semiconductor materials.

Electronic excitations and alternation of conjugated polymers

The electronic excitations and fluorescence of conjugated polymers are related to large or small alternation ? of the transfer integrals t(1 ± ?) along the backbone. The fluorescence of polysilanes (PSs) and poly (para-phenylenevinylene (PPV) is linked to large ?, which places the one-photon gap Eg below the lowest two-photon gap Ea and reduces distortions due to electron-phonon (e-p) coupling. In contrast to small ? not, vert, similar 0.1 in ?-conjugated polymers, such as polyacetylene (PA), para-conjugated phenyls lead to an extended ?-system with increased alternation, to states localized on each ring and to charge-transfer excitations between them. Surprisingly good agreement is found between semiempirical parametric method 3 (PM3) bond lengths and exact Pariser-Parr-Pople (PPP) ?-bond orders for trans-stilbene, where the PPV bipolarons are confined to two phenyls. Stilbene spectra are consistent with increased alternation and small e-p distortions.

Spectral approach for the soliton and periodic solutions of the nonlinear wave equation

A spectral method that obtains the soliton and periodic solutions to the nonlinear wave equation is presented. The results show that the nonlinear group velocity is a function of the frequency shift as well as of the soliton power. When the frequency shift is a function of time, a solution in terms of the Jacobian elliptic function is obtained. This solution is periodic in nature, and, to generate such an optical pulse train, one must simultaneously amplitude- and frequency-modulate the optical carrier. Finally, we extend the method to include the effect of self-steepening.

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]

Work fluctuations in an elastic dumbbell model of polymers in planar elongational flow

We use a path-integral approach to calculate the distribution P(w, t) of the fluctuations in the work W at time t of a polymer molecule (modeled as an elastic dumbbell in a viscous solvent) that is acted on by an elongational flow field having a flow rate (gamma) over dot. We find that P(w, t) is non-Gaussian and that, at long times, the ratio P(w, t)/ P (-w, t) is equal to expw/(k(B)T)], independent of (gamma) over dot. On the basis of this finding, we suggest that polymers in elongational flows satisfy a fluctuation theorem.

Synthesis and characterization of photo-crosslinkable main-chain liquid-crystalline polymers containing bis(benzylidene)cycloalkanone units

A series of new photo-crosslinkable main-chain liquid-crystalline polymers containing bis(benzylidene)cycloalkanone units have been studied. These units in the polymers function as mesogens as well as photoactive centres. Polyesters with three different bis(4-hydroxybenzylidene)cycloalkanones corresponding to three cycloalkanones, namely cyclopentanone, cyclohexanone and cycloheptanone, have been prepared. Three dicarboxylic acids with ether linkages, which were derived from oligoethylene oxides, namely triethylene glycol, tetraethylene glycol and pentaethylene glycol, have been used as spacers in these polymers. Polymerization was carried out by both solution and interfacial polycondensation; the latter method gave high-molecular-weight polymers. Structural characterizations were done by ultra-violet, infra-red and H-1 nuclear magnetic resonance spectroscopy. Liquid-crystalline properties were studied by differential scanning calorimetry and polarized-light optical microscopy. These polymers show a nematic mesophase. Liquid-crystalline transition temperatures were correlated with polymer structure. The decrease in transition temperature with increase in cycloalkanone ring size was explained in terms of the change in geometrical anisotropy of bis(benzylidene)cycloalkanone units. MNDO (modified neglect of differential overlap) calculations were performed on the model compounds, bis(4-acetyloxybenzylidene)cycloalkanone to elucidate the geometrical variation of the mesogenic units with cycloalkanone ring size. Studies of photolysis reveal the two kinds of photoreactions that proceed in these polymer systems, namely photoisomerization and photo-crosslinking. The former reaction disrupts the parallel stacking of the chromophores and is reflected as an increase in the ultra-violet spectral intensity. The favourability of these two reactions depends on the mobility of the polymer chains. When the photolysis was done below T-g, photo-crosslinking dominates over photoisomerization. Above T-g, photoisomerization is followed by photo-crosslinking. The photosensitivity of the polymers decreases with increase in size of the cycloalkanone ring.

Artifact-free reconstruction from off-axis digital holograms through nonlinear filtering

We present experimental investigation of a new reconstruction method for off-axis digital holographic microscopy (DHM). This method effectively suppresses the object auto-correlation, commonly called the zero-order term, from holographic measurements, thereby suppressing the artifacts generated by the intensities of the two beams employed for interference from complex wavefield reconstruction. The algorithm is based on non-linear filtering, and can be applied to standard DHM setups, with realistic recording conditions. We study the applicability of the technique under different experimental configurations, such as topographic images of microscopic specimens or speckle holograms.

Optical and Electronic Properties of Conjugated Polymer -Nanocluster Semiconductor Hybrid Systems

Conjugated polymers are intensively pursued as candidate materials for emission and detection devices with the optical range of interest determined by the chemical structure. On the other hand the optical range for emission and detection can also be tuned by size selection in semiconductor nanoclusters. The mechanisms for charge generation and separation upon optical excitation, and light emission are different for these systems. Hybrid systems based on these different class of materials reveal interesting electronic and optical properties and add further insight into the individual characteristics of the different components. Multilayer structures and blends of these materials on different substrates were prepared for absorption, photocurrent (Iph), photoluminescence (PL) and electroluminscence (EL) studies. Polymers chosen were derivatives of polythiophene (PT) and polyparaphenylenevinylene (PPV) along with nanoclusters of cadmium sulphide of average size 4.4 nm (CdS-44). The photocurrent spectral response in these systems followed the absorption response around the band edges for each of the components and revealed additional features, which depended on bias voltage, thickness of the layers and interfacial effects. The current-voltage curves showed multi-component features with emission varying for different regimes of voltage. The emission spectral response revealed additive features and is discussed in terms of excitonic mechanisms.

Optical and structural properties of TiO2 films deposited by Sol-Gel technique

TiO2 thin films were prepared by sol gel method. The structural investigations performed by means of X- ray diffraction (XRD) technique, Scanning electronic microscopy (SEM) showed the shape structure at T=600°C. The optical constants of the deposited film were obtained from the analysis of the experimental recorded transmittance spectral data over the wavelengths range 200-3000 nm. The values of some important parameters (refractive index n, dielectric constant ε ∞ and thickness d), and the third order optical nonlinear susceptibility χ(3) of TiO2 film are determined from these spectra. It has been found that the dispersion data obey the single oscillator relation of the Wemple-DiDomenico model, from which the dispersion parameters and high – frequency dielectric constant were determined. The estimation of the corresponding band gap Eg , χ (3) and ε ∞ are 2.57 eV, 0.021 × 10-10 esu and 5.20,respectively.

Ultrasound-modulated optical tomography: recovery of amplitude of vibration in the insonified region from boundary measurement of light correlation

We address a certain inverse problem in ultrasound-modulated optical tomography: the recovery of the amplitude of vibration of scatterers [p(r)] in the ultrasound focal volume in a diffusive object from boundary measurement of the modulation depth (M) of the amplitude autocorrelation of light [phi(r, tau)] traversing through it. Since M is dependent on the stiffness of the material, this is the precursor to elasticity imaging. The propagation of phi(r, tau) is described by a diffusion equation from which we have derived a nonlinear perturbation equation connecting p(r) and refractive index modulation [Delta n(r)] in the region of interest to M measured on the boundary. The nonlinear perturbation equation and its approximate linear counterpart are solved for the recovery of p(r). The numerical results reveal regions of different stiffness, proving that the present method recovers p(r) with reasonable quantitative accuracy and spatial resolution. (C) 2011 Optical Society of America

Linearization and Reconstruction of Non-linear Diffuse Optical Tomographic Image

Diffuse optical tomography (DOT) is one of the ways to probe highly scattering media such as tissue using low-energy near infra-red light (NIR) to reconstruct a map of the optical property distribution. The interaction of the photons in biological tissue is a non-linear process and the phton transport through the tissue is modelled using diffusion theory. The inversion problem is often solved through iterative methods based on nonlinear optimization for the minimization of a data-model misfit function. The solution of the non-linear problem can be improved by modeling and optimizing the cost functional. The cost functional is f(x) = x(T)Ax - b(T)x + c and after minimization, the cost functional reduces to Ax = b. The spatial distribution of optical parameter can be obtained by solving the above equation iteratively for x. As the problem is non-linear, ill-posed and ill-conditioned, there will be an error or correction term for x at each iteration. A linearization strategy is proposed for the solution of the nonlinear ill-posed inverse problem by linear combination of system matrix and error in solution. By propagating the error (e) information (obtained from previous iteration) to the minimization function f(x), we can rewrite the minimization function as f(x; e) = (x + e)(T) A(x + e) - b(T)(x + e) + c. The revised cost functional is f(x; e) = f(x) + e(T)Ae. The self guided spatial weighted prior (e(T)Ae) error (e, error in estimating x) information along the principal nodes facilitates a well resolved dominant solution over the region of interest. The local minimization reduces the spreading of inclusion and removes the side lobes, thereby improving the contrast, localization and resolution of reconstructed image which has not been possible with conventional linear and regularization algorithm.