Single-resolution and multiresolution extended-Kalman-filter-based reconstruction approaches to optical refraction tomography

The problem of reconstruction of a refractive-index distribution (RID) in optical refraction tomography (ORT) with optical path-length difference (OPD) data is solved using two adaptive-estimation-based extended-Kalman-filter (EKF) approaches. First, a basic single-resolution EKF (SR-EKF) is applied to a state variable model describing the tomographic process, to estimate the RID of an optically transparent refracting object from noisy OPD data. The initialization of the biases and covariances corresponding to the state and measurement noise is discussed. The state and measurement noise biases and covariances are adaptively estimated. An EKF is then applied to the wavelet-transformed state variable model to yield a wavelet-based multiresolution EKF (MR-EKF) solution approach. To numerically validate the adaptive EKF approaches, we evaluate them with benchmark studies of standard stationary cases, where comparative results with commonly used efficient deterministic approaches can be obtained. Detailed reconstruction studies for the SR-EKF and two versions of the MR-EKF (with Haar and Daubechies-4 wavelets) compare well with those obtained from a typically used variant of the (deterministic) algebraic reconstruction technique, the average correction per projection method, thus establishing the capability of the EKF for ORT. To the best of our knowledge, the present work contains unique reconstruction studies encompassing the use of EKF for ORT in single-resolution and multiresolution formulations, and also in the use of adaptive estimation of the EKF's noise covariances. (C) 2010 Optical Society of America

Demonstration of a microfluidic polarimeter

Traditional methods of detecting chiral molecules, such as optical rotation are not suitable for miniaturization, since, the magnitude of the rotation of polarization scales down linearly with the optical path length of the device. Since the origin of optical activity is due to difference of refractive indices between the two circularly polarized states of light, it is possible to detect chiral media by measuring the dependence of the angles of refraction on the polarization state of the incident light. This however is a weak effect and hence requires sensitive optical detection schemes, based on novel polarization modulation techniques. The device can be scaled down for applications involving small sample volumes. Fabrication details of a prototype microfluidic device are described.

Membrane-based deformable mirror: intrinsic aberrations and alignment issues

A deformable mirror (DM) is an important component of an adaptive optics system. It is known that an on-axis spherical/parabolic optical component, placed at an angle to the incident beam introduces defocus as well as astigmatism in the image plane. Although the former can be compensated by changing the focal plane position, the latter cannot be removed by mere optical realignment. Since the DM is to be used to compensate a turbulence-induced curvature term in addition to other aberrations, it is necessary to determine the aberrations induced by such (curved DM surface) an optical element when placed at an angle (other than 0 deg) of incidence in the optical path. To this effect, we estimate to a first order the aberrations introduced by a DM as a function of the incidence angle and deformation of the DM surface. We record images using a simple setup in which the incident beam is reflected by a 37 channel micro-machined membrane deformable mirror for various angles of incidence. It is observed that astigmatism is a dominant aberration, which was determined by measuring the difference between the tangential and sagittal focal planes. We justify our results on the basis of theoretical simulations and discuss the feasibility of using such a system for adaptive optics considering a trade-off between wavefront correction and astigmatism due to deformation. (C) 2015 Optical Society of America

Localization of light in coherently amplifying random media

We derive and analyze the statistics of reflection coefficient of light backscattered coherently from an amplifying and disordered optical medium modeled by a spatially random refractive index having a uniform imaginary part in one dimension. We find enhancement of reflected intensity owing to a synergy between wave confinement by Anderson localization and coherent amplification by the active medium. This is not the same as that due to enhanced optical path lengths expected from photon diffusion in the random active medium. Our study is relevant to the physical realizability of a mirrorless laser by photon confinement due to Anderson localization.

Two-axis force sensing and control of a Reorientable scanning probe

This paper presents the design and implementation of a reorientable scanning probe that is capable of two-axis force sensing and control in the 2-D scanning (X-Z) plane. The probe is comprised of three major components, namely a compliant manipulator, laser measurement system, and magnetic actuation system. Control of the position and orientation of the probe tip is realized by means of magnetic actuation combined with a novel structural design. The design of the manipulator's compliance and that of the optical path of the laser measurement system together enable achieving sensitivity to lateral (X) forces that is nearly identical to that of normal (Z) forces. The achieved sensitivity ratio, of about 0.6, is significantly higher than that of conventional scanning probe systems. The theoretical bases for the structural design and the sensitivity of the two-axis force sensing system are presented. Subsequently, fabrication of the manipulator is described and the result of experimental evaluation of the scanning probe's features is discussed. The scanning probe is used to access the vertical and re-entrant features on the two sides of a cylindrical micropipette, which are subsequently scanned by regulating the lateral force of tip-sample interaction.

