Probing the chirality of a single microsphere trapped by a focused vortex beam through their orbital period
Kainã Diniz Tanja Schoger Arthur L. Fonseca Rafael S. Dutra Diney S. Ether Jr Gert-Ludwig Ingold Felipe A. Pinheiro Nathan B. Viana Paulo A. Maia Neto


fancy
plain
plain



















 


iblabel[1]#1            
akefntext[1]
[0pt][r]thefnmark #1




1.125 



*

§ 
0pt4pt4pt
*

 §.§ 
0pt15pt1pt




[LO,RE]

    < g r a p h i c s >

[CO]

    < g r a p h i c s >

[CE]

    < g r a p h i c s >

[RO]1–LastPage  
[LE]  1–LastPage

 







 
 










[
  \begin@twocolumnfalse

    < g r a p h i c s >
0pt[0pt][0pt]
    < g r a p h i c s >



    < g r a p h i c s >





    < g r a p h i c s >
     
Probing the chirality of a single microsphere trapped by a focused vortex beam through their orbital period
 

     
 


     Kainã Diniz,^a,b^† Tanja Schoger,^c † Arthur L. Fonseca,^a,b, 
 Rafael S. Dutra,^d
 Diney S. Ether Jr,^a,b Gert-Ludwig Ingold,^c 
 Felipe A. Pinheiro,^a
 Nathan B. Viana,^a,b 
 and Paulo A. Maia Neto^∗ a,b 



    < g r a p h i c s >
     When microspheres are illuminated by tightly focused vortex beams, they can be trapped in a non-equilibrium steady state where they orbit around the optical axis. 
By using the Mie-Debye theory for optical tweezers, we demonstrate that the orbital period strongly depends on the particle's chirality index. 
Taking advantage of such sensitivity, we put forth a method to experimentally characterize with high precision the chiroptical response of individual optically trapped particles. The 
method allows for an enhanced precision at least one order of magnitude larger than that of similar existing enantioselective approaches. It is particularly suited to probe the chiroptical response of individual particles, for which light-chiral matter
interactions are typically weak. 




 \end@twocolumnfalse 


]





§ 






^† These authors contributed equally to this work.

^a Instituto de Física, Universidade Federal do Rio de Janeiro, Caixa Postal 68528, Rio de Janeiro,  Rio de Janeiro, 21941-972, Brazil; Email: dinizkaina@gmail.com

^b CENABIO - Centro Nacional de Biologia Estrutural e Bioimagem, Universidade Federal do Rio de Janeiro,
Rio de Janeiro, Rio de Janeiro, 21941-902, Brazil

^c Institut für Physik, Universität Augsburg, 86135 Augsburg, Germany

^d LISComp-IFRJ, Instituto Federal de Educação, Ciência e Tecnologia, Rua Sebastão de Lacerda, Paracambi, Rio de Janeiro, 26600-000, Brasil















§ INTRODUCTION



Chiral discrimination plays a crucial role in many areas of science such as Chemistry, Molecular Biology and Pharmaceutics (see e.g. Ref.  for review).
Over the years, various methods to separate molecules and particles based on their chiral properties were developed. 
There exist chemical processes to separate enantiomers from each other (see  e.g. Refs.  for reviews), which, however, have the disadvantage that they are usually developed for specific chiral particles and tend to be invasive. In addition, they usually probe only the average chiral response of an ensemble of chiral particles or molecules, rather than that of individual particles, for which such response is typically small. <cit.> To circumvent this limitation, plasmonic nanostructures have been used in enantioselective schemes due to their ability to enhance chiroptical properties based on localized surface plasmon resonance. <cit.>
Recently, all-optical chiral discrimination methods have received significant attention due to their potential as noninvasive alternatives<cit.>, and because they are particularly suited to characterize the chiral response of single, isolated chiral nanoparticles <cit.>.
These methods are possible because chiral particles respond differently to left- and right-circularly polarized light.<cit.> This has been exploited, for instance, in the context of optical tweezers, <cit.> with several methods being introduced in recent years to trap and characterize single chiral particles using tightly focused beams. <cit.>
The proposal in Refs.  is based on an optical torque which the particle experiences when displaced from its equilibrium position on the optical axis by an external force. 
Due to focusing, the spin angular momentum associated with polarization can be exchanged with the trapped particle as orbital angular momentum<cit.>, generating a torque that is sensitive to the chirality of the particle.

