Simultaneous multilayer scanning and detection for multiphoton fluorescence microscopy

Fast three-dimensional (3D) imaging requires parallel optical slicing of a specimen with an efficient detection scheme. The generation of multiple localized dot-like excitation structures solves the problem of simultaneous slicing multiple specimen layers, but an efficient detection scheme is necessary. Confocal theta detection (detection at 90 degrees to the optical axis) provides a suitable detection platform that is capable of cross-talk-free fluorescence detection from each nanodot (axial dimension approximate to 150 nm). Additionally, this technique has the unique feature of imaging a specimen at a large working distance with super-resolution capabilities. Polarization studies show distinct field structures for fixed and fluid samples, indicating a non-negligible field-dipole interaction. The realization of the proposed imaging technique will advance and diversify multiphoton fluorescence microscopy for numerous applications in nanobioimaging and optical engineering.

An inverse problem for Helmholtz equation

This article is concerned with subsurface material identification for the 2-D Helmholtz equation. The algorithm is iterative in nature. It assumes an initial guess for the unknown function and obtains corrections to the guessed value. It linearizes the otherwise nonlinear problem around the background field. The background field is the field variable generated using the guessed value of the unknown function at each iteration. Numerical results indicate that the algorithm can recover a close estimate of the unknown function based on the measurements collected at the boundary.

A novel variable resolution global spectral method on the sphere

A variable resolution global spectral method is created on the sphere using High resolution Tropical Belt Transformation (HTBT). HTBT belongs to a class of map called reparametrisation maps. HTBT parametrisation of the sphere generates a clustering of points in the entire tropical belt; the density of the grid point distribution decreases smoothly in the domain outside the tropics. This variable resolution method creates finer resolution in the tropics and coarser resolution at the poles. The use of FFT procedure and Gaussian quadrature for the spectral computations retains the numerical efficiency available with the standard global spectral method. Accuracy of the method for meteorological computations are demonstrated by solving Helmholtz equation and non-divergent barotropic vorticity equation on the sphere. (C) 2011 Elsevier Inc. All rights reserved.

A numerical study of variable coefficient elliptic Cauchy problem via projection method

In this paper, we investigate a numerical method for the solution of an inverse problem of recovering lacking data on some part of the boundary of a domain from the Cauchy data on other part for a variable coefficient elliptic Cauchy problem. In the process, the Cauchy problem is transformed into the problem of solving a compact linear operator equation. As a remedy to the ill-posedness of the problem, we use a projection method which allows regularization solely by discretization. The discretization level plays the role of regularization parameter in the case of projection method. The balancing principle is used for the choice of an appropriate discretization level. Several numerical examples show that the method produces a stable good approximate solution.

High resolution multiple excitation spot optical microscopy

We propose fundamental improvements in three-dimensional (3D) resolution of multiple excitation spot optical microscopy. The excitation point spread function (PSF) is generated by two interfering counter-propagating depth-of-focus beams along the optical axis. Detection PSF is obtained by coherently interfering the emitted fluorescent light (collected by both the objectives) at the detector. System PSF shows upto 14-fold reduction in focal volume as compared to confocal, and almost 2-fold improvement in lateral resolution. Proposed PSF has the ability to simultaneously excite multiple 3D-spots of sub-femtoliter volume. Potential applications are in fluorescence microscopy and nanobioimaging. Copyright 2011 Author(s). This article is distributed under a Creative Commons Attribution 3.0 Unported License. [doi:10.1063/1.3598413]

High resolution facies record on late Holocene flood plain sediments from lower reaches of Narmada Valley, Western India

A high resolution quantitative granulometric record for site Uchediya 21A degrees 43'2.22aEuro(3) N, 73A degrees 6'26.22aEuro(3) E; 10 m a. s. l.] gives understanding towards accretion history of the late Holocene flood plain in the lower reaches of Narmada River. Two sediment facies (sandy and muddy) and seven subfacies (sandy subfacies: St(MS+FS+CS), SmFS+MS, Sl(FS+VFS), and St(MS + CS); muddy subfacies: FmSILT+VFS+FS, FmSILT+VFS (O) and FmSILT+VFS (T)) are identified based on cluster analysis supplemented with sedimentary structures observed in field and other laboratory data. Changes in hydrodynamics are further deduced based on various sedimentological parameters and their ratios leading to arrive at a depositional model.

Finite Element Method For A Nonlocal Problem Of Kirchhoff Type

The nonlocal term in the nonlinear equations of Kirchhoff type causes difficulties when the equation is solved numerically by using the Newton-Raphson method. This is because the Jacobian of the Newton-Raphson method is full. In this article, the finite element system is replaced by an equivalent system for which the Jacobian is sparse. We derive quasi-optimal error estimates for the finite element method and demonstrate the results with numerical experiments.

A Numerical Scheme for Three-Dimensional Front Propagation and Control of Jordan Mode

As an example of a front propagation, we study the propagation of a three-dimensional nonlinear wavefront into a polytropic gas in a uniform state and at rest. The successive positions and geometry of the wavefront are obtained by solving the conservation form of equations of a weakly nonlinear ray theory. The proposed set of equations forms a weakly hyperbolic system of seven conservation laws with an additional vector constraint, each of whose components is a divergence-free condition. This constraint is an involution for the system of conservation laws, and it is termed a geometric solenoidal constraint. The analysis of a Cauchy problem for the linearized system shows that when this constraint is satisfied initially, the solution does not exhibit any Jordan mode. For the numerical simulation of the conservation laws we employ a high resolution central scheme. The second order accuracy of the scheme is achieved by using MUSCL-type reconstructions and Runge-Kutta time discretizations. A constrained transport-type technique is used to enforce the geometric solenoidal constraint. The results of several numerical experiments are presented, which confirm the efficiency and robustness of the proposed numerical method and the control of the Jordan mode.

