Optical Conveyor Belts for Chiral Discrimination: Influence of De-Phasing Parameter

A numerical analysis is carried out of the influence of the de-phasing parameter of an optical conveyor belt in the enantiomeric separation. The optical conveyor belt is obtained by the interference of a Laguerre Gaussian and a Gaussian beam with different beam waists, which are temporally de-phased. In order to obtain the maximum separation distance between enantiomers, we calculate the optimum range of values of the de-phasing parameter.


Introduction
One important task in biology is the separation of enantiomers, which can be achieved by using different techniques such as gas chromatography, high performance liquid chromatography, or capillary electrophoresis [1].Different groups [2][3][4][5] have recently analyzed the use of optical fields to generate enantiomeric separation.In this sense, a study has been made of lateral optical forces (LOF) generated by interfering plane waves [6] in order to obtain chirality sorting, and the LOF effect on paired chiral particles has also been discussed [7].Furthermore, another way of generating LOF is given by evanescent waves, which transport spin angular momentum and could be used for chirality sorting [3].Thus, micro-metric particles have been sorted in a fluidic medium by using the spin-dependent optical forces of chiral particles [8].Moreover, our group has demonstrated that it is possible to obtain chiral resolution of nanoparticles by using an optical conveyor belt [9][10][11][12] through the superposition of two orthogonal counter-propagating temporally-de-phased beams [13,14].In this paper, we are going to analyze the effects of the temporal de-phasing parameter on the chiral resolution of the optical conveyor belt, showing that the maximum separation distance between enantiomers can be obtained by selecting the correct range of the de-phased parameter.Finally, the dynamics of enantiomers are studied in their final trap region.

Optical Forces
In order to analyze particle dynamics, we have solved the overdamped Langevin equation [15]: where F( R(t)) is the optical force, R(t) is the position vector of the particle at time t, −γ d R dt is the frictional force of a particle, and F (t) is a random function force with time that we have assumed null (the Brownian effect is not taken into account).Coefficient γ = 6πη v a, where η v is the viscosity of the medium (water in our case) and parameter a is the particle radius.We consider small chiral particles with radius a much lower than the wavelength of the electromagnetic field E, H, inside a dielectric medium of refractive index n.The optical force on the chiral sphere (CS) can be deduced by using the dipole approximation for which CS is modeled as an induced magnetic ( m) and electric dipole ( p) given by [16]: where H and E are the magnetic and electric fields, α ee , α mm , and α em are the electric, magnetic, and chiral polarizability, respectively, η is the vacuum impedance, and 0 is the permittivity of free space.
In the case of Rayleigh particles with a small chirality parameter δ (it is assumed that the refractive index of the chiral particle is given by n LR = np ± δ, where n L = np + δ defines the refractive index for the left-handed enantiomer and n R = np + δ for the right-handed one), the electric, magnetic, and chiral polarizability can be written (in absence of absorption) as [3,16]: These expressions have been obtained by taking into account that the chirality parameter |δ| << 1, which implies that in the dipolar limit, all magnetic parts (α mm is proportional to δ 2 ) of force can be disregarded [2].The electric energy density (U e ) and the helicity density (U h ) can be defined as: where and * denote the imaginary part and the complex conjugate respectively, r is the relative permittivity of the medium in which the particle is immersed, and k 0 is the wavenumber.Due to the electric density energy U e and helicity density U h , the particle is subject to a potential energy W given by: In the Rayleigh regime, the optical force F( R(t)), described by Equation ( 1), exerted upon the chiral spheres by the electromagnetic field, can be approximated to a gradient force [4,17]:

Electromagnetic Field
Let us consider two counter-propagating Gaussian and Laguerre-Gaussian beams (LGB) [13], whose electric fields ( e 0 , e 1 ) are elliptically polarized, slightly temporally de-phased (τ t), and with different beam waist: The total electric field may be written as: introducing cylindrical coordinates r = x 2 + y 2 and θ = arctan( y x ), the amplitude and phase of the fields ( e 0 , e 1 ) given by Equations ( 8) and ( 9) can be expressed (assuming a (1-1) Laguerre Gaussian beam) as: where is the wave number inside of a dielectric medium of refractive index n, and A 0 , w q (z), ζ q (z), and R q (z) are the electric field amplitude, beam waist, the Gouy phase, and the curvature radius of Gaussian wave fronts for q ∈ {0, 1}, respectively [18].These parameters are given by: where λ 0 is the vacuum wavelength.The magnetic field according to the Maxwell equations is given by: ω 0 being the electromagnetic field frequency.Under the condition b 0 = −b −1 1 , the polarization of the counter-propagating beams will be right-handed and left-handed, respectively, the degree of circular polarization being different from one if b 1 = 1.In this case, Equation ( 6) can be written as: where: Taking into account that φ 0 − φ 1 ≈ −2kz − θ, Equation (20) can therefore be written as: where φ(z, θ) = τ t − 2k z − θ.Equation ( 24) describes a chiral-dielectric optical trap given by W 0 (r, z) and a chiral conveyor given by W s (r, z) sin By introducing Equation (24) into Equations ( 7) and ( 1), we obtain the following dynamic equations for chiral particles: , so the dynamic equations of chiral particles can be approximated to: ) Equation (28) shows that chiral particles are trapped in the radial positions r p where the function 26) and (27) imply that chiral particles describe helical trajectories similar to those obtained in the case of dielectric particles [19].

