Trap profiles of projector based optoelectronic tweezers ( OET ) with HeLa cells

In this paper we present trap profile measurements for HeLa cells in Optoelectronic Tweezers (OET) based on a data projector. The data projector is used as a light source to illuminate amorphous Si creating virtual electrodes which are used to trap particles through dielectrophoresis. We show that although the trap stiffness is typically greater at the edges of the optical spot it is possible to create a trap with constant trap stiffness by reducing the trap’s size until it is similar to the object being trapped. We have successfully created a trap for HeLa cells with a constant trap stiffness of 3x10 -6 Nm -1 (capable of moving the cell up to 50μms -1 ) with a 12μm diameter trap. We also calculate the depth of the potential well that the cell will experience due to the trap and find that it to be 1.6x10 -16 J (4x10 4 kBT). 2009 Optical Society of America OCIS codes: (350.4855) Optical tweezers or optical manipulation References and links 1. P. Y. Chiou, A. T. Ohta, and M. C. Wu, "Massively parallel manipulation of single cells and microparticles using optical images," Nature 436, 370-372 (2005). 2. S. L. Neale, M. Mazilu, J. I. B. Wilson, K. Dholakia, and T. F. Krauss, "The resolution of optical traps created by light induced dielectrophoresis (LIDEP)," Opt. Express 15, 12619-12626 (2007). 3. P. Y. Chiou, A. T. Ohta, and M. C. Wu, “Toward all optical lab-on-a-chip system: optical manipulation of both microfluid and microscopic particles,” Proceedings of the SPIE, Volume 5514, 73-81 (2004). 4. Y.S. Lu, Y.P. Huang, J.A. Yeh, C. Lee and Y.H. Chang, “Controllability of non-contact manipulation by image dielectrophoresis” Opt. Quantum Electron. 37, 1385-1395 (2005). 5. H. Hwang, Y. J. Choi, W. Choi, S. H. Kim, J. Jang, and J. K. Park, "Interactive manipulation of blood cells using a lens-integrated liquid crystal display based optoelectronic tweezers system," Electrophoresis 29, 1203-1212 (2008). 6. A. T. Ohta, P. Y. Chiou, H. L. Phan, S. W. Sherwood, J. M. Yang, A. N. K. Lau, H. Y. Hsu, A. Jamshidi, and M. C. Wu, "Optically controlled cell discrimination and trapping using optoelectronic tweezers," IEEE J. Sel. Top. Quantum Electron. 13, 235-243 (2007). 7. A. T. Ohta, P. Y. Chiou, T. H. Han, J. C. Liao, U. Bhardwaj, E. R. B. McCabe, F. Yu, R. Sun, and M. C. Wu, "Dynamic Cell and Microparticle Control via Optoelectronic Tweezers," J. Microelectromech. Syst. 16, 491-499 (2007). 8. C. Hsien-Chang, C. Chao-Hung, I. F. Cheng, and L. Chi-Chang, "Manipulation of Bioparticles on Electrodeless Dielectrophoretic Chip Based on AC Electrokinetic Control," in Nano/Micro Engineered and Molecular Systems, 2007. NEMS '07. 2nd IEEE International Conference on (2007), pp. 1175-1178. 9. A. Jamshidi, P. J. Pauzauskie, P. J. Schuck, A. T. Ohta, P. Y. Chiou, J. Chou, P. D. Yang, and M. C. Wu, "Dynamic manipulation and separation of individual semiconducting and metallic nanowires," Nature Photonics 2, 85-89 (2008). 10. M. Hoeb, J. O. Radler, S. Klein, M. Stutzmann, and M. S. Brandt, "Light-induced dielectrophoretic manipulation of DNA," Biophys. J. 93, 1032-1038 (2007). 11. C. F. Chou, J. O. Tegenfeldt, O. Bakajin, S. S. Chan, E. C. Cox, N. Darnton, T. Duke, and R. H. Austin, "Electrodeless dielectrophoresis of singleand double-stranded DNA," Biophys. J. 83, 2170-2179 (2002). 12. A. T. Ohta, A. Jamshidi, J. K. Valley, H.-Y. Hsu, and M. C. Wu, "Optically actuated thermocapillary movement of gas bubbles on an absorbing substrate," Appl. Phys. Lett. 91, 074103-074103 (2007). 13. S. Grilli, and P. Ferraro, "Dielectrophoretic trapping of suspended particles by selective pyroelectric effect in lithium niobate crystals," Appl. Phys. Lett. 92, 3 (2008). #105882 $15.00 USD Received 5 Jan 2009; revised 12 Mar 2009; accepted 14 Mar 2009; published 18 Mar 2009 (C) 2009 OSA 30 March 2009 / Vol. 17, No. 7 / OPTICS EXPRESS 5232 14. J. K. Valley, A. Jamshidi, A. T. Ohta, H. Hsan-Yin, and M. C. Wu, "Operational regimes and physics present in optoelectronic tweezers," J. Microelectromech. Syst. 17, 342-350 (2008). 15. M.P. Hughes, Nanoelectromechanics in engineering and biology, CRC Press (2003). 16. A. R. Minerick, R. H. Zhou, P. Takhistov, and H. C. Chang, "Manipulation and characterization of red blood cells with alternating current fields in microdevices," Electrophoresis 24, 3703-3717 (2003). 17. M. Hollinshead, A. Vanderplasschen, G. L. Smith, and D. J. Vaux, "Vaccinia virus intracellular mature virions contain only one lipid membrane," J. Virol. 73, 1503-1517 (1999).


