Search

Lab-on-a-chip (LOC) devices capable of manipulating micro/nano-sized samples have spurred advances in biotechnology and chemistry. Designing and analyzing new and more advanced LOCs require accurate modeling and simulation of sample/particle dynamics inside such devices. In this work, we present a generalized computational physics model to simulate particle/sample trajectories under the influence of dielectrophoretic or optical forces inside LOC devices. The model takes into account time varying applied forces, Brownian motion, fluid flow, collision mechanics, and hindered diffusion caused by hydrodynamic interactions. We develop a numerical solver incorporating the aforementioned physics and use it to simulate two example cases: first, an optical trapping experiment, and second, a dielectrophoretic cell sorter device. In both cases, the numerical results are found to be consistent with experimental observations, thus proving the generality of the model. The numerical solver can simulate time evolution of the positions and velocities of an arbitrarily large number of particles simultaneously. This allows us to characterize and optimize a wide range of LOCs. The developed numerical solver is made freely available through a GitHub repository so that researchers can use it to develop and simulate new designs.

The motion of a nanoparticle in the vicinity of a near-field optical trap is modeled using the Fokker-Planck equation. A plasmonic C-shaped engraving on a gold film is considered as the optical trap. The time evolution of the position probability density of the nanoparticle is calculated to analyze the trapping dynamics. A spatially varying diffusion tensor is used in the formulation to take into account the hydrodynamic interactions. The steady-state position distribution obtained from the Fokker-Planck equation is compared with experimental results and found to be in good agreement. Computational cost of the proposed method is compared with the conventionally used Langevin equation based approach. The proposed method is found to be computationally efficient (requiring 35 times less computation time) and scalable to more complex lab-on-a-chip systems.

The optical forces generated by a right-handed plasmonic Archimedean spiral (PAS) have been mapped and analyzed. By changing the handedness of the circularly polarized excitation, the structure can switch from a trapping force profile to a rotating force profile. The Helmholtz-Hodge decomposition method has been used to separate the solenoidal component and the conservative component of the force and quantify their relative magnitude. It is shown that the for right-hand circularly polarized excitation, the PAS creates a significant amount of solenoidal forces. Using the decomposed force components, an intuitive explanation of the motion of micro- and nanoparticles in the force field is presented. Vector field topology is used to visualize the force components. The analysis is found to be consistent with numerical and experimental results. Due to the intuitive nature of the analysis, it can be used in the initial design process of complex laboratory-on-a-chip systems where a rigorous analysis is computationally expensive.

The force-field generated by a near-field optical trap is analyzed. A C-shaped engraving on a gold film is considered as the trap. By separating out the conservative component and the solenoidal component of the force-field using Helmholtz-Hodge decomposition, it was found that the force is non-conservative. Conventional method of calculating the optical potential from the force-field is shown to be inaccurate when the trapping force is not purely conservative. An alternative method is presented to accurately estimate the potential. The positional statistics of a trapped nanoparticle in this non-conservative field is calculated. A model is proposed that relates the position distribution to the conservative component of the force. The model is found to be consistent with numerical and experimental results. In order to show the generality of the approach, the same analysis is repeated for a plasmonic trap consisting of a gold nanopillar. Similar consistency is observed for this structure as well.

CsBr photocathodes have 10 times higher quantum efficiency with only 3 times larger intrinsic transverse emittance than copper. They are robust and can withstand 80 MV/m fields without breaking down or emitting dark current. They can operate in 2×10^{−9} torr vacuum and survive exposure to air. They are well suited for generating high pulse charge in rf guns without a photocathode transfer system.

Conference name: Proceedings of SPIE, Three-Dimensional Imaging, Visualization, and Display 2011 Date of Conference: 27–28 April 2011 The class of two-dimensional non-separable linear canonical transforms is the most general family of linear canonical transforms, which are important in both signal/image processing and optics. Application areas include noise filtering, image encryption, design and analysis of ABCD systems, etc. To facilitate these applications, one need to obtain a digital computation method and a fast algorithm to calculate the input-output relationships of these transforms. We derive an algorithm of NlogN time, N being the space-bandwidth product. The algorithm controls the space-bandwidth products, to achieve information theoretically sufficient, but not redundant, sampling required for the reconstruction of the underlying continuous functions. © 2011 SPIE.

A fast and accurate algorithm is developed for the numerical computation of the family of complex linear canonical transforms (CLCTs), which represent the input-output relationship of complex quadratic-phase systems. Allowing the linear canonical transform parameters to be complex numbers makes it possible to represent paraxial optical systems that involve complex parameters. These include lossy systems such as Gaussian apertures, Gaussian ducts, or complex graded-index media, as well as lossless thin lenses and sections of free space and any arbitrary combinations of them. Complex-ordered fractional Fourier transforms (CFRTs) are a special case of CLCTs, and therefore a fast and accurate algorithm to compute CFRTs is included as a special case of the presented algorithm. The algorithm is based on decomposition of an arbitrary CLCT matrix into real and complex chirp multiplications and Fourier transforms. The samples of the output are obtained from the samples of the input in approximately N log N time, where N is the number of input samples. A space-bandwidth product tracking formalism is developed to ensure that the number of samples is information-theoretically sufficient to reconstruct the continuous transform, but not unnecessarily redundant.

We report a fast and accurate algorithm for numerical computation of two-dimensional non-separable linear canonical transforms (2D-NS-LCTs). Also known as quadratic-phase integrals, this class of integral transforms represents a broad class of optical systems including Fresnel propagation in free space, propagation in graded-index media, passage through thin lenses, and arbitrary concatenations of any number of these, including anamorphic/astigmatic/non-orthogonal cases. The general two-dimensional non-separable case poses several challenges which do not exist in the one-dimensional case and the separable two-dimensional case. The algorithm takes approximately N log N time, where N is the two-dimensional space-bandwidth product of the signal. Our method properly tracks and controls the space-bandwidth products in two dimensions, in order to achieve information theoretically sufficient, but not wastefully redundant, sampling required for the reconstruction of the underlying continuous functions at any stage of the algorithm. Additionally, we provide an alternative definition of general 2D-NS-LCTs that shows its kernel explicitly in terms of its ten parameters, and relate these parameters bidirectionally to conventional ABCD matrix parameters.