In addition to spin, light can also carry intrinsic orbital angular momentum, which is associated with the field's phase distribution in space.<cit.> Paraxial beams that carry this type of angular momentum are called vortex beams.
An important class of vortex beams are the Laguerre-Gaussian modes, usually denoted by LG_pℓ, where p is a positive integer which determines the number of radial nodes, and ℓ is an integer called topological charge. 
In addition to spin angular momentum associated with polarization, such modes carry an orbital angular momentum of ℓħ per photon related to their helix-shaped wavefront, with the sign of ℓ determining the direction of the twist of the helices.
Upon interaction with such paraxial fields, a chiral dipole cannot discriminate between different topological charges. <cit.> 
An experiment with tightly focused vortex beams showed also no response of chiral molecules on beams with different topological charges. <cit.>
However, more recent studies revealed that chiral materials indeed respond in a discriminatory way to the handedness and magnitude of light's orbital angular momentum because of quadrupole contributions.<cit.> 
If the field becomes strongly focused, the spin and orbital degrees of freedom become coupled <cit.>, and a chiral particle’s response will be different for different topological charges.<cit.>

In the context of optical trapping, focusing of vortex beams with ℓ≠ 0 leads to a ring-shaped focal spot. 
If a particle is small compared to the diameter of the ring, it can be trapped in a non-equilibrium steady state where it orbits around the optical axis.<cit.>
For brevity, we refer to this type of state as the ring-trapping regime in the following.
Li et al. <cit.> found that, for a particle confined to the focal plane, both the radius of the orbit and the optical torque that drives the particle depend on its chirality. Also, it has recently been shown that optical tweezers with vortex beams with ℓ≠ 0 exert an enhanced torque upon trapped objects, and that this effect can be used to characterize material properties of microspheres. <cit.>
Here we propose to use the period of particles in the ring regime as a probe for their chirality. Beyond the usual discussion about enantioselectivity, we present a proposal to quantify microsphere's chirality while estimating the resolution that could be achieved.
Additionally, by calculating the radius of the orbit and its location along the axis from the conditions of 
vanishing axial and radial force components, we provide a more realistic model when compared to the ones which consider the azimuthal force only in the focal plane. 
We also demonstrate that, in our scenario, analyzing the period yields a higher chiral resolution than doing so with just the orbital radius. This result is particularly suited for enantioselection of individual particles, where chiroptical response is typically small, and for this reason our proposal singles out with respect to other existing enantioselective methods for single chiral particles.



§ MIE-DEBYE THEORY FOR CHIRAL NANOSPHERES TRAPPED BY A VORTEX BEAM


To describe the response of a chiral particle to an electromagnetic field, we use the following set of constitutive equations<cit.>

    [ 𝐃; 𝐁 ] = 
      [  ϵ_0 ϵ  i κ/c; - iκ/c  μ_0 μ ][ 𝐄; 𝐇 ] ,

where ϵ and μ are the relative permittivity and permeability,
c = 1/√(ϵ_0 μ_0) is the vacuum speed of light, and κ is a pseudo-scalar known as the chirality parameter. 
Although these equations assume a homogeneous and isotropic response, particles whose chirality arises from their geometry can also be considered in terms of an effective chirality parameter.<cit.>
Notice that κ accounts for an electro-to-magnetic and magneto-to-electric coupling.

To describe the trapping of a chiral spherical particle of radius R by a tightly focused vortex beam, we have developed a version of the Mie-Debye theory for optical tweezers with Laguerre-Gaussian modes <cit.> that includes chiral scatterers <cit.>. 
The field before focusing is assumed to be a circularly polarized (σ = ±  1) Laguerre-Gaussian beam  LG_0ℓ with one intensity node and topological charge ℓ.
The angular spectrum representation of the electric field resulting from the focusing of such a beam by an objective is given by <cit.>

    𝐄^(σ, ℓ)(𝐫)    = -ikf E_0 e^-ikf/2π(√(2)f/w_0)^|ℓ|∫_0^2π dφ  e^iℓφ
         ×∫_0^θ_0 dθ sinθ√(cosθ)sin^|ℓ|(θ) 	
    	e^- (fsin(θ)/w_0)^2
       ×
    	e^i𝐤·𝐫ϵ̂_σ (θ, φ) .

