Chiral forces in longitudinally invariant dielectric photonic waveguides
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.
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
Fig. 1. (a) Schematic showing the combined action of the translation of the cloud of enantiomers due to the optical force ( ) and the increase in the cloud size due to Brownian motion ( ). After a sorting time ( ), the clouds are separated by a distance of . Schematic adapted from Ref. [28]. (b) Sorting time of the enantiomer clouds as a function of the modulus of the separating optical force. Particles of different radii display different sorting time curves, from 1 to 1000 nm.
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].
Fig. 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 –direction that a quasi-TE mode in a SiN strip waveguide exerts on a particle, depending on its size. This is calculated for three different SiN waveguides, each one operating at a different wavelength: 405 nm, 780 nm, or 1310 nm. The SiN cross-section sizes ( ) of the waveguides are: (for ), (for ), and (for ). (d) Net transversal chiral and achiral forces that a quasi-CP mode in a strip waveguide exerts on a particle depending on the particle size. (e) Cross-section of the strip waveguide showing the position where the forces are evaluated for (c) and (d): at a vertical distance over the top of the waveguide equal to the radius of the particle, and at a horizontal distance equal to the fourth of the waveguide width from the center of the waveguide. The calculated forces are expected to be more accurate for particles whose size is smaller than the wavelength. This upper bound is , 290 nm, and 490 nm for , 780 nm, and 1310 nm, respectively.
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.
5.2 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 .
Fig. 3. Field properties ( , , and ) that generate the dominant forces of the quasi-TE mode in a strip waveguide ( ) at . The force stemming from the as well as the total optical force is represented for both chiralities of the particle ( ). All forces are calculated for a particle of 80 nm radius. The axes of all graphs refer to the – and –coordinates measured in μm units. The arrow map represents the transversal components of the vectorial quantities being plotted, and the colormap represents the scalar quantity or the –component of the vectorial quantity being plotted. The particle tracking graph shows the last position of 500 particles per enantiomer inside the microchannel ( ) after 2 s of motion given the total optical force field produced by the waveguide. The initial –coordinate of each particle was uniformly randomized between and 0.5 μm, and the initial –coordinate was uniformly randomized with any value within the microchannel. The orange circle represents the particle size to scale, whereas the magenta and green dots represent the center of mass of the enantiomeric particles.
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.
Fig. 4. Total optical force (transversal in arrow map, and longitudinal in colormap) that a quasi-TE mode in a wide strip waveguide ( ) at exerts on an 80 nm radius particle is shown for both chiralities of particle . The particle tracking graph shows the last position of 500 particles per enantiomer inside the microchannel ( ) after 2 s of motion given the total optical force field shown in the same figure. The initial –coordinate of each particle was uniformly randomized between and 0.5 μm, and the initial –coordinate was uniformly randomized with any value within the microchannel. The orange circle represents the particle size to scale, whereas the magenta and green dots represent the center of mass of the enantiomeric particles.
5.3 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 ).
Fig. 5. Field properties ( and ) that generate the dominant forces of the quasi-CP mode in a strip waveguide ( ). The force stemming from the gradient of as well as the total force is represented for both chiralities of the particle ( ). All forces are calculated for a particle of 52 nm radius. The arrow map represents the transversal components of the vectorial quantities being plotted, and the colormap represents the scalar quantity or the –component of the vectorial quantity being plotted. The particle tracking graph shows the last position of 500 particles per enantiomer inside the microchannel ( ) after 2 s of motion given the total optical force field produced by the waveguide. The initial –coordinate of each particle was uniformly randomized between and 0.5 μm, and the initial –coordinate was uniformly randomized with any value within the microchannel. The orange circle represents the particle size to scale, whereas the magenta and green dots represent the center of mass of the enantiomeric particles.
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.
Fig. 6. Net chiral and achiral forces exerted on a particle depending on its chirality parameter for (a) quasi-TE mode ( ) along the –direction and particle size and (b) quasi-CP mode ( ) and particle size of . The limit where the achiral force becomes stronger than the chiral force is at for the quasi-TE mode and for the quasi-CP mode. The forces were evaluated in the same position as in Figs. 2 (c) and 2 (d).
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.
[1]
[14]
[19]
Article Outline
Josep Martínez-Romeu, Iago Diez, Sebastian Golat, Francisco J. Rodríguez-Fortuño, Alejandro Martínez. Chiral forces in longitudinally invariant dielectric photonic waveguides[J]. Photonics Research, 2024, 12(3): 431.