Introduction
Optoelectronic tweezers (OET) are a relatively new micromanipulation technique that complement the established techniques of Optical Tweezers (OT) and Dielectrophoresis (DEP).OET uses the dielectrophoretic force to manipulate particles, however, the metal electrodes, that are usually used in DEP, are replaced by 'virtual electrodes' created by selectively illuminating areas of a thin film of photoconductive material (see Fig. 1).These virtual electrodes can then be moved giving the real time continuous control of particles.This continuous control is one of the features of OT that has made it a popular technique however OT requires very high light intensities.The reduced intensity necessary for OET give us the ability to perform massively parallel manipulation of microparticles, with the creation of up to 15,000 traps at once [1], and to create traps with stiffness similar to OT whilst using 500 times less optical power [2].It has been shown that OET can be used to manipulate live cells [3][4][5][6][7][8], nanowires [9], DNA [10,11] and even gas bubbles [12].Another method for DEP trapping with no metal electrodes is by using the pyroelectric effect [13].
Previous work has shown that the properties of the trapped objects are among the factors determining the size and shape of the OET trap.In this report we will specifically look at the properties of cell traps.Using positive DEP, where particles are attracted to the high field region, increasing the power of the incident light has been shown to increase both the strength and the size of an OET trap [3].In this work we will keep the illumination intensity constant whilst changing the size of the illuminated region and study the effect on the trapped cell.The maximum lateral distance at which cells are affected by the trap has been measured [4] however the forces experienced by cells within the trap have not previously been studied.To study these forces we exert drag forces on the trapped cell by moving the OET chamber with respect to the light pattern.This has been shown to allow the plotting of force versus position profiles of the trap, allowing the calculation of trap stiffness [2].The light patterns used in these experiments are created by focusing the pattern from a data projector through a 20x microscope objective.The optical pattern created is of constant intensity across the illuminated region which is very different from the Gaussian profile of a laser spot that has previously been studied [2].The intensity profile of the optical pattern is an important parameter as gradients in the optical image create gradients in the electrical field which in turn produce the force (see Section 3).The optical gradients used here are much lower than would be required for an OT which demonstrates the power of OET as a low intensity technique, not requiring tight focusing of the light, but retaining full control of the position and shape of the trap.Fig. 2. A trapped cell.Video frames from an experiment showing the light pattern, outlined in red, being used to trap the HeLa cell, outlined in white, at velocities A) and D) 10 µms -1 , B) and E) 30 µms -1 , C) and F) 50 µms -1 .A) to C) show a spot diameter of 73µm, where the optical spot is large compared to the size of the cell it moves to the edge of the optical spot even at low velocities.D) to F) show a 12µm diameter spot, here the trap is similar size to the cell and it stays central at low velocities.X is the lateral distance between the center of the cell and the center of the trap, which are indicated by vertical lines.