The integral covers the direction of all wave vectors 𝐤 = 𝐤(k, θ, φ) within the medium of refractive index n_w surrounding the sphere, up to a maximal angle defined by sin(θ_0) = NA/n_w, where NA is the numerical aperture of the objective that performs the focusing. 
The wave number k= 2π n_w/λ_0 is defined in terms of the vacuum wavelength λ_0 of the beam. E_0 denotes the field amplitude, while f defines the focal length and w_0 the beam waist at the entrance of the objective. 
The polarization unit vector is given by ϵ̂_σ(θ, φ) = e^iσφ(θ̂
 + iσφ̂)/√(2) where θ̂ and φ̂ refer to the unit vectors in spherical coordinates.
We obtain the scattered field by applying Mie theory.
The optical force 𝐅 exerted by the total field can be calculated by integrating the time-averaged Maxwell's stress tensor over a closed surface around the spherical scatterer.
In the context of the Mie-Debye theory, rather than working directly with the force, it is convenient to define the dimensionless quantity 𝐐 called efficiency factor<cit.>

    𝐐 = 𝐅/(n_w/c)P ,

where P is the power on the sample. 
The efficiency factor quantifies the force exerted by the field upon the particle per unit power. 
Due to the axial symmetry of the vortex beam, it is convenient to express the optical force in cylindrical coordinates 𝐐 = Q_ρρ̂ + Q_ϕϕ̂ + Q_z ẑ. 
The component Q_z defines the axial force along the propagation direction of the beam, while Q_ρ and Q_ϕ are the transverse force components in the radial and azimuthal direction, respectively.
Furthermore, we also define the position (ρ, z, ϕ) of the sphere with respect to the focus in cylindrical coordinates.
The explicit force expressions for a trapped dielectric sphere can be found in Refs. . 
For a chiral sphere, the electric and magnetic Mie scattering coefficients a_j and b_j of multipole order j have to be replaced by

    a_j → a_j + iσ d_j ,     
    b_j → b_j - iσ c_j  .

The scattering coefficients for a size parameter x=kR are given by 

    a_j(x)    =  Δ^-1_j(x) [V_j^(x) A_j^(x) + V_j^(x) A_j^(x)]  , 
    
    b_j(x)    =  Δ^-1_j(x) [W_j^(x) B_j^(x) + W_j^(x) B_j^(x) ]  , 
    
    c_j(x)    = - d_j(x) =   iΔ^-1_j(x) [W_j^(x) A_j^(x) - W_j^(x) A_j^(x)]  ,

where we used the following auxiliary functions

    Δ_j(x) = W_j^ (x) V_j^ (x) + W_j^ (x) V_j^ (x)  ,


    W_j^ (x)    = Mψ_j(N_ x) ξ'_j(x) -       ξ_j(x)ψ'_j(N_ x) , 
    
    V_j^(x)    =   ψ_j(N_ x) ξ'_j(x) - M ξ_j(x)ψ'_j(N_ x)  , 
    
    A_j^(x)    = Mψ_j(N_ x) ψ'_j(x) -      ψ_j(x)ψ'_j(N_ x)  , 
    
    B_j^(x)    =   ψ_j(N_ x) ψ'_j(x) - Mψ_j(x)ψ'_j(N_ x)
 
with the relative impedance M= n_w√(μ/ϵ) and the relative refractive index N_ = (n ±κ)/n_w, n = √(ϵμ) for left (L) and right (R) polarized waves. 
Note that we adapted the notation from Ref. , where similar Mie coefficients were obtained, but for different constitutive equations than the ones presented in Eq. (<ref>).
The Mie coefficients for the scattered field are expressed in terms of the Riccati-Bessel functions ψ_j(z) = z j_j(z) and ξ_j(z) = z h_j^(1)(z), where j_j(z) and h_j^(1)(z) are the spherical Bessel and Hankel functions of the first kind, respectively.
Due to the reciprocity of chiral materials, the polarization-mixing coefficients fulfill c_j = -d_j. 
If the chirality parameter vanishes, the coefficients reduce to the usual Mie coefficients a_j = A_j/W_j, b_j = B_j/V_j and c_j =0 = d_j.



§ RESULTS AND DISCUSSION



Using the Mie-Debye theory for chiral particles outlined above, we examine the period of stably trapped objects.
A particle in a steady-state orbit around the optical axis is in equilibrium in the axial and radial directions, which means that the optical force components in those two directions must vanish, as illustrated in Fig. <ref>.
To find the coordinates of the circular orbit ρ_ and z_, we use the Mie-Debye theory to simultaneously solve the equations

    Q_z(ρ_, z_) = 0 ,


    Q_ρ(ρ_, z_) = 0 .