Diffusion ordered spectroscopy for resolution of double bonded cis, trans-isomers

NMR spectroscopic separation of double bonded cis- and trans-isomers, that have different molecular shapes but identical mass have been carried out using Diffusion Ordered Spectroscopy (DOSY). The mixtures of fumaric acid and maleic acid, that have similar hydrodynamic radii, have resolved been on the basis of their diffusion coefficients arising due to their different tendencies to associate with micelles or reverse micelles. Sodium dodecyl sulfate (SDS) and Dioctyl sulfosuccinate sodium salt (AOT) have been used as the media to mimic the chromatographic conditions, modify the average mobility and to achieve differential diffusion rates. The best separation of the components has been achieved by Dioctyl sulfosuccinate sodium salt (AOT) in D2O solution. (C) 2012 Elsevier B.V. All rights reserved.

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.

An invariant subspace-based approach to the random eigenvalue problem of systems with clustered spectrum

The repeated or closely spaced eigenvalues and corresponding eigenvectors of a matrix are usually very sensitive to a perturbation of the matrix, which makes capturing the behavior of these eigenpairs very difficult. Similar difficulty is encountered in solving the random eigenvalue problem when a matrix with random elements has a set of clustered eigenvalues in its mean. In addition, the methods to solve the random eigenvalue problem often differ in characterizing the problem, which leads to different interpretations of the solution. Thus, the solutions obtained from different methods become mathematically incomparable. These two issues, the difficulty of solving and the non-unique characterization, are addressed here. A different approach is used where instead of tracking a few individual eigenpairs, the corresponding invariant subspace is tracked. The spectral stochastic finite element method is used for analysis, where the polynomial chaos expansion is used to represent the random eigenvalues and eigenvectors. However, the main concept of tracking the invariant subspace remains mostly independent of any such representation. The approach is successfully implemented in response prediction of a system with repeated natural frequencies. It is found that tracking only an invariant subspace could be sufficient to build a modal-based reduced-order model of the system. Copyright (C) 2012 John Wiley & Sons, Ltd.

High resolution methyl selective C-13-NMR of proteins in solution and solid state

New C-13-detected NMR experiments have been devised for molecules in solution and solid state, which provide chemical shift correlations of methyl groups with high resolution, selectivity and sensitivity. The experiments achieve selective methyl detection by exploiting the one bond J-coupling between the C-13-methyl nucleus and its directly attached C-13 spin in a molecule. In proteins such correlations edit the C-13-resonances of different methyl containing residues into distinct spectral regions yielding a high resolution spectrum. This has a range of applications as exemplified for different systems such as large proteins, intrinsically disordered polypeptides and proteins with a paramagnetic centre.

Data-resolution based optimization of the data-collection strategy for near infrared diffuse optical tomography

Purpose: To optimize the data-collection strategy for diffuse optical tomography and to obtain a set of independent measurements among the total measurements using the model based data-resolution matrix characteristics. Methods: The data-resolution matrix is computed based on the sensitivity matrix and the regularization scheme used in the reconstruction procedure by matching the predicted data with the actual one. The diagonal values of data-resolution matrix show the importance of a particular measurement and the magnitude of off-diagonal entries shows the dependence among measurements. Based on the closeness of diagonal value magnitude to off-diagonal entries, the independent measurements choice is made. The reconstruction results obtained using all measurements were compared to the ones obtained using only independent measurements in both numerical and experimental phantom cases. The traditional singular value analysis was also performed to compare the results obtained using the proposed method. Results: The results indicate that choosing only independent measurements based on data-resolution matrix characteristics for the image reconstruction does not compromise the reconstructed image quality significantly, in turn reduces the data-collection time associated with the procedure. When the same number of measurements (equivalent to independent ones) are chosen at random, the reconstruction results were having poor quality with major boundary artifacts. The number of independent measurements obtained using data-resolution matrix analysis is much higher compared to that obtained using the singular value analysis. Conclusions: The data-resolution matrix analysis is able to provide the high level of optimization needed for effective data-collection in diffuse optical imaging. The analysis itself is independent of noise characteristics in the data, resulting in an universal framework to characterize and optimize a given data-collection strategy. (C) 2012 American Association of Physicists in Medicine. http://dx.doi.org/10.1118/1.4736820]

Fast image reconstruction for fluorescence microscopy

Real-time image reconstruction is essential for improving the temporal resolution of fluorescence microscopy. A number of unavoidable processes such as, optical aberration, noise and scattering degrade image quality, thereby making image reconstruction an ill-posed problem. Maximum likelihood is an attractive technique for data reconstruction especially when the problem is ill-posed. Iterative nature of the maximum likelihood technique eludes real-time imaging. Here we propose and demonstrate a compute unified device architecture (CUDA) based fast computing engine for real-time 3D fluorescence imaging. A maximum performance boost of 210x is reported. Easy availability of powerful computing engines is a boon and may accelerate to realize real-time 3D fluorescence imaging. Copyright 2012 Author(s). This article is distributed under a Creative Commons Attribution 3.0 Unported License. http://dx.doi.org/10.1063/1.4754604]