The importance of cell manipulation and cultivation is increasing rapidly in various fields, such as drug discovery, regenerative medicine, and investigation of new energy sources. This paper presents a method to transport cells in a microfluidic chip without contact. A local vortex was generated when high-frequency oscillation of a microtool was induced in a microfluidic chip. The vortex was controlled by tuning the tool's oscillation parameters, such as the oscillation amplitude and frequency. The cells were then transported in the chip based on the direction of the tool's movement, and their position, posture, and trajectories were controlled. Bovine oocyte manipulations, that is, transportation and rotation, were conducted to demonstrate the capability of the proposed method, without any contact by the microrobot with high-frequency oscillation.

In biology, the study of cells attracts considerable attention because the cell is the basic structural, functional, and biological unit of all known living organisms. The importance of cell manipulation and cultivation has been increasing rapidly in various fields, such as drug discovery, regenerative medicine,1 and investigation of new energy sources. Biomanipulations widely used in many applications show great commercial value, with over tens of billions of dollars invested into this field worldwide.2–6 There are two types of methods for cell manipulation: contact7–9 and noncontact manipulations. To achieve much more flexible actuation in a closed space, researchers have recently begun studying the noncontact manipulation of robots, especially for biological applications. Noncontact methods show great merits on non-contaminations and remote control. Therefore, many optional methods without contact, such as electric force,10,11 optical force,12,13 acoustic force,14,15 and magnetic force,16–18 have been well studied. For example, a rotating magnetic field was applied to drive Ni nanowires with a tumbling motion.19–21 The local flows around the nanowire can be used to manipulate micro-objects without contact. However, the applied force on the object was small, and thus insufficient for the transportation of some large objects, such as mammalian oocytes.

In this paper, we propose a method for noncontact cell transportation by using permanent magnets to oscillate a microtool at a high frequency in a microfluidic chip. A local vortex could be generated when the tool is oscillated at a high frequency, and the fluid stream transports cells without tool contact. Figure 1 shows the concept of noncontact manipulation through the generated vortex. The vortex form varies depending on the tool oscillation form, and can thus be controlled by changing several parameters, such as the oscillation amplitude, frequency, and relative position between the tool and microchannel wall. Owing to their attraction to the vortex, the cells transport through the microtool, which oscillates at a high frequency in a microfluidic chip. The tuning of the oscillation parameters of the tool changes the viscous fluid dynamics, and the cells can be manipulated away from the tool without contact. As this is a noncontact operation, there are no concerns regarding cell adhesion to the tool, unlike for other methods in which there are occasional issues when the manipulator handles adherent cells. In addition, the tool can manipulate cells through contact as well because it has enough physical strength to handle cells.

FIG. 1.

Concept of noncontact cell manipulation through high-frequency oscillation in a microfluidic chip.

FIG. 1.

Concept of noncontact cell manipulation through high-frequency oscillation in a microfluidic chip.

Close modal

A steady stream, which is a function of oscillation, is generated when an object is oscillated in the viscous fluid.22–24 However, the Reynolds number is generally small at a microscale, and a high frequency is required to generate the vortex. Assuming that the steady stream occurs two-dimensionally when a plate oscillates in a microfluidic chip, the local fluid velocity (Ux,Uy) can be calculated using the model proposed by Wang et al. and the modified Reynolds number Re = ωa2/μ, where ω is the frequency of the plate oscillation, a is the gap between the oscillating plate and channel wall, and μ is the kinematic viscosity. Then,

(1)
(2)
(3)

where fy,λ is the function that considers the oscillation effects and A is oscillation amplitude. The imaginary number i is defined solely based on the property that its square is −1, and Re stands for Reynolds number.25 

Figure 2 shows the simulation results when the microtool was oscillated using the permanent magnets with a ±0.5 mm amplitude in the Y direction only, with a sine wave oscillating frequency of 30 Hz. Under this condition, a regular vortex was generated, and the speed of the vortex reached to tens of micrometers per second. The vortex direction regularly changed with the movement of the oscillating microtool. This result implies that a tool with high-frequency oscillation in a microfluidic chip can transport cells without contact when the cells are attracted by the vortex, and the vortex can be controlled by configuring the oscillation.

FIG. 2.

Simulation results of velocity vector field in a microfluidic chip (arrow represents the velocity vector, and the contour is the magnitude). The microrobot oscillated in the (a) positive direction of the Y-axis and (b) negative direction of the Y-axis.

FIG. 2.

Simulation results of velocity vector field in a microfluidic chip (arrow represents the velocity vector, and the contour is the magnitude). The microrobot oscillated in the (a) positive direction of the Y-axis and (b) negative direction of the Y-axis.

Close modal

By using the previous calculation results, the drag force from the generated vortex acting on a single cell can be estimated by