We also require that the derivatives ∂_ρ Q_ρ and ∂_z Q_z at
(ρ_,z_) are negative to ensure that the orbit is stable.



As the particle is typically immersed in some fluid, it experiences a drag force proportional to its speed. <cit.> 
The particle will perform a uniform circular motion whose speed v_ϕ is such that the drag force and the azimuthal component of the optical force cancel each other. Thus, using the definition (<ref>), we find the following relation for the orbiting speed

    v_ϕ = n_wP/c γ Q_ϕ(ρ_, z_)  ,

where γ is the Stokes drag coefficient, i.e., the proportionality constant between the particle speed and the drag force.
Together with the relation v_ϕ = ρ_ω between the velocity and angular velocity ω, the period T=2π/ω can be expressed as

    T = 2πρ_γ/(n_w P/c) Q_ϕ(ρ_, z_) .

We characterize the liquid by a viscosity η and account for the influence of the walls of the sample by applying the Faxén correction to the drag coefficient of a spherical particle 
<cit.>

    γ = 6πη R/1 - 9/16R/h 
        + 1/8(R/h)^3 
        - 45/256(R/h)^4
        - 1/16(R/h)^5  ,

where h is the distance of the sphere's center from the interface.

We analyze the period, as given by Eq. (<ref>), and its dependence on the chirality index of a sphere for vortex beams with various topological charges.
The beam is assumed to be left-circularly (σ =1) polarized. 
For all numerical results discussed below, we assume an objective with numerical aperture NA = 1.2 and back aperture radius R_ obj = 2.8   mm, values that are typical for commercially available objectives. 
To make fair comparisons between different modes, one must ensure that they have similar filling conditions at the objective entrance port. 
Thus, for each Laguerre-Gaussian mode, except when ℓ=0, we compute the beam waist such that the ratio between the radius of the ring of maximal intensity and the objective equals 0.8, as it is described in detail in Ref. . 
This implies that the beam waist is given by  w_0(ℓ) = 0.8 R_ obj√(2/|ℓ|) for ℓ≠ 0.
We note that this type of dynamic waist control can be performed with light modulation devices<cit.>, the same that can generate vortex beams. 
For the Gaussian mode, we set w_0(ℓ=0)=2.2 mm. 
The microsphere center of mass is always taken to be h= 2 above the coverslip. 
We consider a non-magnetic scatterer with a refractive index n = √(ϵ) = 1.57 for a vacuum wavelength λ_0 = 1064 nm, so as to emulate a polystyrene microsphere. <cit.> 
The refractive index of water is n_w=1.32. <cit.> 
For the calculation of the Stokes drag force, we use the viscosity of water at 20  ^∘C, which is η = 1.0016 m Pa · s.
The power at the sample is set to P=10 mW. 
As we will see later, the choice of P in our method is a matter of experimental convenience, and the theoretical resolution for κ measurements is not directly affected by it.



Fig. <ref>(a) depicts the period of chiral and non-chiral spheres as a function of their radius for Laguerre-Gaussian beams of topological charges ℓ = ±4  and ℓ = ±5. These values were chosen so that the period could be shown for a variety of sphere radii. If the topological charge is too small, larger particles will be trapped on the beam axis. <cit.>
For all cases, the period for spheres with chirality index κ = 0.01 (dashed lines) and κ = -0.01 (dotted lines) is shown, as well as the period for a non-chiral sphere (solid lines). 
The shaded area between the dotted and dashed curves accounts for the period of spheres with a chirality index in the interval -0.01< κ < 0.01.
Notice that the curves for topological charges with the same absolute value but different signs are not the same. 
This happens because we are considering a left-circularly polarized beam before focusing (σ = 1), thus breaking the symmetry between the ±ℓ cases even for achiral spheres. 
Indeed, when the topological charge is positive, the orbital angular momentum has the same sign as the spin angular momentum, while for negative topological charges, the sign is opposite. 

