Computational modelling of optical tweezers with many degrees of freedom using dynamic simulation: cylinders, nanowires, and multiple particles

Computational tasks such as the calculation and characterization of the optical force acting on a sphere are relatively straightforward in a Gaussian beam trap. Resulting properties of the trap such as the trap strength, spring constants, and equilibrium position can be easily determined. More complex systems with non-spherical particles or multiple particles add many more degrees of freedom to the problem. Extension of the simple methods used for single spherical particles could result in required computational time of months or years. Thus, alternative methods must be used. One powerful tool is to use dynamic simulation: model the dynamics and motion of a particle or particles within the trap. We demonstrate the use of dynamic simulation for non-spherical particles and multi-particle systems. Using a hybrid discrete dipole approximation (DDA) and T-matrix method, we find plausible equilibrium positions and orientations of cylinders of varying size and aspect ratio. Orientation landscapes revealing different regimes of behaviour for micro-cylinders and nanowires with different refractive indices trapped with beams of differing polarization are also presented. This investigation provides a solid background in both the function and properties of micro-cylinders and nanowires trapped in optical tweezers. This method can also be applied to particles with other shapes. We also investigate multiple-particle trapping, which is quite different from single particle systems, as they can include effects such as optical binding. We show that equilibrium positions, and the strength of interactions between particles can be found in systems of two and more particles.