Experimental results
The optical setup consisted of a data projector (2400MP, DELL) focused into a 20x objective (N.A. 0.42) within a microscope constructed from Thorlab's 30mm cage system.The OET device consisted of two ITO plates spaced 120µm apart, the bottom plate being coated with 1µm of hydrogenated amorphous silicon.An A.C. voltage of 20V peak to peak was applied at 100kHz.When the OET chamber is moved with respect to the light pattern the drag force on the trapped particle is countered by the DEP force.Due to the cell's low mass the force required to accelerate the cell is insignificant compared to the viscous drag forces e.g. a force of just 1.67x10 -15 N could accelerate the cell to 50µms -1 in the time between video frames if there were no viscous forces present.This compares to viscous forces in the pN range.By increasing the velocity whilst measuring the particle's position within the trap until the particle falls out of the trap, a force profile can be found [2].However, as will be discussed later, the force profile of the trap varies with height above the amorphous silicon (a-Si) photoconductor.As the particle rests on the lower surface this means that different size particles will experience not just different magnitude forces but different shape traps.Therefore, the force profile should be designed for the particle of interest.One of the main proposed applications for OET is in the manipulation of biological cells so in this work we use a standard cell type, HeLa cells, as analyte.Figure 2 shows HeLa cells being trapped by a 73µm diameter trap and a 12µm trap at increasing velocities.It can be seen that when the trap is much larger than the cell, the cell does not experience any force in the centre of the trap but moves to the edge of the trap even at low velocities.When the trap is similar in size to the cell, the particle is only moved from the centre of the trap at higher velocities.The force on the cell can be calculated from Stokes drag force.As the cell is close to a surface we must also use Faxen's correction [2].The cell is assumed to be spherical and touching the lower surface.The viscosity of the liquid is assumed to be 1mPas.The position of the particle is measured by recording a video of the particle being trapped at a constant velocity and then using particle tracking software to measure the displacement from the centre of the trap.We can the plot the force verses position profile of each trap, as shown in Fig. 3.
Figure 3 shows that the cells could be trapped at up to 60µms -1 (26.2 pN) with the 31µm diameter trap and 50µms -1 (21.8 pN) with the other spots.For small displacements the larger traps exhibit zero force on the particle.This is not good for fine control of the cell as the trap must be moved far before the cell feels any force.The 12µm diameter trap shows that it is possible to create a trap with constant trap stiffness, which will allow fine control over the position of the cell, the stiffness created was 3x10 -6 Nm -1 .Figure 3 also shows the results of simulations where the cell is assumed to be 3.5µm from the surface, the reasons for this distance are discussed in section 3. Fig. 3. Trap profile results.The experimental results are plotted as crosses and compared to numerical simulations plotted as solid lines for optical spot diameters 73µm (black), 49µm (red), 31µm (green), 12µm (blue).