Independently of the topological charge, all the curves exhibit the same general behavior. 
For radii R smaller than about 0.35m, the period monotonically increases as the radius decreases, while for larger radii, it exhibits oscillations. 
This can be understood in terms of a decomposition of the optical force into a conservative and a non-conservative component. 
The conservative component, usually called the gradient force, pulls the particle towards the region of maximum intensity. 
On the other hand, the non-conservative component, usually called the scattering force, arises from radiation pressure and from the field's non-uniform helicity <cit.>.
Since the beam before focusing is circularly polarized, the azimuthal component Q_ϕ
does not depend on ϕ, by azimuthal symmetry. 
Thus, the line integral of the optical force along a closed circle around the optical axis is proportional to Q_ϕ, 
showing that this component is non-conservative.
When the particle radius is small compared to the wavelength of the light (R≪λ), 
its scattering can be well described in the Rayleigh limit. In this limit, the conservative component is proportional to the gradient of the electric energy density and dominates the non-conservative one, which explains the strong suppression of Q_ϕ and the resulting increase in the period as the radius decreases. 

On the other hand, the non-conservative contribution builds up as the radius increases and becomes comparable to the wavelength (R≈λ_0/n_ w) in the Mie scattering regime, giving rise to an azimuthal force component that drives the particle on its orbit. 
The oscillations shown in Fig. <ref>(a) are a consequence of interference effects inside the sphere.<cit.>

In spite of the overall similar behavior discussed above, Fig. <ref>(a) shows a clear split between the curves corresponding to κ=-0.01 (dotted) and κ=0.01 (dashed).
The rotation period decreases monotonically with the chirality index as is exemplified in Figure <ref>(b) for a sphere of radius R = 0.3. 
An approximately linear relationship exists between the chirality parameter and the period for a fixed topological charge.
Variations in the chirality index δκ_ℓ are thus directly proportional to variations of the period δ T_ℓ, i.e.

    δκ_ℓ = |b_ℓ| δ T_ℓ  ,

where b_ℓ is the slope of the linear fit of the κ(T)-curves as illustrated in Fig. <ref>(b).
Notice that we have added an index ℓ to the error in period measurements δ T_ℓ. Since the period monotonically increases with |ℓ| <cit.>, a fixed uncertainty would mean that the precision at higher topological charges is greater than that at smaller ones. Then, any gain in the resolution δκ_ℓ could be considered as an artifact of assuming a progressively smaller relative uncertainty. In order to allow for a fair comparison between different modes, we assume that the period is measured with the same  relative uncertainty ξ for all modes and define

    δ T_ℓ = ξT_ℓ  ,

where T_ℓ is the average period for the mode ℓ in the considered κ interval from -0.01 to 0.01.


Using the definition (<ref>), we have investigated quantitatively the κ-resolution that could be achieved through period measurements.
Figure <ref> displays δκ_ℓ as a function of the topological charge ℓ for beads of radii 0.15, 0.25 and 0.35 scaled by the relative error of the period ξ. 
In each case, we plot all the values of ℓ for which a well-defined on-ring position exists.
For values of ℓ below those displayed in the set of points corresponding to each radius, the respective particle would be trapped on the optical axis. 
On the other hand, for values of ℓ greater than those presented, no point in space satisfies Eqs. (<ref>) and (<ref>) simultaneously for the given radius values, and no optical trapping is possible.
Notice that this upper limit for the available ℓ values can only exist in the Mie scattering regime, where the scattering component of the optical force plays an important role.
Indeed, in the Rayleigh regime the gradient force will necessarily pull the particle towards the ring of maximum intensity.
This is in accordance with the fact that the number of available topological charges decreases as the particle becomes larger.



For R = 0.15, the precision in κ measurements monotonically increases with |ℓ|. 
The lowest δκ_ℓ value achieved was δκ_-14≈ 0.24 ξ. 
On the other hand, the results for R = 0.25 and R = 0.35 show that δκ_ℓ finds its minima for much lower values of the topological charge. 
This can be advantageous, since the period for these modes is much smaller at the same power, allowing for good statistics with less acquisition time. 
The lowest values of δκ_ℓ obtained for R = 0.25 and R = 0.35 were δκ_2≈ 0.16 ξ and δκ_1≈ 0.19 ξ, respectively. 
Hence, for a relative uncertainty ξ=10^-3 of the period measurement<cit.>, we would find a precision of the order of 10^-4 for the chirality measurement. 
Due to the linear relation (<ref>), improving the precision of period measurements would lead to a proportional enhancement in the chiral resolution of our method.