(4)

where Fd is the drag force, ρ is the fluid density, v is the speed of the object relative to the fluid, Cd is the drag coefficient, and A is the reference area. From the equation, we can easily find that the drag force is relative to the fluid velocity. In our experimental condition, we used the culture medium, the kinematic viscosity of which approached that of water so that ρ=1g/cm3, and the frequency is 30 Hz. When the microrobot moved in the microchip, the average gap between the oscillating plate and channel wall was 35 mm. Therefore, the Reynolds number was approximately equal to 2×105. The drag coefficient can be approximated as follows:26 

(5)

The reference area is estimated as A=π×R23.14×75×1062=1.77×108m2 (R is the cross-sectional radius). Thus,

(6)

According to Eqs. (1), (2), and (6), the drag force is proportional to the oscillated frequency, thus when the frequency increases the force also increases.

The fabrication process of microrobot is described as follows.27 The silicon-based microrobots can be made into diverse shapes through deep reactive-ion etching (DRIE). First, the photoresist OFPR (Tokyo Ohka Kogyo, Co.) was coated on the substrate. After exposure on the OFPR side, an OFPR pattern was developed. Next, DRIE was conducted from the OFPR side, and the silicon was etched to a depth of 125 μm. By exposing both sides, we coated the opposite side of the silicon with OFPR and performed the patterning and exposure after the DRIE process. Then, we finished creating the silicon-based microrobot. The images obtained using the scanning electron microscope show the two-time DRIE-fabricated microrobot. Last, four neodymium (Nd2Fe14B) magnets (Magfine, Co., Japan), that is, φ0.5 mm × 0.5 mm permanent magnets, were assembled to shape the microrobot.

To control the microrobot oscillation in a closed microfluidic chip, the magnetic control had many obvious advantages, such as non-contact, high-speed actuation28 (up to 90 Hz), and precise positioning accuracy29 (micrometer order) in the XY plane. Therefore, a magnetically driven microtool18 (MMT) was the optimal option for the oscillating tool. The experiments showed the velocity vector field on the microfluidic chip when the microrobot was moving with high-frequency oscillation. Fluorescent beads of 2 μm diameter were added to the fluid to visualize the actual distribution of the vortex. The shape of the microrobot made of silicon and nickel is shown in Fig. 1, in which four permanent magnets are fixed at four corners. The size of the oscillation tool was 0.3 mm (width) × 0.2 mm (height), and the width and height of the microfluidic chip were 3 cm and 0.5 mm, respectively.

When the microrobot oscillated at the Y-axis by permanent magnets, objects were pulled to the edge of the microrobot. Then, only the object was transported by moving the microrobot on the XY plane. Therefore, by moving the microrobot, the fluorescent beads can easily be arranged in the plane and the microfluid distribution can be visible, as shown in Fig. 3 (Ref. 30) (Multimedia view). Figure 3(a) shows the high-frequency oscillating microrobot. Figure 3(b) shows the movement of the microrobot from left to right, whereas Fig. 3(c) shows the movement in the opposite direction. To move the microrobot from the rear side to the front side, the fluorescent beads can be arranged as shown in Fig. 3(d). Figure 3 also shows the measured velocity vector of the fluorescent beads under a set condition: amplitude of ±0.5 mm and frequency of 30 Hz. The object velocities were estimated using a method based on image processing. In this study, we used the block matching theory, which estimates motion between two images or two video frames using “blocks” of pixels. Although the velocity vector field was slightly different from the simulation results because of the limitations of the robotic shape, the vortex around the edge was based on the simulation analysis. According to velocity vector analysis, the highest velocity of the fluid was at the edge of the tool. Therefore, the cell transportation and rotation could be achieved by these vortexes. In addition, in the process of moving the microrobot with high-frequency oscillation, the vector field slightly changed, and the vortex changed with change in the oscillating conditions. Of course, these are only some simple experiments to prove that the microrobot with high-frequency oscillation can produce a certain flow field distribution. However, more flow fields can also be achieved by adjusting the fluid factors, such as oscillation amplitude and frequency.

FIG. 3.