Simulations
It has been shown that the electrical field gradients created in OET can be numerically modeled by using the Finite Element Method (FEM) [2,3].Here the electric fields within the OET chamber are modeled and the resulting gradient of the square of the electric field is used to find the DEP force using the well-know equation [6]; where m ε is the permittivity of the medium, Re[k(ω)] is the Clausius-Mossotti factor and 2 E ∇ is the gradient of the electrical field squared [6].The Clausius-Mossotti factor is given by; ε are the complex permittivities of the particle and medium respectively and ω is the angular frequency.The relation between the real permittivity and complex permittivity is given by ) / ( * ω σ ε ε j − = with ε the permittivity and σ is the conductivity.Whilst this approach gives good agreement between experimentally measured forces and the simulated forces for small homogeneous particles the modeling of cells is more complicated as the structure of a cell is more complicated.A reasonable approximation can be made by modeling the cell as a single shell of low conductivity material surrounding a highly conductive core [15,16]. Dielectrophoretic force can either be positive (pDEP), towards areas of high electric field, or negative (nDEP), away from these areas depending on the sign of the Clausius-Mossotti factor.A cell will experience nDEP at low frequencies and pDEP at high frequencies crossing over at a certain frequency.This crossover frequency is dependent on the physical properties of the suspending medium and the cell.Some of the properties can be deduced by measuring this frequency in liquid of varying conductivity [16].Cultured HeLa cells were removed from a substrate with trypsin before being re-suspended in an isotonic sugar solution [1].By adding varying amounts of cell culture medium to the solution the conductivity was varied from 3.2 to 2.2x10 -2 Sm -1 and the crossover frequency was measured as shown in Fig. 4. The range of conductivities that can be explored is limited by the response of the OET chamber with higher conductivities giving a reduction in the force [14].The trap stiffness experiments shown earlier in this paper were performed in 1x10 -2 Sm -1 media at 100kHz giving a positive force.Fig. 4. Simulations of DEP.The experimentally measured crossover frequency is plotted as points along side simulated lines for a conductive core (radius 7.5µm, relative permittivity 50, and conductivity 0.53Sm -1 ) surrounded by a single shell of; A. Conductivity 0.87µSm -1 , thickness varying from 20 to 1 nm from top to bottom curves.B. Thickness 5nm and conductivity varying from 0.1 to 3µSm -1 from top to bottom curves.The relative permittivity of the medium was assumed to be 78.
These experimental results were compared to a single shell model in MATLAB.The effect of the cell's insulating shell is accounted for by defining an effective permittivity for the core-shell structure that is given by [15]; Where * 1 ε and * 2 ε are the complex permittivities of the core and shell, respectively, and 1 r and 2 r are their radii.This effective permittivity can then be used in place of the particle's and thickness of 5nm gave the best agreement.This is consistent with value from the literature where the membrane thickness has been measured as 4.97±0.2nmby scanning electron microscopy [17].These values give us a Clausius-Mossotti factor of 0.8 at 1x10 5 Hz, the frequency the experiments were performed at.To calculate the force it is now just necessary to find 2 E ∇ .3D numerical simulations were performed in COMSOL Multiphysics (COMSOL) with the light spot being modeled as a saturated Gaussian that matches the profiles measured from the optical spots from the experiments.These saturated Gaussians are entered into COMSOL as profiles in conductivity of the a-Si varying from 1x10 -6 Sm -1 to 1.5x10 -4 Sm -1 the dark and illuminated conductivities of the a-Si, respectively (see Fig. 5).These are the conductivities measured for an illumination power of 2.5Wcm -2 , the intensity produced by our data projector when reduced through a 20x objective.
The size and shape of the optical profiles were measured by recording video of the trap and analyzing a frame in a graphics package (National Instruments Vision Assistant).As the images of the traps were projected in red, the red pixel value is used as a measure of the projected light intensity.The optical profiles are shown in Fig. 5(a) and it can be seen that the measured optical profiles match well with saturated Gaussian curves.It has been shown that the conductivity of the photoconductor is proportional to the illumination intensity [14].Fig. 5. Measured light patterns and simulated forces.A) The profile of the optical spots are plotted as crosses and compared to saturated Gaussian curves plotted as solid lines for optical spot diameters 73µm (blue), 49µm (red), 31µm (green), 12µm (black).B) Simulated force profile at 1µm (green curve, reduced by a factor of 10), 3.5µm (red curve, reduced by a factor of 4) and 7.5µm (blue curve, reduced by a factor of 4/3) above the surface of the OET chamber.The experimental points for the 31µm trap are shown as crosses.
Figure 5(b) shows the experimental results for the 31µm diameter trap plotted against the simulated forces at different heights from the a-Si surface.As the cells have a diameter of 15µm, one might expect the simulated forces at 7.5µm from the surface to match the experimental data.However it was found that the shapes of the experimental and simulated traps at this height did not match well.
Using 2 E ∇ at 7.5µm above the a-Si surface is effectively assuming that this value does not vary over the volume of the cell or that the higher forces near the OET surface are cancelled by the lower forces further from the a-Si so that the value half way up the cell is a good average.This is not a valid assumption in this case as the cell is large compared to the distance over which the gradients vary and the force profile closer to the surface is different from the profile further from the a-Si so that the two effects do not cancel each other.From Fig. 5(b) it can be seen that the force profile at 7.5µm from the surface is close to an ideal trap profile where the force is proportional to the displacement, a trap that could be represented with a constant spring constant over the central part of the trap.This means that once the trap is moved a small amount the trapped particle would immediately feel a restoring force pushing it towards the center of the trap.The experimental results however show a trap that is completely flat in the centre and the particle does not experience a force until it is several microns from the centre of the trap.This is similar to the force profiles lower down in the trap, closer to the a-Si.Here 2 E ∇ is greater so to compare the trap profiles the magnitude of the gradients has to be reduced.It was found that the profile at 3.5µm above the a-Si matched the measured trap profiles surprisingly well when reduced by a factor of four for all four traps, as shown in Fig. 3.As the cell is not visibly deformed the extra force can't be due to the cell being pulled closer to the surface although this may have a small effect.This factor of four is simply the amount that the force at 3.5µm from the a-Si is greater than the average force over the volume of the cell.The simplicity of this modeling method gives us an incredibly useful tool for future experiments allowing us to accurately predict the force on the cell and hence design light patterns that will translate into the patterns of force we desire.One such desirable pattern is a simple ideal trap with a constant spring constant as mentioned earlier.Figure 3 shows that this can be produced by reducing the size of the optical spot to 12µm.Although reducing the illuminated area to this extent does reduce the maximum force that can be exerted by 17%, the stiffness is the same.In this case the trap stiffness is 3x10 -6 Nm -1 and the maximum velocity that can be achieved is 50µms -1 with the 12µm trap.An analytical expression can be fitted to the numerically simulated force profile.It was found that a 3 rd order polynomial could be fitted with good agreement.This expression can then be integrated to find the work done by the trap on the particle which can then be used to find the potential energy.Figure 6 shows the potential energy landscape that the particle experiences due to the trap.The depth of the potential well is found to be 1.6x10 -16 J (4x10 4 k B T).