It is worth noting that a change in radius appears to displace the points globally, i. e., it either enhances or diminishes precision across all ℓ values for the examples shown in Fig. <ref>. 
Also, there is no monotonic relationship between δκ_ℓ and R. 
The general κ-resolution enhances from R = 0.15 to R = 0.25, but worsens from the latter to R = 0.35. 
This suggests the existence of an optimal radius for which the period is most sensitive to κ. 
In order to find such a value, we have performed a calculation of δκ_ℓ/ξ as a function of the particle radius for fixed topological charges ℓ = ± 4 and ℓ=±5, and the results are depicted in Fig. <ref>. 
The curves exhibit global minima near R ≈ 0.28 , which 
seems to be the radius allowing for the most sensitive chirality measurement. 
In the region R < 0.28, the resolution progressively worsens as the microsphere becomes smaller. 
On the other hand, in the region R > 0.28, δκ_ℓ exhibits an oscillatory behavior, meaning that the resolution of the measurement is of the same order of magnitude for particles within that region.



We would like to highlight an interesting aspect of Eq. (<ref>). 
Since the
period scales with power as T∼ 1/P, 
the slope b_ℓ = (∂κ/∂ T)_κ = 0 goes as b_ℓ∼ P,
and then δκ_ℓ/ξ, as defined by
Eqs. (<ref>) and (<ref>),
is power-independent, and so are the arguments developed throughout this section. 
Hence, in a real implementation, the power can be chosen such that the precision of the period measurement is maximized.
For example, when using large values of |ℓ|, one may freely increase the power in order to reduce the period and thus reduce the data acquisition time necessary to perform a good statistical analysis. 
In addition, increasing the power also reduces the effect of the particle's Brownian fluctuations, allowing for more precise determinations of periods, and thus providing a greater chiral resolution.


In Ref. , the authors show that, for particles confined to the focal plane, the average orbital radius depends on κ. 
Inspired by their work, we have also investigated the possibility of characterizing a particle's chirality through the orbital radius, rather than using the period.
In Fig. <ref>(a), the orbit coordinates ρ_ and z_ as well as the azimuthal force efficiency in the ring regime Q_ϕ = Q_ϕ(ρ_, z_) are shown as functions of κ, normalized by their value at κ = 0. In the represented interval, all quantities exhibit linear behavior, but it can be observed that the azimuthal force varies more rapidly than the orbital radius. This fact is not just a particularity of the chosen radius, as can be seen in Fig. <ref>(b), where we plot the same relative quantities as functions of the sphere radius for different chirality indices. From Eq. <ref>, we see that the period is proportional to ρ_ and inversely proportional to Q_ϕ, and therefore the dependence of the period on κ is mainly due to Q_ϕ. Thus, a measure of κ based solely on the orbital radius, even if done with the same precision as period measurements, would necessarily have lower chiral resolution than a measurement made through the period. The stronger variation of the azimuthal force with the chirality index also explains why the period shown in Fig. <ref> decreases with increasing κ.
Moreover, it should be noted that highly precise measurements are easier to perform for the period than for the radius of the orbit.
By measuring the back<cit.> or forward-scattered <cit.> light and Fourier transforming the signal to produce a power spectrum, one can extract the orbital frequency.
In contrast with the radius, the period depends on externally tunable parameters like the beam power, the Stokes drag coefficient, and the topological charge, thus allowing for optimizing the measurement.








§ CONCLUSIONS


In conclusion, we have introduced a method to measure the chirality index of micro-sized particles with a precision up to 10^-4. 
The method is based on measuring the period of a particle trapped by a focused vortex beam. 
The resolution does not depend on the power of the beam.
Compared to similar existing proposals <cit.> our method offers a gain of at least one order of magnitude in precision. This result is of particular interest for probing and characterizing the chiroptical response of single, individual particles that typically exhibit weak light-chiral matter interactions. The precision of this method can be further improved if period measurements with a relative error below 10^-3 can be achieved. Our findings may have applications in enantioselection of particles with very small chiral indexes, such as particles made of naturally occurring materials. 





§ CONFLICTS OF INTEREST

The authors declare no conflicts of interest.



§ ACKNOWLEDGEMENTS


We are grateful to Cyriaque Genet, Guilherme Moura, Leonardo Menezes, and Paula Monteiro for fruitful discussions. 
P. A. M. N.,  N. B. V., and F. A. P. acknowledge funding from the Brazilian agencies Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq–Brazil), Coordenação de Aperfeiçamento de Pessoal de Nível Superior (CAPES–Brazil),  Instituto Nacional de Ciência e Tecnologia de Fluidos Complexos  (INCT-FCx), and the Research Foundations of the States of Rio de Janeiro (FAPERJ) and São Paulo (FAPESP).




rsc 