Conformational Stability and Intramolecular Hydrogen Bonding in 1,2-Ethanediol and 1,4-Butanediol

The gas-phase infrared spectra of 1,2-ED and 1,4-BD have been, recorded at three different temperatures using a multipass gas cell of 6 m optical path length. DFT calculation has also been carried out using 6-311++G** and aug-cc-pVDZ basis sets to look for the existence of intramolecular hydrogen bonding, in them from the red shift and infrared absorption intensity enhancement of the bonded O-H band compared to that of the free O-H hand. Equilibrium population analysis With 10 conformers of 1,2-ED and 1,4-BD at experimental temperatures were-carried out for the reconstruction of the Observed vibrational spectra at that temperature,using standard statistical relationships. The most abundant conformer at experimental temperatures, was identified. In 1,2-ED a red shift of 45 cm(-1) in the intramolecularly interacting O-H stretching vibrational band position and no significant intensity enhancement compared to that of-the free O-H have been observed. On the contrary, in one of the hydrogen bonded conformers of 1,4-BD, a 124 cm(-1) red shift in the O-H stretching frequency and 8.5 times-intensity enhancement for the ``bonded'' O-H compared to that of the ``free'' O-H is seen. On the basis of this comparative study, we have concluded that strong intramolecular hydrogen bonding exists in 1,4-BD. But there appears, to be weak intramolecular hydrogen bonding in 1,2-ED at temperatures of 303, 313, and 323 Km the gas phase We have found that most stable hydrogen-bonded conformers of 1,4-BD are less populated than some of the non-hydrogen-bonded conformers. Even for the 1,4-BD, the relative population of the g'GG'Gt conformer, which has a strong intramolecular,hydrogen bond, is less than what is predicted. Perhaps the intramolecular hydrogen bond plays a less Significant role in the relative stability,of the various Conformers than what has been predicted from calculations and prevails in the literature.

Non-interferometric methods of phase estimation for application in optical tomography

We describe here two non-interferometric methods for the estimation of the phase of transmitted wavefronts through refracting objects. The phase of the wavefronts obtained is used to reconstruct either the refractive index distribution of the objects or their contours. Refraction corrected reconstructions are obtained by the application of an iterative loop incorporating digital ray tracing for forward propagation and a modified filtered back projection (FBP) for reconstruction. The FBP is modified to take into account non-straight path propagation of light through the object. When the iteration stagnates, the difference between the projection data and an estimate of it obtained by ray tracing through the final reconstruction is reconstructed using a diffraction tomography algorithm. The reconstruction so obtained, viewed as a correction term, is added to the estimate of the object from the loop to obtain an improved final refractive index reconstruction.

Ultrasound Assisted Optical Tomography: Estimation of phase Shift experienced by photon on transit through US insonified region for detection of breast tumor

A Monte Carlo model of ultrasound modulation of multiply scattered coherent light in a highly scattering media has been carried out for estimating the phase shift experienced by a photon beam on its transit through US insonified region. The phase shift is related to the tissue stiffness, thereby opening an avenue for possible breast tumor detection. When the scattering centers in the tissue medium is exposed to a deterministic forcing with the help of a focused ultrasound (US) beam, due to the fact that US-induced oscillation is almost along particular direction, the direction defined by the transducer axis, the scattering events increase, thereby increasing the phase shift experienced by light that traverses through the medium. The phase shift is found to increase with increase in anisotropy g of the medium. However, as the size of the focused region which is the region of interest (ROI) increases, a large number of scattering events take place within the ROI, the ensemble average of the phase shift (Delta phi) becomes very close to zero. The phase of the individual photon is randomly distributed over 2 pi when the scattered photon path crosses a large number of ultrasound wavelengths in the focused region. This is true at high ultrasound frequency (1 MHz) when mean free path length of photon l(s) is comparable to wavelength of US beam. However, at much lower US frequencies (100 Hz), the wavelength of sound is orders of magnitude larger than l(s), and with a high value of g (g 0.9), there is a distinct measurable phase difference for the photon that traverses through the insonified region. Experiments are carried out for validation of simulation results.

Ray tracing based path-length calculations for polarized light tomographic imaging