Experimental result using fluorescent beads (Φ2μm): (a) microrobot oscillation; (b) microrobot moving from left to right; (c) microrobot moving from right to left; and (d) microrobot moving from rear to front (the red arrow represents the microrobot movement direction, and the blue arrow represents the oscillation direction. Multimedia view: https://doi.org/10.1063/1.5009545.1

FIG. 3.

Experimental result using fluorescent beads (Φ2μm): (a) microrobot oscillation; (b) microrobot moving from left to right; (c) microrobot moving from right to left; and (d) microrobot moving from rear to front (the red arrow represents the microrobot movement direction, and the blue arrow represents the oscillation direction. Multimedia view: https://doi.org/10.1063/1.5009545.1

Close modal

Figure 4 shows the results of the hydrodynamic force before and during oscillation. The oscillating conditions included amplitude of ±0.5 mm, oscillation frequency of 30 Hz. After the velocity of the fluid was analyzed, the hydrodynamic force distribution was calculated using the equation of drag force [Eq. (6)]. The simulation shows that the hydrodynamic force dragging a single cell could reach the order of micro-Newton. According to Fig. 4, before the oscillation, only a slight force was generated away from the tool because of the external interference. Moreover, when the microrobot was oscillating, the drag force around the tool edge was significantly higher than that at other areas, while the hydrodynamic force was the highest at the tip area of the microrobot. Thus, this hydrodynamic force can drag cells when the microrobot with high-frequency oscillation moves.

FIG. 4.

Hydrodynamic force distribution in a chip: (a) before the microrobot oscillated and (b) during microrobot oscillation.

FIG. 4.

Hydrodynamic force distribution in a chip: (a) before the microrobot oscillated and (b) during microrobot oscillation.

Close modal

The basic experiments and simulations of the velocity vector field showed that with microrobot with high-frequency oscillation was able to transport cells without contact. Therefore, by using the bovine oocytes as the manipulating object, we demonstrated its actual transport state. Except for the diameter of the oocyte, which was 150 μm, the other conditions were the same as those in the experiments of fluorescent beads, that is, the oscillating amplitude was 0.5 mm and frequency was 30 Hz. Figure 5 shows that by using a microrobot, the cell can be easily transported to the designated spot (Multimedia view). Figure 5(a) shows that the cell was stationary and away from the microrobot before its oscillation. The oocyte can be transported from the front side to the rear side as shown in Fig. 5(b). Figure 5(c) shows the transportation of the oocyte from left to right, while Fig. 5(d) shows the transportation in the opposite direction. If the relative position of the cell and microrobot was very close, cells would be occasionally adsorbed into the tool. However, this had no significant effect on the process of transporting other cells without contact. If necessary, the microrobot can also perform cell manipulation, such as cutting and injection, through contact, as presented in our previous works.31,32

FIG. 5.

Bovine oocyte manipulation using the microrobot: (a) oocyte kept stationary; (b) oocyte was transported to the rear side; (c) oocyte was transported from left to right; and (d) oocyte was transported from right to left (the red arrow represents the microrobot movement direction, and the blue arrow represents the oocyte transporting direction. Multimedia view: https://doi.org/10.1063/1.5009545.2

FIG. 5.

Bovine oocyte manipulation using the microrobot: (a) oocyte kept stationary; (b) oocyte was transported to the rear side; (c) oocyte was transported from left to right; and (d) oocyte was transported from right to left (the red arrow represents the microrobot movement direction, and the blue arrow represents the oocyte transporting direction. Multimedia view: https://doi.org/10.1063/1.5009545.2

Close modal

In summary, the movement of the microrobot with high-frequency oscillation in the microfluidic chip generated a velocity vector field, the distribution of which depends on the conditions of the oscillation and movement trajectory. Cells could be attracted by the vortex flowing around the edge of the oscillating microrobot and follow the tool without contact when the microrobot moved in the microfluid chip. This noncontact manipulation has many advantages, including larger control of the area and the movement of the cell in a predetermined trajectory. Although sometimes cells accidentally adhered to the tool, other noncontact manipulations were not disturbed, and the tool could conduct contact operation on the cell.31,32

In this paper, we demonstrated the noncontact transportation of a single cell. It is also possible to simultaneously manipulate several cells by using the experimental results of fluorescent beads. In addition, the tuned movement trajectory or oscillation parameters of the microrobot can change the transporting route of the cell, thus achieving the manipulation requirements. Our future work would involve the achievement of high-precision movement of noncontact manipulations in a three-dimensional microfluidic chip.

This work was supported in part by the National Natural Science Foundation of China (61602022) and ZG216S1751, ZG216S1745.

1.
V. S.
Nirmalanandhan
and
G. S.
Sittampalam
,
J. Biomol. Screening
14
,
755
(
2009
).
2.
I.
Economidis
, “An overview of the biotechnology research activities in the European Community,” in
Modern Biochemical Engineering
, Advances in Biochemical Engineering/Biotechnology, vol 46 (
Springer
,
Berlin, Heidelberg
,
1992
), pp.
225
234
.
3.
C.
Juma
and
V.
Konde
,
Environment
44
,
22
(
2002
).
4.
S.
Liyanage
and
P.
Gluckman
,
Int. J. Biotechnol.
6
,
281
(
2004
).
7.
W.
Hu
,
K. S.
Ishii
, and
A. T.
Ohta
,
Appl. Phys. Lett.
99
,
094103
(
2011
).
8.
Y.
Sun
,
J. S.
Evans
,
T.
Lee
,
B.
Senyuk
,
P.
Keller
,
S.
He
, and
I. I.
Smalyukh
,
Appl. Phys. Lett.
100
,
241901
(
2012
).
9.
D. H.
Kim
,
U. K.
Cheang
,
L.
Kohidai
,
D.
Byun
, and
M. J.
Kim
,
Appl. Phys. Lett.
97
,
173702
(
2010
).
10.
M.
Barbic
,
J. J.
Mock
,
A. P.
Gray
, and
S.
Schultz
,
Appl. Phys. Lett.
79
,
1399
(
2001
).
11.
S.
Grilli
and
P.
Ferraro
,
Appl. Phys. Lett.
92
,
232902
(
2008
).
12.
S.
Sato
,
M.
Ishigure
, and
H.
Inaba
,
Electron. Lett.
27
,
1831
(
1991
).
13.
M.
Ichikawa
,
K.
Kubo
,
K.
Yoshikawa
, and
Y.
Kimura
,
J. Biomed. Opt.
13
,
10503
(
2008
).
14.
M.
Cecchini
,
S.
Girardo
,
D.
Pisignano
,
R.
Cingolani
, and
F.
Beltram
,
Appl. Phys. Lett.
92
,
104103
(
2008
).
15.
L.
Meng
,
F.
Cai
,
J.
Chen
,
L.
Niu
,
Y.
Li
,
J.
Wu
, and
H.
Zheng
,
Appl. Phys. Lett.
100
,
173701
(
2012
).
16.
M. S.
Sakar
,
E. B.
Steager
,
D. H.
Kim
,
M. J.
Kim
,
G. J.
Pappas
, and
V.
Kumar
,
Appl. Phys. Lett.
96
,
43705
(
2010
).
17.
L.
Zhang
,
J. J.
Abbott
,
L.
Dong
,
B. E.
Kratochvil
,
D.
Bell
, and
B. J.
Nelson
,
Appl. Phys. Lett.
94
,
64107
(
2009
).
18.
M.
Hagiwara
,
T.
Kawahara
,
Y.
Yamanishi
, and
F.
Arai
,
Appl. Phys. Lett.
97
,
013701
(
2010
).
19.
K. E.
Peyer
,
L.
Zhang
, and
B. J.
Nelson
,
Appl. Phys. Lett.
99
,
174101
(
2011
).
20.
L.
Zhang
,
T.
Petit
,
Y.
Lu
,
B. E.
Kratochvil
,
K. E.
Peyer
,
R.
Pei
,
J.
Lou
, and
B. J.
Nelson
,
ACS Nano
4
(
10
),
6228
6234
(
2010
).
21.
L.
Zhang
,
T.
Petit
,
K. E.
Peyer
, and
B. J.
Nelson
,
Nanomed.: Nanotechnol., Biol., Med.
8
,
1074
1080
(
2012
).
22.
C. Y.
Wang
,
J. Fluid Mech.
32
,
55
68
(
1968
).
23.
T.
Sarpkaya
,
J. Fluid Mech.
165
,
61
71
(
1986
).
24.
M.
Tatsuno
and
P. W.
Bearman
,
J. Fluid Mech.
211
,
157
182
(
1990
).
25.
C. Y.
Wang
and
B.
Drachman
,
Appl. Sci. Res.
39
,
55
68
(
1982
).
26.
J. P.
Owen
and
W. S.
Ryu
,
Eur. J. Phys.
26
,
1085
1091
(
2005
).
27.
L.
Feng
,
P.
Di
, and
F.
Arai
,
Int. J. Rob. Res.
35
,
1445
1458
(
2016
).
28.
M.
Hagiwara
,
T.
Kawahara
,
Y.
Yamanishi
,
T.
Masuda
,
L.
Feng
, and
F.
Arai
,
Lab Chip
11
,
2049
2054
(
2011
).
29.
M.
Hagiwara
,
T.
Kawahara
,
T.
Iijima
,
Y.
Yamanishi
, and
F.
Arai
, in
Proceedings of the IEEE International Conference on Robotics and Automation (ICRA)
,
St. Paul Minnesota, USA
,
2012
, pp.
2517
2522
.
30.
R. J.
Dijkink
,
J. P.
Dennen
,
C. D.
Ohl
, and
A.
Prosperetti
,
J. Micromech. Microeng.
16
,
1653
1659
(
2006
).
31.
L.
Feng
,
M.
Hagiwara
,
A.
Ichikawa
, and
F.
Arai
,
Micromachines
4
,
272
285
(
2013
).
32.
L.
Feng
,
Y.
Sun
,
C.
Ohsumi
, and
F.
Arai
,
Biomicrofluidics
7
,
054113
(
2013
).