Abstract
1. INTRODUCTION
Separation of enantiomers from racemic mixtures is essential in fields such as chemistry and pharmaceutics since the required performance is only exhibited by an enantiomer with a certain handedness (or chirality). Usually, methods based on chemical interactions, such as chiral high-performance liquid chromatography (HPLC), are employed to separate enantiomers starting from racemic mixtures. However, such chemical methods are usually slow, expensive, and molecule-dependent [1]. An interesting alternative would be the use of optical chiral forces [2]: since the chiral force exerted by light changes its sign when acting on enantiomers of different handedness [3], it could ultimately lead to light-driven enantiomeric separation overcoming many of the limitations of chemical methods.
There have been many recent theoretical and simulation works addressing the separation of enantiomers using light [4–10]. Indeed, there have been several experiments demonstrating optically induced separation, though for relatively large nanoparticles and nanostructures [11–13]. In all cases, free-space propagating beams are used, either being reflected at dielectric interfaces [5] or interfering with other beams [7] to produce the transverse optical spin that generates the required forces. Alternatively, one may think of using guided light for separation purposes, taking advantage of the enhancement of chiral interaction because the light is confined in subwavelength cross-sections over relatively long distances (ideally infinite for a lossless waveguide). One possibility is the use of optical nanofibers with cylindrical cross-sections, as recently noticed by Golat
In this work, we analyze the chiral separation properties of the simplest photonic integrated structure: a lossless dielectric waveguide with a rectangular cross-section on a lower-index substrate. We consider silicon nitride (SiN) as the material to build the waveguide core since it is transparent from telecom wavelengths down to the ultraviolet. The refractive index of SiN is large enough to ensure tight guiding when the core lies on a silicon dioxide substrate and is surrounded by water. Remarkably, waveguides with very low propagation loss () can be fabricated using mature tools and processes [20]. We first calculate the order of magnitude of the required chiral forces to perform the enantiomeric separation of particles under a reasonable time assuming Brownian diffusion of the target particles in a liquid environment. Then, we calculate numerically the electric and magnetic fields at different wavelengths from 405 to 1310 nm to obtain the optical forces using well-established equations [14]. We show that at short wavelengths (405 nm) lateral chiral forces arising from the transversal spin of the guided quasi-TE mode [15] can overcome the achiral forces and be used for separation of chiral particles of 80 nm radius. At longer wavelengths (1310 nm), this lateral force becomes much smaller but we combine the quasi-TE and quasi-TM modes of the waveguide with a proper 90° phase shift between them to generate a quasi-circularly polarized (quasi-CP) mode [16]. This mode produces a large transversal chiral gradient force that could separate nanoparticles of 52 nm radius and potentially molecules with a size of the order of 1 nm. Our results suggest that long SiN waveguides are simple but realistic structures towards achieving enantiomeric separation within seconds using guided light in PICs.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
2. OPTICAL FORCES EXERTED ON SMALL CHIRAL PARTICLES
The electromagnetic field of light carries momentum that can be transferred to a particle through the action of an optical force, and consequently, cause its motion. In this work, we restrict to study the motion of small chiral particles (whose size is smaller than the wavelength of light) subjected to optical forces. A small particle is fully characterized by its electric dipole moment, , which can be thought of as the separation of positive and negative charges, and its magnetic dipole moment, , which represents the overall current loop within the particle. The time-averaged force that the electromagnetic field exerts on a small particle is [5,14,21,22]
In this expression, and are, respectively, the electric and the magnetic fields at the position of the particle, is the impedance of the surrounding medium, is the electric permittivity of the medium, is the magnetic permeability of the medium, is the wavenumber, and is the wavelength of light. The dipole moments of a particle arise due to its interaction with the electric and magnetic fields of light, and are obtained as follows:
With Eqs. (1)–(4) we are ready to compute the time-averaged optical force due to a monochromatic electromagnetic field on a small particle. However, we gain more insight by developing the expression of the force from Eq. (1). The force can be split into several terms that depend on the following six time-averaged local field properties: electric energy density , magnetic energy density , helicity density , electric spin density , magnetic spin density , and complex Poynting vector :
The fully developed expression of the optical force acting on a particle is shown below, which has been split into the chiral and achiral terms, depending on whether the terms are a function of the chiral polarizability or not, respectively [14]:
3. OPTICALLY DRIVEN MOTION OF SMALL CHIRAL PARTICLES IN A FLUID
To study under what circumstances the optical enantioseparation is possible, we first need to examine the motion of particles in a fluid under the influence of an external chiral optical force field. To this end, we consider the following assumptions for this system: there is no net fluid flow, and the mass of the particles is negligible so that the viscous forces dominate the inertial forces. Under these considerations, the motion of the particle is driven by a combination of the external optical force, , the friction or drag force due to the viscosity of the fluid (which is opposite to the movement of the particle), and the force arising from the stochastic collisions of the smaller fluid molecules with the particle (Brownian motion). The variation of the particle position, , with time is governed by the overdamped Langevin equation [28–30]
The magnitude of the chiral optical forces required for enantiomeric separation and the needed sorting time can be estimated upon further assumptions. For that we examine the movement of a cloud of particles under the influence of an external optical force field in the bulk, i.e., in an infinite system with no boundaries. We followed the derivation from Kravets
Figure 1.(a) Schematic showing the combined action of the translation of the cloud of enantiomers due to the optical force (
In the integrated waveguide system, the particles flow along a microfluidic channel, and the force field is generated by a waveguide. Thus, the actual force field is not spatially uniform and there are boundaries that enclose the region of space where the particles can move and that modify the mobility coefficient value. However, Eq. (9) provides a quite accurate estimation of the sorting time for a longitudinally invariant dielectric waveguide-microchannel system. We assume the initial extension of the racemic mixture cloud is . This is achievable with the well-known microfluidics technique called hydrodynamic flow focusing [34], where two lateral flows can control the width of the middle fluid channel where the particles are suspended. The fluid is composed mainly of water at temperature so the dynamic viscosity is [28].
Figure 1(b) shows the sorting time for particles of different radii, calculated with Eq. (9) for , which marks the starting point of separation. Stronger forces are needed to sort larger particles for the same value of sorting time. Moreover, larger particles take longer to be sorted for the same value of the force. We emphasize that this sorting time is obtained for a force field that is constant in time and uniform in space, which is not the actual situation of the force field generated by a waveguide system. In the latter case, the forces are stronger at distances closer to the waveguide and decay with the distance due to the evanescent field of the mode in the fluid. The usefulness of this graph is to know what range of optical force magnitude the waveguide system needs to generate to sort particles under a reasonable time: 1 ms to few hours. We must therefore look for integrated waveguides that generate optical chiral forces within the range of to for sorting particles of radii between 1 and 1000 nm. Once we have designed those waveguide systems, we use the particle tracking algorithm [Eq. (8)] to test the actual enantiomeric separation that our waveguides can reach.
4. DESCRIPTION OF THE INTEGRATED PHOTONIC WAVEGUIDE
We consider a photonic strip waveguide made with a SiN core (refractive index ) on top of a substrate, and surrounded by water as the system that produces the force field responsible for the enantiomeric sorting, as shown in Fig. 2(a). The use of SiN has several practical advantages such as transparency at visible and near-infrared wavelengths, a relatively large refractive index to ensure tight localization of the fields in the waveguide core, and its processing with a mature silicon technology to produce low-loss waveguides [20]. Moreover, SiN is particularly appropriate for applications requiring immersion in fluid, such as photonic biosensing [35].
Figure 2.(a) Schematic of the waveguide cross-section representing the opposite action of the total optical force exerted by the guided mode (in red) onto particles with opposed chirality. (b) Transversal electric field intensity (in color map) and polarization (in arrow or ellipse map) of the quasi-TE mode and of the quasi-CP mode. The handedness of the polarization ellipses does not change its sign throughout the cross-section. (c) Net chiral and achiral force along the
In order to obtain the electric and magnetic fields of the guided modes, the Maxwell’s equations are solved in the waveguide using the finite element method implemented by the FemSIM solver in the commercial software RSoft (Synopsis). The software computes the eigenmode of the cross-section of the waveguide system, which is assumed to be invariant along the longitudinal direction (translational symmetry along the optical axis, i.e., –axis). The resulting electric and magnetic fields are plugged into Eq. (6) to obtain the optical forces per amount of power guided by the mode. We assume 20 mW of power in our simulations, which, despite being a high power level for integrated optics, can be attained using a standard semiconductor continuous wave laser without causing material damage. Noticeably, the forces depend on the characteristics of the particle, which are modeled by the polarizabilities [Eqs. (3) and (4)]. The particle’s material is modeled with a relative permittivity of , relative permeability , and chirality parameter . These values have been widely used to characterize chiral nanoparticles in the literature [5,36]. The medium (water) is modeled with the values and . We take into account the following two design criteria for the waveguide system towards favoring enantiomeric separation: (i) , to ensure that chiral forces dominate achiral ones (although it is not a necessary condition for achieving separation); and (ii) to , according to our calculations from Fig. 1.
After the total optical force field is computed, we use the particle tracking algorithm to track the position of an individual particle for each enantiomer for 2 s throughout a hypothetical microfluidic channel surrounding the waveguide. Due to the stochastic nature of the Brownian motion, we repeat this tracking 500 times to do a statistical analysis of the enantioseparation process. In addition, the initial positions of the particles were randomized to account for the realistic scenario of not being able to precisely control their starting positions within the microchannel. The initial –position was uniformly randomized at , which is achievable with hydrodynamic flow focusing as commented before [34], and the initial –position was uniformly randomized throughout the microchannel height. From the final position of the particles we calculate the enantiomer fraction (EF) for each enantiomeric cloud. More details about the particle tracking algorithm and statistical analysis can be found in Appendix A.
5. RESULTS
We consider two different approaches to sorting chiral nanoparticles throughout the transversal plane (–plane) around the waveguide: the fundamental quasi-transverse electric mode (quasi-TE mode) for horizontal sorting and a quasi-circularly polarized mode (quasi-CP mode) for attractive-repulsive sorting. The electric field intensity and polarization of the quasi-TE mode and the quasi-CP mode are plotted in Fig. 2(b). A more detailed decomposition of the electric and magnetic fields for the quasi-TE, quasi-TM, and quasi-CP modes can be found in Appendix C.
A. Quasi-TE Mode
The fundamental quasi-TE mode (or ) in a strip waveguide is characterized by the transverse electric field mainly pointing along the horizontal plane (–axis) and the transverse magnetic field mainly polarized along the vertical plane (along –axis). This guided mode exhibits a non-zero longitudinal component of the electric and magnetic fields due to confining the wave inside a waveguide. The longitudinal component oscillates out-of-phase in comparison with the transversal components. This results in a transverse spin [15], which can be interpreted as the quantum spin Hall effect of light [37] and gives rise to transversal chiral forces.
The optical forces that a quasi-TE mode exerts on a particle are computed for a wavelength of 405 nm and a particle radius of 80 nm. The field properties that are responsible for the dominant transversal forces in this system are the electric energy density , the magnetic energy density , and the magnetic spin . These field properties, together with their respective forces, are plotted in Fig. 3. The electric gradient force attracts any particle toward the sidewalls of the waveguide, whereas the magnetic gradient force repels any particle from the top of the waveguide. The magnetic spin force moves the -particles towards the left and the -particles towards the right. Over the top of the waveguide, the achiral gradient electric and gradient magnetic forces have opposite directions, thus reducing the strength of the total achiral force. That is why the chiral force magnitude () is of the same order as the total achiral force (), resulting in the total force along the –axis changing sign for opposite enantiomers over the top of the waveguide, thus, pushing enantiomers to opposite sides along the –direction. The motion of 500 particles per enantiomer is tracked individually for 2 s given the total optical force field shown in Fig. 3. The microchannel dimensions () are . The resulting final positions of the particles are plotted in Fig. 3 for both enantiomers. From among the 500 ()-particles 66.2% end up on the left side () and 65.6% of -particles end up on the right side (), yielding an enantiomer fraction (EF) of 65.8% and 66.0%, respectively. The ()-EF is calculated within the region and the -EF is calculated within the region .
Figure 3.Field properties (
A way to increase the strength of chiral forces over the achiral forces is by exploiting the wavelength dependency of the forces. Some of the forces depend on the ratio , as previously discussed for Eq. (6). This dependency is studied in Fig. 2(c) for the quasi-TE mode, where the total chiral and achiral forces along –direction, which is the sorting direction, are represented with respect to the particle radius for three wavelengths (405 nm, 780 nm, and 1310 nm). The calculated forces are expected to be more accurate for particles whose size is smaller than the wavelength: due to the dipolar approximation assumption [38]. This upper bound is , 290 nm, and 490 nm for , 780 nm, and 1310 nm, respectively. In all the cases, the forces are evaluated at a point positioned at a vertical distance equal to the particle’s radius over the top of the waveguide and at a horizontal distance equal to 1/4 of the waveguide’s width from the center, as shown in the inset in Fig. 2(e). That vertical distance is the minimum distance the particle can be placed due to its size. Figure 2(c) shows that there is a particle size range over which chiral forces become larger than achiral forces along the –direction: 65–132 nm for , 123–257 nm for , and 204–378 nm for . In fact, there is a specific radius that maximizes the ratio of chiral force over achiral force within those intervals. In these ranges, the spin magnetic force is stronger than the other achiral forces along the –direction. However, for smaller radii, the achiral magnetic gradient becomes dominant, and for larger radii, the sum of the achiral flow force due to the imaginary part of the Poynting vector and the achiral magnetic gradient becomes dominant instead. In addition, these intervals suggest that operating at shorter wavelengths is more suitable for sorting smaller chiral particles.
Another way to facilitate the horizontal chiral separation is to lower the strength of the achiral gradient forces along the – axis. This can be achieved by making the waveguide wider. The power is thus spread over a larger area, thereby reducing the gradient of the fields and their respective achiral gradient forces along the –axis. In addition, the smallest particle size for which chiral and achiral forces have the same magnitude reduces down to 28 nm for . This comes at the expense of reducing the strength of the chiral optical force too (), as shown in Fig. 4 for a waveguide , since the power of the mode has been distributed over a larger cross-section. The particle tracking simulation inside a microchannel of the same dimensions as before () yielded values of EF of 60.6%/60.2% for the -enantiomers, which are somewhat smaller when compared with the narrower waveguide despite the magnitude of the forces being one order of magnitude less. This might be because the wider waveguide allows the interaction between the particle and optical force for a longer time because the waveguide width is larger. In spite of yielding similar enantiomeric separation capability, this wider configuration should be easier to implement experimentally due to the larger area of interaction between the mode and the particles.
Figure 4.Total optical force (transversal in arrow map, and longitudinal in colormap) that a quasi-TE mode in a wide strip waveguide (
B. Quasi-CP Mode
We refer to the quasi-CP mode as a guided mode that is obtained by the superposition of the mode and the mode delayed by a phase shift of 90°. As the electric field is predominantly horizontally polarized in the mode and predominantly vertically polarized in the mode, the combination originates a guided mode with an effective circular polarization and, therefore, local helicity [16]. The waveguide width and thickness are chosen so that the mode and the mode are degenerate, i.e., both modes exhibit the same effective refractive index () at the target wavelength. This degeneracy allows the circular polarization of the mode to be maintained along the waveguide.
We compute the dominant achiral and chiral forces of this system for a particle of 52 nm radius and a wavelength of 1310 nm (see Fig. 5). The field of the quasi-CP mode exhibits an intrinsic local helicity , which naturally diminishes along the evanescent field beyond the waveguide core. This is a robust way of achieving a helicity gradient, which is in fact the predominant chiral force in this system. This force is particularly interesting because it depends only on the chiral polarizability, which means that it grows with the radius of the particle as instead of , and does not depend explicitly on the wavelength. The total optical force exhibits opposite directions for opposite chiralities, being attractive towards the waveguide for the ()-particles (with value ) and repulsive for the (–)-particles (with value ).
Figure 5.Field properties (
In order to identify the possible range of sizes that might be easier to sort with the quasi-CP compound mode, we performed the study of the chiral force’s dependence on the particle’s size at the wavelength of 1310 nm. As shown in Fig. 2(d), the net transversal chiral force (along – and –axis) dominates the net transversal achiral force, up to a maximum size (378 nm) where the combination of the forces due to the real and imaginary parts of the Poynting vector and the rotational of the spin overtake the helicity gradient force. Crucially, this means that the size range has no minimum radii where the chiral force is smaller than the achiral force, as it occurs for the quasi-TE mode. This means that we can expect to reach enantioseparation for smaller particles when compared to the TE mode, even down to radii around 1 nm (molecular size), making this force the most promising for sorting molecules. The particle tracking results inside a microchannel () for the system in Fig. 5 show that 33.2% of the 500 ()-particles become “attached” to the waveguide within a radius of 425 nm measured from the center of the waveguide, and that 90.2% of the 500 (–)-particles are repelled. The ()-EF within a radius of 425 nm is 77.2% and the (–)-EF outside the radius of 425 nm is 57.5%. Particle tracking simulations assuming optical power above 20 mW, and therefore stronger forces, suggest that the ()-EF and (–)-EF increase with power injected into the mode yielding 94.5% and 77.9% for 50 mW and 99.6% and 97.3% for 100 mW.
Enantioseparating forces depend on the chirality parameter , as shown in Eq. (6), which we have considered to be . However, the chiral parameter might be much lower for realistic particles such as molecules. We have computed the modulus of the chiral and achiral forces depending on for both the quasi-TE [Fig. 6(a)] and quasi-CP [Fig. 6(b)] modes. The results suggest that particles with () can be sorted under quasi-TE (quasi-CP) mode operation. Other configurations are needed to achieve the separation of particles with lower chirality. This will be explored in future work.
Figure 6.Net chiral and achiral forces exerted on a particle depending on its chirality parameter
6. CONCLUSIONS
In conclusion, we have identified an opportunity for transversal enantioseparation via optical forces in photonic integrated waveguides of high-chirality particles utilizing distinct mechanisms depending on particle size and operating wavelength. To this end, we have used the most straightforward waveguide structure: a strip SiN waveguide placed on a silica substrate. At short wavelengths (405 nm), the spin magnetic force arising from a quasi-TE mode is strong enough to sort particles of 80 nm radius in less than 2 s. At longer wavelengths (1310 nm), the helicity gradient force stemming from a quasi-CP mode can separate particles of radius as small as 52 nm radius under time spans below 2 s. In contrast to other approaches using optical waveguides [17–19], our waveguides are longitudinally invariant, meaning that the optical chiral forces could be exerted over long distances (cm-scale), thus facilitating practical enantioseparation with realistic optical powers in the chip (20 mW). These findings highlight the potential of optical forces generated in integrated waveguides in facilitating enantioseparation within the specified parameters of high chirality, particle size, and wavelength.
Appendix A: PARTICLE TRACKING ALGORITHM
The motion for a Brownian particle immersed in a fluid and subjected to a (optical) force field is governed by the Langevin equation [
This equation can be integrated by iteratively calculating the displacements over successive time intervals . If the time step is larger than the momentum relaxation time , i.e., , the previous equation simplifies to the so-called overdamped Langevin equation [
Our choice of fulfils this requirement. The mass of the particles is so small that in practice, one can consider that they lack inertia and the velocity is constant over . The upper bound on the time step is chosen so that the particle moves one-pixel distance as maximum. In our case, the spatial resolution was , and the maximum value of the force for the quasi-TE mode was 93 fN (20 mW) for an 80 nm particle; therefore, this upper bound is
The mobility coefficients are modified from the bulk value to account for the hydrodynamic interaction between the particle and the enclosing boundaries. The modifications are different depending on whether the particle moves along the directions parallel () or perpendicular () to a non-slip planar boundary [
The force field used in the algorithm is obtained for a system that is not enclosed by a microchannel. However, we assume the presence of the microchannel does not modify the guided modes, and thus, the force field produced by the waveguide throughout the surrounding medium (water, ). This approximation is valid due to the higher refractive index of the waveguide core (SiN, ) with respect to that of the microchannel material (, ), which ensures the guidance of the mode along the core, and because the microchannel walls (boundaries) are separated by a distance from the waveguide walls.
Once the particle tracking simulation is finished, we do statistical analysis with the last positions of the 500 particles for both enantiomers, to compute the enantiomer fraction. Given the number of ()-particles, , and the number of (–)-particles, , inside a region of space, we define the enantiomer fraction (EF) as [
Since the chiral forces separate opposite enantiomers in opposite directions, the ()-EF and the (–)-EF are calculated in different regions of space. For instance, for the quasi-TE mode, the ()-EF is calculated for the region , whereas the (–)-EF is obtained for the region . And, for the quasi-CP mode, the ()-EF is calculated for the region whereas the (–)-EF is obtained for , where is an arbitrary radius that defines a circular region from the center of the waveguide.
APPENDIX B: SORTING TIME FOR A CLOUD OF ENANTIOMERS IN A FLUID
We derive the expression for calculating the time needed for a chiral optical force, , to separate two clouds of opposite enantiomers with a distance . This has been previously done by Kravets
This equation can be solved as a second-degree equation on the variable , whose solution is
Note that no dipolar approximation is needed to derive this equation since we do not require the expression of the force.
APPENDIX C: ELECTRIC AND MAGNETIC FIELDS OF THE GUIDED MODES
The decomposition of the electric and magnetic fields for the quasi-TE, quasi-TM, and quasi-CP modes can be found in Fig.
Figure 7.Electric and magnetic field patterns of the guided modes discussed in the main text for a strip waveguide. The fields are split into the transversal (
References
[1] J. Jacques, A. Collet, S. H. Wilen. Enantiomers, Racemates and Resolutions(1981).
[31] P. E. Kloeden, E. Platen. Numerical Solution of Stochastic Differential Equations(1992).
Set citation alerts for the article
Please enter your email address