Conclusions
In this paper we have shown that the profile of an OET trap varies with height from the a-Si surface.This is important since different size particles will experience different trap profile from highly non-linear traps to ideal traps that can be described by a single spring constant.We find that the behavior of HeLa cells can be predicted through numerical modeling with surprising accuracy by simply considering the force at 3.5µm from the surface, reduced by a factor of four.Using this information, we show that it is possible to create an ideal trap where the force is proportional to the distance from the trap centre for HeLA cells in a 12µm trap.The resulting trap has a spring constant of 3x10 -6 Nm -1 , and is capable of moving the cell up to 50µms -1 .We also calculate the potential energy well that the cell experiences due to the trap and find the depth of the well to be 1.6x10 -16 J.

Fig. 1 .
Fig. 1.The OET device.Schematic showing the OET chamber consists of a top contact and a bottom photoconductive substrate with an illuminated region.

#
105882 -$15.00USD Received 5 Jan 2009; revised 12 Mar 2009; accepted 14 Mar 2009; published 18 Mar 2009 (C) 2009 OSA permittivity in Equation 2[6], allowing us to calculate the Clausius-Mossotti factor.Using values for the permittivity and conductivity of the cell's cytoplasm and membrane from the literature crossover frequency for varying medium conductivities can be modeled.It was found that to fit the simulations with the experimentally measured crossover frequencies a A B #105882 -$15.00USD Received 5 Jan 2009; revised 12 Mar 2009; accepted 14 Mar 2009; published 18 Mar 2009 (C) 2009 OSA

Fig. 6 .
Fig. 6.Potential energy well.The potential energy well experienced by the trapped HeLa cell is calculated from the force verses position profile for the 12 µm diameter trap.The depth of this potential energy well is 1.6x10 -16 J.