Discussion
In the numerical calculation, we have used the values of λ = 1070 nm, a = λ/15, n p = 1.45, n = 1.33, w 0 = 5λ, w 1 = 4λ, b 1 = 0.9, δ = ±0.06,and A 0 = 2 × 10 8 V/m assuming that the initial position of the enantiomers is r 0 = 3λ, θ 0 = 0, z 0 = 10λ.Figure 1 shows the potential energy W(r, z) given by Equation (20).As can be observed, W shows a maximum at radial position r p = 1.6λ (green dashed line), which is the radial coordinate where particles will be initially trapped according to Equation (28).Thus, by using the indicated fixed values, the chiral conveyor is mainly controlled by the de-phase τ, which governs the conveyor velocity.In fact, if τ = 0, there will not be any temporal dependence in Equations ( 29) and (30), and there will not be a chiral conveyor; consequently, W will act like a chiral trap.In this case, both enantiomers will be trapped in the same region, as shown in Figure 2a, and as a result, there is no effective enantiomeric separation.Figure 2b-f shows the enantiomeric separation produced by the action of the chiral optical conveyor described by Equations ( 1), (7), and ( 25)-( 27).
As can be observed, the particle trajectories are spiral in all cases, the pitch and velocity of optical conveyors being strongly dependent on de-phase parameter τ [19].In this sense, the behavior of chiral conveyors is very similar to dielectric conveyors, and the main difference from the results obtained in [19] is that the movement of particles in the z direction is limited due to the axial dependence of parameters W 0 (r, z) and W s (r, z).Furthermore, as can be seen in Figure 2, on increasing τ, the enantiomeric separation also increases until it reaches a maximum separation that slowly decreases for high values of de-phase τ.This result can be clarified in Figure 3, where the behavior of enantiomeric separation (∆, defined as the axial distance between trapped enantiomers) is represented as a function of τ.As previously mentioned, there is a fast increment of ∆ when τ rises, obtaining a maximum value of 18λ for τ = 10 Hz, which remains nearly constant (see the inset graphics in Figure 3), until reaching a value of τ = 100 Hz to decrease slowly subsequently.Figure 4 shows the particles position at different times for a de-phase parameter τ = 10 Hz.As can be observed in Figure 4a, enantiomers are trapped at r = r p in different z positions in a short time (t = 0.01 s). Figure 4b-d shows that positive enantiomers (δ > 0) are trapped more quickly than negative ones.It can also be observed in Figure 4b-e that there are two traps, but at the end (t = 10 s), only one of them is stable for each enantiomer.It is interesting to note that when enantiomers are finally trapped at different z positions (see Figures 2 and 4f), they describe circular trajectories with angular frequencies ω, which are shown in Figure 5.As can be observed in Figure 5, at low and high τ values, the frequency of the circular trajectories ω shows similar values for both enantiomers (1 kHz < ω < 1 Hz).However, for τ values between 1 and 1000 Hz, the frequency ω is nearly constant and equal to 2 Hz for enantiomers with δ < 0 and 3 Hz in the case of δ > 0. One explanation for this result could be that, at the final time, particles describe a limit cycle in which z(t) reaches a limit value z τ ±δ (each enantiomer goes to a final z position that depends on their value ±δ and de-phase τ) at the fixed value r = r p (see Figure 2).In this case (when the limit cycle is reached z = z τ ±δ , r = r p ), Equations ( 28) and (30) can be approximated to: Introducing Equation (33) into (32), we finally obtain that: As can be observed in Figure 5, the approximated function (34) qualitatively explains the final frequency values ω, although quantitatively, the predictions are overestimated be approximately 1 Hz.
A possible explanation for this observed overestimation is that we have used the approximations: (∂ z W s (r, z) + 2kW c (r, z))sin(φ(z, θ)) ≈ 2kW c (r, z)sin(φ(z, θ)) (∂ z W c (r, z) − 2kW s (r, z))cos(φ(z, θ)) ≈ −2kW s (r, z)cos(φ(z, θ)) in order to obtain the analytical approximated Equation (34).Finally, Figure 6 shows the particle separation of particles with different chiral parameter δ when we use the de-phasing parameter τ = 10 Hz, taking in all cases a final time of 50 s.As can be observed, the particle separation is null for low values of δ, increasing for higher values of δ until reaching a saturation zone for a fixed time of 50 s.It is important to note that at this time, particles with δ > 0.7 have not reached their limit cycle.

Figure 3 .
Figure 3. Enantiomeric separation versus τ.The inset figure shows enantiomeric separation versus τ in the logarithmic scale.

Figure 5 .
Figure 5. Angular frequency of the last circular trajectory of trapped enantiomers (δ > 0 blue and δ < 0 red) versus τ (the continuous line shows the predicted results from Equation (34)).

Figure 6 .
Figure 6.Particle separation of particles with different chiral parameter δ, when we use the de-phasing parameter τ = 10 Hz, taking in all cases a final time of 50 s.