A ray tracing based path length calculation is investigated for polarized light transport in a pixel space. Tomographic imaging using polarized light transport is promising for applications in optical projection tomography of small animal imaging and turbid media with low scattering. Polarized light transport through a medium can have complex effects due to interactions such as optical rotation of linearly polarized light, birefringence, diattenuation and interior refraction. Here we investigate the effects of refraction of polarized light in a non-scattering medium. This step is used to obtain the initial absorption estimate. This estimate can be used as prior in Monte Carlo (MC) program that simulates the transport of polarized light through a scattering medium to assist in faster convergence of the final estimate. The reflectance for p-polarized (parallel) and s-polarized (perpendicular) are different and hence there is a difference in the intensities that reach the detector end. The algorithm computes the length of the ray in each pixel along the refracted path and this is used to build the weight matrix. This weight matrix with corrected ray path length and the resultant intensity reaching the detector for each ray is used in the algebraic reconstruction (ART) method. The proposed method is tested with numerical phantoms for various noise levels. The refraction errors due to regions of different refractive index are discussed, the difference in intensities with polarization is considered. The improvements in reconstruction using the correction so applied is presented. This is achieved by tracking the path of the ray as well as the intensity of the ray as it traverses through the medium.

A critical analysis of the shear difference method

The shear difference method which is commonly used for the separation of normal stresses using photoelastic techniques depends on the step-by-step integration of one of the differential equations of equilibrium. It is assumed that the isoclinic and the isochromatic parameters measured by the conventional methods pertain to the state of stress at the midpoint of the light path. In practice, a slice thin enough for the above assumption to be true and at the same time thick enough to give differences in the shear-stress values over the thickness is necessary. The paper discusses the errors introduced in the isoclinic and isochromatic values by the conventional methods neglecting the variation of stresses along the light path. It is shown that while the error introduced in the measurement of the isochromatic parameter may not be serious, the error caused in the isoclinic measurement may lead to serious errors. Since the shear-difference method involves step-by-step integration the error introduced will be of a cumulative nature.

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.

Exact path-integral evaluation of the heat distribution function of a trapped Brownian oscillator

Using path integrals, we derive an exact expression-valid at all times t-for the distribution P(Q,t) of the heat fluctuations Q of a Brownian particle trapped in a stationary harmonic well. We find that P(Q, t) can be expressed in terms of a modified Bessel function of zeroth order that in the limit t > infinity exactly recovers the heat distribution function obtained recently by Imparato et al. Phys. Rev. E 76, 050101(R) (2007)] from the approximate solution to a Fokker-Planck equation. This long-time result is in very good agreement with experimental measurements carried out by the same group on the heat effects produced by single micron-sized polystyrene beads in a stationary optical trap. An earlier exact calculation of the heat distribution function of a trapped particle moving at a constant speed v was carried out by van Zon and Cohen Phys. Rev. E 69, 056121 (2004)]; however, this calculation does not provide an expression for P(Q, t) itself, but only its Fourier transform (which cannot be analytically inverted), nor can it be used to obtain P(Q, t) for the case v=0.

Analysis of the evanescent field coupling between a fiber and a planar waveguide using the explicit finite difference method with a nonuniform mesh (Journal Paper)

We investigate an optical waveguide system consisting of an unclad fiber core suspended at a constant distance parallel to the surface of a planar waveguide. The coupling and propagation of light in the combined system is studied using the three-dimensional explicit finite difference beam propagation method with a nonuniform mesh configuration. The power loss in the fiber and the field distribution in the waveguide are studied as a function of various parameters, such as index changes, index profile, and propagation distance, for the combined system.

Quality-of-Transmission-Aware Manycasting Over Optical Burst-Switched Networks

Many next-generation distributed applications, such as grid computing, require a single source to communicate with a group of destinations. Traditionally, such applications are implemented using multicast communication. A typical multicast session requires creating the shortest-path tree to a fixed number of destinations. The fundamental issue in multicasting data to a fixed set of destinations is receiver blocking. If one of the destinations is not reachable, the entire multicast request (say, grid task request) may fail. Manycasting is a generalized variation of multicasting that provides the freedom to choose the best subset of destinations from a larger set of candidate destinations. We propose an impairment-aware algorithm to provide manycasting service in the optical layer, specifically OBS. We compare the performance of our proposed manycasting algorithm with traditional multicasting and multicast with over provisioning. Our results show a significant improvement in the blocking probability by implementing optical-layer manycasting.

Influence of Raman scattering on the cross phase modulation in optical fibers

The effect of Raman scattering on co-propagation of two short optical pulses is considered. The intra pulse Raman scattering causes the self-frequency shift of each pulse. The effect of the inter pulse Raman scattering is to enhance the frequency shift while the stimulated Raman scattering (SRS) term suppresses (enhances) the frequency shift if the center frequency difference between the optical pulses falls to the right (left) of the Raman gain peak. An expression for the frequency shift as a function of the propagation distance is obtained.