Stokes trap for multiplexed particle manipulation and assembly using fluidics
- aDepartment of Mechanical Science and Engineering, University of Illinois at Urbana–Champaign, Urbana, IL 61801;
- bDepartment of Chemical and Biomolecular Engineering, University of Illinois at Urbana–Champaign, Urbana, IL 61801;
- cCenter for Biophysics and Computational Biology, University of Illinois at Urbana–Champaign, Urbana, IL 61801
See allHide authors and affiliations
Edited by David A. Weitz, Harvard University, Cambridge, MA, and approved March 1, 2016 (received for review December 20, 2015)

Significance
In recent years, techniques for trapping and manipulating small particles and single molecules have catalyzed major advances in several fields of science. Here, we describe a new method for confining, manipulating, and assembling multiple particles in free solution using the sole action of fluid flow. The Stokes trap provides a fundamentally new method for particle manipulation and assembly using a simple microfluidic device and nonlinear feedback control in the absence of optical, electric, or magnetic fields. This method holds strong potential to enable fundamental studies in soft materials, colloidal science, and biophysics.
Abstract
The ability to confine and manipulate single particles and molecules has revolutionized several fields of science. Hydrodynamic trapping offers an attractive method for particle manipulation in free solution without the need for optical, electric, acoustic, or magnetic fields. Here, we develop and demonstrate the Stokes trap, which is a new method for trapping multiple particles using only fluid flow. We demonstrate simultaneous manipulation of two particles in a simple microfluidic device using model predictive control. We further show that this approach can be used for fluidic-directed assembly of multiple particles in solution. Overall, this technique opens new vistas for fundamental studies of particle–particle interactions and provides a new method for the directed assembly of colloidal particles.
In recent years, particle-trapping methods have enabled a multitude of fundamental studies in the natural and applied sciences. To this end, the precise manipulation of small particles and single molecules has opened new windows into viewing biological and physical processes with unprecedented levels of control. Nearly all trapping methods rely on a fundamental set of forces for particle confinement, with techniques including optical traps (1⇓–3), magnetic traps (4), acoustic traps (5, 6), electrokinetic traps (7, 8) and the hydrodynamic trap (9⇓–11). Using these methods, researchers have studied phenomena ranging from the cooling of atoms to low temperatures (12) to the mechanistic behavior of single biomolecules (13).
Particle-trapping methods can be characterized by the nature of the underlying forces and the action of confinement. Trapping techniques can be broadly classified as passive (non–control-based) or active (control-based) traps. Passive traps allow for particle confinement over long periods of time without the need for active feedback control (e.g., optical traps), wherein a stable potential well is used to generate a restoring force to confine a particle at or near a target location. Active traps use a feedback-control algorithm for particle trapping, wherein disturbances such as thermal fluctuations can be corrected using active feedback control (e.g., electrokinetic or hydrodynamic traps). Despite the pervasive use of particle trapping techniques for fundamental investigations, several methods require complex experimental setups and impose specific restrictions on the properties of the trapped particle and the surrounding medium (e.g., dielectric constants or ionic strength). On the other hand, hydrodynamic traps impose no specific constraints on a particle’s surface chemistry, morphology, or refractive index (14). From this perspective, hydrodynamic methods possess several advantages for particle trapping and manipulation, and yet this approach has not yet been extended to multiple particles.
In 1934, G. I. Taylor developed the first active trap based on hydrodynamic flow (15). The four-roll mill was able to confine single macroscopic particles or immiscible droplets by continuous rotation of four cylinders, thereby generating a stagnation point flow. The original four-roll mill relied on manual (human) feedback control for particle confinement, which limits precision and controllability. In 1986, Bentley and Leal developed an automated four-roll mill by controlling the rotational speeds of the rollers using a computer (16). The automated four-roll mill, however, was based on a “macrofluidic” experimental flow cell, which complicates the confinement of small particles because of finite system response times. In recent years, the advent of microfluidics has allowed for fabrication of flow devices with small feature sizes, thereby enabling dynamic fluidic control with correspondingly small timescales. In 2003, a microfluidic cross-slot device was used to trap single DNA molecules near the stagnation point of a planar extensional flow, albeit using manual (tedious) human control over hydrostatic pressure in the outlet flow channels (17). A microfluidic four-roll mill was also developed and shown to generate extensional, rotational, and linear mixed flows (18, 19); however, these devices lack feedback control, which precludes long-term confinement of single particles and molecules.
In 2010, Schroeder and coworkers built and demonstrated an automated hydrodynamic trap capable of confining single particles at a target position for several minutes in a polydimethylsiloxane (PDMS)-based microfluidic device (9). In this device, an on-chip membrane valve was used to modulate flow rate in one outlet channel, thereby enabling particle trapping near a stagnation point for long times (10). By incorporating multiple on-chip valves, it was shown that single particles can be precisely manipulated in two dimensions using the sole action of fluid flow (11). In all cases, the automated hydrodynamic trap relies on a simple proportional-integral-derivative (PID) controller for particle manipulation (10, 20), which was sufficient for precise confinement of small particles in solution (e.g., 500-nm particles confined to within 180 nm of a set point position) (11). Despite these advances, however, the original automated hydrodynamic trap is strictly limited to trapping and manipulation of single particles. From this perspective, there is a strong need for new fluidic-based methods for precise confinement of multiple particles in solution. The ability to manipulate multiple particles in solution using simple microfluidics will enable a wide array of fundamental studies of particle–particle interactions and particle assembly. In particular, flow-based control over multiple particles will allow for detailed studies of soft particle interactions, including vesicle fusion, drop coalescence, and multiphase fluid interactions in small drops or particles.
In this work, we describe the Stokes trap, which is a new method for the manipulation of multiple particles using the sole action of fluid flow. The Stokes trap relies on precise control over particle position in viscous-dominated (low-Reynolds number) flows. Although prior work considered an open-loop algorithm for microfluidic assembly (21), this method was solely focused on computational studies of non-Brownian particles. Here, we develop and experimentally demonstrate a highly robust and scalable closed-loop control strategy that enables the precise manipulation of multiple particles and fluidic directed assembly of multiple particles.
Modeling and Design of the Stokes Trap
Fluidic Model and Governing Equations.
Consider the problem of manipulating P particles independently in a microfluidic device in which N channels intersect to form an N-sided polygonal flow device (Fig. 1A shows cross-slots with
Overview of the Stokes trap and control algorithm. (A) Schematics of four-channel and six-channel microfluidic devices for manipulating one and two particles, respectively. (B) Overview of the experimental setup. Inlet/outlet channels in the microfluidic device are connected to fluidic reservoirs that are pressurized by regulators controlled by a custom LabVIEW program. In this way, the fluidic reservoirs drive fluid flow in the microdevice. (C) Block diagram of the control loop for particle manipulation and trapping, with typical time constants of
We first consider the fluid dynamics within the microfluidic device. At low Reynolds number (Stokes flow conditions), we can neglect fluid inertia which yields the following conservation equations (22):
Using this model, the fluid velocity at a point inside the polygonal device is completely determined by the imposed flow rates
Model Predictive Control.
Given the governing equations for fluid flow, the next challenge lies in implementing a control algorithm to determine experimentally feasible flow rates
Steering a single particle from an initial position
We briefly discuss the MPC strategy followed by experimental implementation. First, a finite time horizon
For the general case of trapping P particles, we redefine the state vector
By solving the optimization problem, we obtain flow rates Q that move a particle from its current position to a target position at a minimal cost, as determined by the weighting factors. However, in our approach, we only apply the first set of flow rates
Results
Stokes Trap: Overview and Design.
A schematic illustration of the four- and six-channel microdevices is shown in Fig. 1A. In these experiments, we use a four-channel device for trapping a single particle and a six-channel device for trapping two particles. For trapping two particles, we found that an additional control variable
Fluidic Model and Streamline Topologies.
Before embarking on particle trapping, we first validated the fluidic model (Eq. 2) by experimentally determining streamline topologies using particle tracking and fluorescence microscopy. In these experiments, we use an aqueous glycerol buffer (viscosity, η = 12.6 cP) containing 2.2-μm-diameter fluorescent beads to characterize the flow field. We use standard bead tracking algorithms [ImageJ with MOSAIC plugin (25)] to determine bead trajectories over at least 30-s durations, which is suitable for quantitative flow field analysis at typical strain rates of 1–10
Streamline topologies in a six-channel microdevice from experiments and computation. (A and B) Schematic of the relative magnitude and direction of the flow rates for generating the streamline topologies in the figures below. Arrows pointing inwards represent flow entering the device, and arrows pointing outwards represent flow exiting the device. The size of the arrows signifies the relative magnitude of the flow rates. (C) Experimental streamlines showing the linked-arms topology, generated when the flow rates have a specific symmetry. Two stagnation points are clearly visible. (D) Experimental streamlines showing the non–linked-arms topology, generated if the symmetry in C is broken. (E and F) Streamline topologies obtained from numerical solution of Eq. 2. For display, streamlines emanating from inlet channels are plotted using distinct colors.
Trapping Single Particles Using the Stokes Trap.
We first characterized the performance of the Stokes trap by confining single particles at a target position and determining the trap stiffness. The Stokes trap can be used to confine particles in either quiescent conditions (i.e., no net flow if the particle is at the target position) or in net-flow conditions (i.e., net flow even if the particle is at target position). In quiescent conditions, fluid flow is only applied to correct for the Brownian motion of the particle, and in this case, the Stokes trap functions in an analogous fashion to electrokinetic traps (7). As an aside, we note that unlike the Stokes trap, the Generation I automated hydrodynamic trap cannot operate under quiescent conditions. The Generation I trap uses an on-chip membrane valve to control flow, and the membrane valve is incapable of influencing particle position in the absence of a net imposed flow.
To quantify trap performance, we used the Stokes trap with MPC (Eq. 5) to confine a single (2.2-μm-diameter) fluorescent bead in a four-channel device, and we recorded the particle trajectory for a duration of 400 s (Fig. 3A). We use a bead-tracking algorithm to localize particle position, followed by determination of the power spectral density (PSD) of the particle position fluctuations (Fig. 3B). We fit the PSD with a Lorentzian and a maximum-likelihood estimator (26, 27), and the corner frequency
Characterizing the performance of the Stokes trap. (A) Trajectory of trapped 2.2-μm-diameter bead over a period of 400 s. (Inset, Upper) Probability distribution of the position of the particle. (Inset, Lower) Schematic of microdevice used for trapping. (B) PSD particle position fluctuations for the trajectory shown in A. The PSD is analyzed to determine the corner frequency
Manipulating Two Particles Using the Stokes Trap.
We next used the Stokes trap to confine and manipulate two particles simultaneously (Fig. 4). Here, we used a six-channel microfluidic device to manipulate two 2.2-μm-diameter particles in an arbitrary scenario. First, the two particles are trapped ∼198 μm apart. Next, the target positions of both particles are instantaneously interchanged, and the MPC given by Eq. 5 is used to control the process of particle interchange. In this experiment, we do not prescribe the intermediate trajectory for either particle; rather, the trajectory is generated by the controller. We found that the controller successfully generates smooth trajectories in real time by minimizing the objective function, as shown in Fig. 4. Using this strategy, the particle-position interchange occurs over a relatively short timescale, with the entire process completing in 22.5 s (Movie S1).
Manipulating two particles using the Stokes trap, where the objective is to switch the center-of-mass positions of both particles. (A) At t = 0 s, two 2.2-μm-diameter fluorescent beads are initially trapped ∼198 μm apart. At t = 0.033 s, we instantaneously interchange the target positions of both particles, after which the controller generates a trajectory for each particle and calculates and applies the flow rates. (B and C) Tracing the two particle trajectories during the experiment, with the yellow line showing the past history of both particles. (D) The process finishes at t = 22.5 s, at which time the particle positions have been interchanged.
We further demonstrate the precise manipulation of two particles along preprogrammed paths (Fig. 5 and SI Appendix, Fig. S4). The experiment begins by tracking two freely suspended particles in the field of view, and at time t = 0 s, the controller is activated. First, the controller traps both particles, and the two particles are smoothly delivered to their respective (initial) target positions. Next, the target position is slowly and repeatedly stepped along a preprogrammed path, during which time the controller works to continuously move the particles synchronously with the target position. In this experiment, the target positions are programmed to trace the letter I, and the process completes after 328 s. As shown in Fig. 5, we are able to accurately and precisely control both particles along a predetermined path that spans hundreds of microns, which is orders of magnitude larger than the particle size (Movie S2).
Manipulating two particles using the Stokes trap, where the objective is to precisely control the paths of two 2.2-μm beads to trace the letter I. (A–D) Snapshots of both particles at various instants of time, with the yellow line showing the past history of both particles.
Fluidic-Directed Assembly of Particles.
We also used the Stokes trap to assemble two small particles in solution (Fig. 6). In this experiment, the goal is to link two “sticky” particles together via directed assembly, in particular by linking a 7.4-μm biotin-coated bead and a 10.6-μm streptavidin-coated bead by the strong noncovalent biotin–streptavidin interaction. At time t = 0 s, both particles are trapped at two separate positions. Next, the target position for the bead on the left (biotin-coated) is slowly moved toward the target position for the bead on the right (streptavidin-coated). As the beads approach each other, both particles begin to deviate from their target positions, which is a consequence of interparticle hydrodynamic interactions (22, 29). Nevertheless, particles follow their target positions at all interparticle distances using the MPC algorithm. For assembly experiments, the controller weight parameter γ in the objective function (Eq. 5) is increased slightly above the values used in two particle manipulation experiments described above, which amounts to a stricter penalization for not reaching the final target position. At t = 95 s, the two beads are firmly linked together by biotin–streptavidin interaction. The interparticle distance during the assembly event is plotted in Fig. 6, with snapshots of both particles shown on the right. For reference, the expected interparticle distance after linking is also plotted on the same figure. After linking the two particles together, we moved the assembled structure around in solution to confirm that the beads were tightly linked (Movie S3).
Fluidic-directed assembly of two particles using the Stokes trap. Interparticle distance is shown as a function of time during the assembly event between a 10.6-μm-diameter streptavidin-coated bead and a 7.4-μm-diameter biotin-coated bead. (Right) Snapshots of both particles as a function of time, where points a to d correspond to the interparticle distance plot. The beads are successfully linked around 95 s.
Discussion
In this paper, we demonstrate a new technique for manipulating and assembling microscale objects in free solution using the sole action of fluid flow. We experimentally implement model predictive flow control, which greatly simplifies the microfluidic device design for trapping by eliminating the need for integrated electrodes and on-chip membrane valves. The MPC algorithm is versatile and allows for the ability to control several particles independently, and as a proof of concept, we directly show the precise manipulation of two particles along arbitrary paths. We further use this technique to bring two particles in close contact, thereby facilitating linkage and directed assembly into a simple structure. Although simultaneous manipulation of multiple particles has been achieved using optical tweezers (2) and electrokinetic traps (8), multiparticle manipulation and assembly mediated only by fluid flow has not yet been achieved.
A major advantage of the MPC algorithm is the ability to correct for model imperfections and system perturbations in a robust manner. For example, our simple fluidic model for fluid velocity (and particle velocity) is based on a linear superposition of 2D point source flows around the perimeter of a polygonal flow cell. This model neglects interparticle hydrodynamic interactions, which give rise to perturbations in the fluid velocity that become significant upon close approach of the particles (e.g., when interparticle distance is less than ∼5 particle diameters). Nevertheless, as shown in the particle assembly experiment, the MPC algorithm is robust and successfully corrects limitations of the fluidic model.
A second advantage of implementing the MPC algorithm with distinct fluidic models is the ability to directly measure deviations of the system behavior from the model predictions. During an experiment, we record the transient and steady-state control input that effectively cancels any perturbations to the flow field. Interestingly, this approach allows for direct measurement of the fluid velocity arising from the flow field disturbances around trapped particles. In this manner, it should be possible to use one particle as a probe particle to record the fluid flow in the vicinity of a second particle. Thus, we believe that this technique can serve as a powerful tool to study interparticle interactions.
We anticipate that the Stokes trap can be used to study fundamental interactions between hard or soft particles, including vesicle fusion, polymersome fusion, droplet coalescence, coacervate droplet interactions and disassembly, red blood cell collisions, and interactions between living cells. Indeed, the Stokes trap places no restrictions on the physical properties of trapped particles (e.g., no restrictions on index of refraction or magnetic susceptibility), aside from the need to image particles. This versatility greatly extends the range of applicability of the method to a wide variety of studies, including manipulation of hard colloids, anisotropic particles, or soft particles. For example, our technique could serve as an alternative to micropipette aspiration in the study of cell adhesion, wherein one cell is typically confined using a small pipette before probing cell–cell interactions, a process that may perturb cell properties via physical confinement. Furthermore, fluidic trapping generally involves application of gentle, viscous-dominated flows that result in nonperturbative forces, thereby enabling trapping of particles ranging from nanometer-sized inorganic particles to giant vesicles, all with using a simple and robust modular setup. Moreover, fluidic confinement forces scale favorably and linearly with particle diameter R compared with
We further use the Stokes trap for the assembly of freely suspended particles in solution, that is, for particles with no external force or torque couples (e.g., no external electric fields or optical fields). In recent years, self-assembly of particles and molecules has been used to generate advanced functional materials with hierarchical structures (31). The process of self-assembly, however, mainly relies on self-organization at thermal equilibrium (32). From this perspective, it might be possible to access new levels of spatial ordering and hierarchical assembly by coupling external fields during the assembly process in a nonperturbative fashion. In the past, directed assembly of colloidal particles has been achieved by using external electrical fields (33, 34), acoustic fields (35), or optical fields (36). Nevertheless, not all techniques are well suited for directed assembly. Particle assembly with electrokinetic traps has been shown to be difficult, where the electrode potentials (voltages) required for bringing two 5-μm particles in contact diverged upon close approach, and as a consequence, these particles could not be brought closer than ∼8 μm (8). Moreover, during any assembly event, the applied force fields should be nonperturbative such that they do not alter the underlying system properties.
Despite its advantages, the Stokes trap has a few disadvantages for particle trapping. In the current implementation, we only control motion in the x-y plane, with no control over the position in the z direction. In addition, trapping occurs at an unstable equilibrium position in flow, which necessitates trapping via active feedback control, rather than passive trapping as achieved in optical traps.
In future work, it may be possible to use the Stokes trap to systematically build higher-order and more complex assemblies of particles. In the current version of the trap, we use n = 6 inlet/outlet channels, though it has been predicted that this can be scaled up to n = 10–15 while retaining controllability with physically reasonable flow rates (21). In the short term, however, the current device with n = 6 control variables could be used to control the 2D center-of-mass position of two particles while simultaneously controlling the orientation of one of the two particles. By adding one additional channel (n = 7), it would be straightforward to control the 2D orientation and center-of-mass position of two anisotropic particles.
Materials and Methods
Microfluidic Device Fabrication.
Standard procedures for PDMS-based soft lithography were followed to fabricate four- and six-channel microdevices (SI Appendix, SI Text, 1.1. Device Design and Fabrication).
Experimental Setup and Sample Preparation.
The microfluidic device was mounted on the stage of an inverted microscope (Olympus IX71) equipped with a 10× objective lens and a CCD camera (AVT Stingray F-033B/Point Gray GS3). The inlet/outlet channels on the microdevice were connected to fluidic reservoirs (Elveflow) through FEP tubing (IDEX Health & Science), and the reservoirs were pressurized using electronic pressure regulators (Proportion Air). In all experiments, the buffer is a glycerol–water solution with a viscosity of 12.6 Pa-s at room temperature (25 °C). For trapping experiments, we used 2.2-μm-diameter fluorescent polystyrene beads (Nile Red; Spherotech). For assembly experiments, we used 10.6-μm and 7.4-μm-diameter polystyrene beads coated with streptavidin and biotin (Spherotech and Bangs Laboratories, respectively).
Scaling and Nondimensionalization.
The optimization problem (Eq. 5) is solved numerically in real-time during experiments, and it is important to scale all quantities to be
Acknowledgments
We thank Prof. Sascha Hilgenfeldt for insightful discussions. We thank Prof. Paul J. A. Kenis and his group for providing access to cleanroom facilities for microdevice fabrication. This work was supported by a FMC Educational Fund Fellowship (to A.S.), a Packard Fellowship from the David and Lucile Packard Foundation, and National Science Foundation (NSF) Faculty Early Career Development (CAREER) Award Chemical, Bioengineering, Environmental, and Transport Systems (CBET) 1254340 (to C.M.S.).
Footnotes
- ↵1To whom correspondence should be addressed. Email: cms{at}illinois.edu.
Author contributions: A.S. and C.M.S. designed research; A.S. performed research; C.V.R. contributed new reagents/analytic tools; A.S. analyzed data; and A.S., C.V.R., and C.M.S. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1525162113/-/DCSupplemental.
References
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Ding X, et al.
- ↵
- ↵.
- Probst R,
- Cummins Z,
- Ropp C,
- Waks E,
- Shapiro B
- ↵
- ↵
- ↵
- ↵
- ↵.
- Selvin PR,
- Ha T
- ↵
- ↵
- ↵.
- Bentley BJ,
- Leal LG
- ↵.
- Schroeder CM,
- Babcock HP,
- Shaqfeh ESG,
- Chu S
- ↵
- ↵
- ↵.
- Shenoy A,
- Tanyeri M,
- Schroeder CM
- ↵
- ↵.
- Batchelor GK
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Leal LG
- ↵
- ↵.
- Hecht S
- ↵.
- Whitesides GM,
- Boncheva M
- ↵
- ↵
- ↵
- ↵
Citation Manager Formats
Article Classifications
- Physical Sciences
- Engineering
- Biological Sciences
- Biophysics and Computational Biology