Computational Models for New Fiber Optic Tweezers

ritbeiro@gmail.com Abstract: This paper discusses the calculation of the trapping forces in optical tweezers using a combination of the finite differences time domain (FDTD) method and the Lorentz force on electric dipoles. The results of 2D simulations of the trapping of a circular particle by a waveguide with a circular tip are presented and


Introduction
This paper presents a simulation method for computing the force generated by an electromagnetic field produced by a steering and focusing waveguide on small dielectric particles. This method combines the calculation of the stationary state of the electromagnetic field in a system composed of a waveguide and a dielectric particle, as shown in Fig. 1, using the electromagnetic equation propagation (MEEP) implementation of the finite differences time domain (FDTD) method, developed at the Massachusetts Institute of Technology (MIT), with the Lorentz equation for the force on an electric dipole. This method can be used to study and optimize trapping in optical tweezers (OT), for example, for applications in single cell manipulation in suspension or in optofluidic platforms [1].
With the advent of microfluidic and optofluidic technologies [2,3] it has become possible to manipulate single cells in a contactless minimally invasive fashion. This can be used to facilitate single cell biosensing, for example, single cell genomics, proteomics, and metabolomics [4]. structure (in grey on the left) and a circular dielectric particle (in grey on the right), both with a relative dielectric constant of 2.88: a point source particle is located along the optical axis of the waveguide at x=1.1, and the spatial units are in micrometers.
Most advanced OT systems rely on bulk elements comprising a trapping laser, beam expansion and steering optics, high numerical aperture (NA) objective, holder, and elements that allow dynamic control over the position and stiffness of the trapping force: typically acousto-optic deflectors or liquid crystals space light modulators (SLM) are used. Such setups are often incorporated in an inverted microscope adding the imaging capability. While these approaches enable advanced features such as multiple point holographic trapping and time multiplexing of trapping points, their implementation is often bulky and expensive, thus preventing their widespread use [5]. Furthermore, operation of bulk OT setups in turbid media or thick samples is often compromised by difficulties in focusing. This hinders the use of this technique in unprocessed samples or in vivo applications [1,5].
In this context, fiber optics based tweezers can provide a miniaturized, versatile solution in a relatively low cost setup. The majority of papers in this field described experimental results [6,7] with different designs of fiber optic OT. Even for other applications of OT, there were very few models and systematic studies which could be used to optimize the design and operation of OT. One of the few exceptions was based on the model of Mie for light scattering [6], however, this has been applied only to spherical and near spherical particles [7]. Therefore, it is necessary to develop new tools capable of addressing scatters with different and complex geometry, such as living cells, which will allow designing more effective fiber OT. This paper is organized as follows: Section 2 describes the theory of calculating the trapping force via the Lorentz force that acts on the dipoles that compose the scattering particle and the simulation method, Section 3 discusses the results, and Section 4 presents the conclusions.

Simulation method
There has been a strong debate in the literature of how to compute the force produced by an electromagnetic wave on a particle. Two main groups of approaches could be identified. A macroscopic approach based on either the Abraham or Minkowski formulation, considered respectively the kinetic and the canonical momentum in the presence of a scattering object. And also, a microscopic approach based on the Lorentz force which acts on each fundamental constituent of the scattering object [8,9].
The Lorentz theory can be used in two ways. Either it considers that the particle is composed of individual particles or microscopic dipoles [10,11]. At the microscopic scale, both these approaches lead to the same results, however, at a mesoscopic level the results differ marginally. This work considers only the approach based on the Lorentz force applied on dipoles [9], though the extensions to the other cases can be easily done.
In this case, the Lorentz force density per unit volume acting on the dipoles is where P is the microscopic polarization, E is the microscopic electric field, and B is the microscopic magnetic field. The microscopic Maxwell equations in the absence of currents and charges are where c is the speed of light in the vacuum.
Combining these equations and (1) yields [12] ( ) where the indexes denote the spatial components of the fields. When the field is stationary, the second term in (2) is null. Introducing the definition of the macroscopic polarization field into (2) allows to compute 2 0 where  is the relative permittivity of the scattering particle, and 0  is the vacuum permittivity.
For monochromatic radiation, it is useful to rewrite (4) in terms of the complex macroscopic electric field, resulting in 2 0 Equation (5) expresses the force density per unit volume at each point of the particle in terms of the stationary electric field used to trap it. Integrating (5) over the entire volume of the particle gives the net force acting on it.
For low intensity light beams, it is necessary to take into account the reflections of light at the surface of the particle, which was done using the MEEP simulation software. MEEP is an open source program, developed at MIT, which uses the FDTD method [13]. This is widely used to simulate electromagnetic wave propagation. This algorithm is based on the Yee grid discretization to stagger the electric and magnetic field in time and space. A grid with a perfectly matched layer (PML) was implemented to minimize the effects of reflections at the boundary.
The computational grid had 12 per 4 units, with a resolution of 16 -1 and a PML thickness of 0.5 units. The simulation considered a 2D system composed of the long waveguide with a circular tip and a circular particle, as seen in Fig. 1. The surrounding medium had a relative dielectric constant of 1, and both the particle and the waveguide had a dielectric constant of 2.88. A continuous wave point source with z-polarization was located in the middle of the waveguide. Figure 2 presents the stationary amplitude distribution of the electric field obtained from the simulations for different positions of the dielectric particles. When the particle is placed exactly in front of the beam, it produces a decrease in the electric field amplitude on the right side of the particle (identifiable in the figure) which results in a momentum transfer from the beam to the particle, producing a force along the optical axis. However, if the particle is placed off the axes, the incident beam is deflected on the scatter, and it produces a kick on the particle in the opposite direction. Both cases exemplify the mechanism of momentum transfer from the beam to the particle, indicating the model validity.  Fig. 3 Field of forces acting on the beam for particle radii respectively of 0.25 micrometers (a) and 0.5 micrometers (b) (the calculations were done using the normalized field intensity, the amplitude of the force is in arbitrary units and will scale linearly with the power of the light source, and the spatial units are in micrometers). Figure 3 represents the field of forces acting on the particle as it moves in front of the beam. Each arrow represents the force acting on the particle if it is placed in that point of space. It can be noticed that, near optical axis of the beam, the particles are attracted to the optical beam axis trapping them while, far from the axis, they are pushed forward. The structure of the field of forces depends on the size of the particle which determines the respective scattering cross sections as well as the profiles of the beam that exits the waveguide.

Conclusions
This paper presents a method of calculating the forces produced on a dielectric particle by an electromagnetic wave. This method is based on a combination of an FDTD solver (implemented in MEEP) of the electromagnetic field and a microscopic model for the Lorentz force acting on each dipole which constitutes the dielectric particle. Using MEEP takes into account the changes in the field produced by the scattering of light through the particle, thus improving the results relative to previous studies [7]. The preliminary results indicated the suitability of this method to be used in the computation of the force for arbitrary shapes of particles becoming a powerful tool and to optimize the design of new fiber optic tweezers.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.