Electroosmotic Flow and DNA Electrophoretic Transport in Micro/Nano Channels
Dissertation
Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School of The Ohio State University
By
Lei Chen, M.S.
Graduate Program in Mechanical Engineering
The Ohio State University
2009
Dissertation Committee: A.T. Conlisk, Advisor Joseph Heremans Vishwanath V. Subramaniam Sandip Muzumder c Copyright by
Lei Chen
2009 ABSTRACT
In Micro/nano fluidic systems, electrokinetic transport is a convenient method to move materials, such as water, ions and particles for fast, high-resolution and low-cost anal- ysis and synthesis. It has wide applications to drug delivery and its control, DNA and biomolecular sensing, manipulation, the manufacture laboratories on a microchip (lab-on- a-chip) and many other areas. In the present work, electrokinetically driven fluid flow and particle transport in micro/nanoscale channels/pores with heterogeneous surface potential or converging shape are investigated theoretically and numerically.
A step change in wall potential is found to induce a recirculation region in the bulk electroosmotic flow and interesting flow structures can be achieved by manipulating the surface heterogeneous patterns. Most of the previous work on this problem is based on the
Debye-Huckel approximation and the validity of Boltzmann distribution for ionic species.
In the present work, ionic species distributions in the electric double layers are found to be different from the Boltzmann distribution and this deviation is more noticeable for higher applied electric field.
A mathematical model is developed to simulate the electroosmotic flow (EOF) and the transport of embedded particles in micro/nano nozzles/diffusers. Results can be used to estimate the mass transport of charged/uncharged species in micro/nano nozzles/diffusers which has a potential application in transdermal drug delivery. The model is extended to investigate the DNA electrophoretic transport through a converging nanopore for the
ii purpose of DNA sequencing. The flow field, the resistive forces acting on the DNA, the
DNA velocity and the ionic current through the nanopore are calculated numerically based on the Poisson-Boltzmann theory and the lubrication approximation. It is found that the electroosmotic flow inside the nanopore plays an important role in the DNA translocation process and the resulting viscous drag decreases the effective driving force acting on the
DNA substantially. Entropic forces, used to be considered as the main resistive forces in previous works, are found to be small and negligible. Modeling and simulations are validated by the good agreement with the experimental data for the tethering force and the
DNA velocity.
iii ACKNOWLEDGMENTS
I am deeply grateful to Professor A. T. Conlisk, whose help, stimulating suggestions and encouragement helped me in all the time of research for and writing of this thesis. I would also like to thank Prof. Vishwanathan Subramaniam, Prof. Sandip Muzumder and
Prof. Joseph Heremans for accepting to be in the my Docotoral Examination Committee and for their valuable comments and advises on my thesis work. I would also like to thank
Subhra for his discussions and supporting on my research.
This work is supported by NSEC, Center for Affordable Nanoengineering of Poly- meric Biomedical Devices (NSF Grant No. EEC-0425656) and the PI is Professor James
Lee from Department of Chemical and Biomolecular Engineering at OSU. The author is grateful for their support for this work.
Especially, I would like to give my special thanks to my parents and my husband Xi- aoyin for their encouraging and support.
iv VITA
April 9, 1979 ...... Born - Beijing, China
July 2001 ...... B.S. Mechanical Engineering Department of Thermal Engineering Tsinghua University, Beijing, China August 2001 - June 2003 ...... Graduate Engineer HVAC Lab Tsinghua University Beijing, China September 2003-Jan 2006 ...... M.S. Mechanical Engineering Department of Mechanical Engineering The Ohio State University Columbus, OH Feb 2006 - present ...... TheOhio State University Columbus, OH
PUBLICATIONS
Research Publications
Lei Chen and A. T. Conlisk, “Modeling of DNA translocation in nanopores”, 47th AIAA Aerospace Sciences Meeting, AIAA paper 2009-1121, 2009
Lei Chen and A. T. Conlisk, “Effect of nonuniform surface potential on electroosmotic flow at large applied electric field strength”, Biomedical Microdevices, Vol. 11, 251-258, 2009
Lei Chen and A. T. Conlisk, “Electroosmotic flwo and particle transport in micro/nano nozzles and diffusers”, Biomedical Microdevices, Vol. 10, 289-298, 2008
v Lei Chen, P. Gnanaprakasam and A. T. Conlisk, “Electroosmotic Flow in Micro/nano Nozzles”, 45th AIAA Aerospace Sciences Meeting and Exhibit, AIAA paper 2007-0931, 2007
Lei Chen and A. T. Conlisk, “Generation of nanovortices in electroosmotic flow in nanochan- nels”, 4th AIAA Theoretical Fluid Mechanics Meeting, AIAA paper 2005-5057, 2005
Lei Chen and A. T. Conlisk, “Electroosmotic flow in annular diatoms”, 43rd AIAA Aerospace Sciences Meeting and Exhibit, AIAA paper 2005-0674, 2005
FIELDS OF STUDY
Major Field: Mechanical Engineering
Studies in Fluid Mechanics, Mass Transer and Electrokinetics: Professor A.T. Conlisk
vi TABLE OF CONTENTS
Page
Abstract ...... ii
Acknowledgments ...... iv
Vita...... v
List of Tables ...... x
List of Figures ...... xii
Chapters:
1. Introduction ...... 1
1.1 Background ...... 1 1.2 Electrokinetic Phenomena ...... 4 1.3 The Molecular Biology of DNA ...... 7 1.4 Applications ...... 11 1.5 Literature Survey ...... 20 1.5.1 Electroosmotic Flow ...... 20 1.5.2 DNA Transport ...... 23 1.6 Present Work ...... 28
2. Governing Equations and General Assumptions ...... 30
2.1 Introduction ...... 30 2.2 Assumptions ...... 30 2.3 Governing Equations ...... 32 2.3.1 Poisson Equation ...... 32 2.3.2 Nernst-Planck Equation ...... 34
vii 2.3.3 Navier-Stokes Equation ...... 37 2.4 Dimensional Analysis ...... 38 2.5 Boltzmann Distribution and the Debye-H¬uckel Approximation ...... 43 2.6 Ionic Current and Current Density ...... 45 2.7 Summary ...... 46
3. Electroosmotic Flow with Nonuniform Surface Potential ...... 47
3.1 Introduction ...... 47 3.2 Governing Equations and Boundary Conditions ...... 48 3.3 Analytical Solution for Potential Based on the Debye-H¬uckel Approxi- mation ...... 51 3.4 Numerical Methods ...... 55 3.5 Numerical Results ...... 60 3.6 Ionic Transport in a Nanopore Membrane with Charged Patches ..... 67 3.7 Summary ...... 74
4. Electroosmotic Flow and Particle Transport in Micro/Nano Nozzles/Diffusers . 76
4.1 Introduction ...... 76 4.2 Governing Equations and Boundary Conditions ...... 80 4.3 Lubrication Solutions for EOF in the Debye-H¬uckel Limit ...... 82 4.3.1 Known Surface Potential ...... 82 4.3.2 Known Surface Charge Density ...... 86 4.4 Results for EOF Based on the Debye-H¬uckel Approximation ...... 88 4.4.1 Results for EOF: Known Surface Potential ...... 88 4.4.2 Results for EOF: Known Surface Charge Density ...... 97 4.5 Analysis Based on the Thin Double Layer Limit ( 1)...... 100 4.6 Numerical EOF Calculation Based on Poisson-Boltzmann Model . . . . 108 4.7 Numerical Solution using COMSOL ...... 115 4.8 Particle Transport ...... 120 4.9 Particle Transport Using Two-Phase Flow Modeling ...... 125 4.10 Modeling the Fabrication of the Nanonozzles ...... 129 4.11 Summary ...... 134
5. The Forces that Affect DNA Translocation ...... 137
5.1 Introduction ...... 137 5.2 Governing Equations and Boundary Conditions ...... 139 5.3 Electroosmotic Flow in an Annulus ...... 141 5.4 Lubrication Solutions for the Flow Field in a Conical Nanopore: The Viscous Drag Force Acting on the DNA inside the Nanopore ...... 146
viii 5.5 Uncoiling/recoiling Forces Acting on the DNA ...... 148 5.6 Viscous Drag on the DNA Outside the Nanopore ...... 154 5.7 Numerical Validation of the Lubrication Model using Tethering Force Data157 5.8 Summary ...... 166
6. DNA Translocation Velocity ...... 169
6.1 Introduction ...... 169 6.2 Governing Equations and Boundary Conditions ...... 171 6.3 Solution for the Flow Field ...... 172 6.4 Asymptotic Solution for potential φ and the EOF velocity ueof ...... 174 6.5 Results for DNA Translocation Velocity ...... 178 6.6 DNA Entering and Leaving the Nanopore ...... 185 6.7 Modeling DNA Delivery into a Cell through the Nanonozzles ...... 198 6.8 Summary ...... 201
7. Summary ...... 204
Appendices ...... 211
Appendices:
A. Discretising the Equations ...... 211
A.1 Introduction ...... 211 A.2 Discritization ...... 211 A.3 The SIMPLE Algorithm ...... 231 A.4 Residuals ...... 232
ix LIST OF TABLES
Table Page
2.1 The scales used to nondimensionalize the equations listed in the previous section...... 39
3.1 The numerical accuracy check for h =20nm channel with over-potential φe =0.2. The results of 400 × 80 grid are compared with the results of 800 × 160 grid and two-digit accuracy has been attained...... 58
4.1 Values of parameters used for simulations. Note that the concentration of species shown here are reservoir concentrations...... 90
4.2 Parameters used for the numerical calculations for known surface charge. Note that the concentration of species shown here are reservoir concentra- tions...... 97
4.3 The numerical accuracy check for the numerical calculation for the elec- troosmotic flow in the nano-diffuser shown in Figure 4.16. The results for 81×161 grid are compared with the results of 161×321 grid and two-digit accuracy has been attained. The numerical results are shown in Figure 4.17. Note that here the pressure p is rescaled on EI,x...... 112
4.4 The numerical accuracy check for the results calculated using COMSOL shown in Figure 4.21. The mesh I has 62, 442 triangular elements and the mesh II is a finer mesh, having 183, 712 triangular elements. The results shown for the electric potential φ∗(mV ), the velocity u∗(mm/s) and pressure p∗(Pa) are in dimensional form...... 118
4.5 The dependence of the particle motion on the wall charge and the particle charge. Here uF is the bulk motion of the fluid flow and uEM is the ve- locity of particle due to electrical migration. The particle velocity is the summation of these two terms...... 123
x 5.1 Comparison between the magnitudes of the uncoiling and recoiling force ∗ ∗ (FU and FR) calculated based on different models and the electrical driving ∗ force (Fe )...... 154
5.2 Comparison between the magnitudes of the viscous drag force acting on the DNA blob-like configuration outside the nanopore and the electrical ∗ driving force (Fe ). The viscous drag due to the approaching of the DNA ∗ blob to the nanopore is defined as Fblob1 and the viscous drag force acting on the DNA blob in Reservoir I due to the flow discharged from the nanopore ∗ is defined as Fblob2...... 157
5.3 The numerical accuracy check for the numerical calculation based on the experimental parameters (Keyser et al., 2006) (results are shown in Figure 5.10). The results for 121×161 grid (Mesh I) are compared with the results of 241 × 321 (Mesh II) grid and two-digit accuracy has been attained. Note that the pressure shown here is rescaled on EI,x)...... 162
5.4 Comparison between the magnitudes of the forces that may affect DNA translocation, calculated for a 16.5kbps double stranded DNA through a nanopore with a radius of 5nm at the small end. The applied voltage drop is 120mV . Note that the viscous drag force acting on the DNA inside ∗ the nanopore (Fd ) listed here is calculated for a DNA immobilized in the ∗ =0 nanopore (uDNA )...... 167
6.1 Comparison for DNA velocity between the numerical results (VNum in ∗ 2 m/s) and the experimental data (VExp in m/s). Here σw (in C/m )is the surface charge density of the nanopore; LDNA is the length of the DNA used in the experiments; CKCl is the concentration of the KCl solution. . . 182
xi LIST OF FIGURES
Figure Page
1.1 The required pressure drop and voltage drop for nanochannels with dif- ferent channel height. The required flowrate is 1µl/min; the length of the channel is 3.5µm and the width of the channel is 2.3µm. The zeta potential of the channel wall is −11mV and the working fluid is 1M NaCl (Conlisk, 2005)...... 2
1.2 Micro and nanonozzles: (a) The Polymethyl Methacrylate (PMMA) mi- cronozzles used to characterize fluid flow and particle transport. (b)The Polymethyl Methacrylate (PMMA) nanonozzles fabricated using sacrifi- cial template imprinting (STI) methods (Wang et al., 2005)...... 3
1.3 An illustration of the ion distribution near a glass surface and the electric double layer (EDL). The dash line denotes the approximate edge of the electric double layer...... 4
1.4 Electric double layer around a charge particle...... 6
1.5 DNA double-helix structure showing base pairs as vertical lines. The di- ameter of the spiral of the helix is about 2nm ...... 8
1.6 Conformation of DNA confined in a pore. Here h∗ is the half height of the channel of the radius of the pore; Rg is the radius of gyration of the DNA and lp is the persistence length of DNA...... 10
1.7 Images of microneedles used for transdermal drug delivery. (a) Solid mi- croneedles (150µm tall) etched from a silicon wafer were used in the first study to demonstrate microneedles for transdermal delivery (Henry et al., 1998). (b) Microneedle arrays (up to 300µm tall) in standard silicon wafer using potassium hydroxide (KOH) wet etching (Wilke et al., 2005). . . . . 12
xii 1.8 Schematic drawing of a transdermal drug delivery system. There is a power source in the back. The foam ring with the drug reservoir is placed under the cathode and a microprojection array is integrated into the patch(Lin et al., 2001)...... 14
1.9 Gene delivery using nanotip array (Fei, 2007) under an applied electric field. 15
1.10 Concept for sequencing DNA by using a single protein pore. A single- stranded DNA (or RNA) molecule moves through the pore in the trans- membrane electric field. As it passes a ”contact site” each base produces a characteristic modulation of the amplitude in the single channel current (Alper, 1999)...... 18
1.11 Nanopores used in the experiments to analyze DNA: (a) natural α-hemolysin nanopore and (b) synthetic nanopore...... 19
1.12 The geometry of the nanopore and the DNA...... 26
2.1 Cartesian coordinates and cylindrical coordinates used for modeling. . . . 32
3.1 The geometry of the channel with a patch with over-potential φp...... 48
3.2 Analytical results for potential calculated from equation (3.23). The zeta potential of the channel wall and the patch is −1.3 mV and 6.5 mV respec- tively. The electrolyte solution is (a) 0.1M NaCl and (b) 0.001M NaCl, and the EDLs are not overlapped in (a) and overlapped in (b). The height of the channel is (a) h =20nm and (b) h =50nm...... 54
3.3 A co-located mesh used for the numerical calculation...... 56
3.4 The flowchart of the numerical calculation using SIMPLE algorithm. . . . . 59
3.5 Numerical results for potential, mole fractions calculated from the Poisson- Boltzmann equations. The imposed electric field is 106V/m. The height of the channel is 20nm and the wall surface potential is −12mV . The electrolyte is 0.1MNaClsolution far upstream and the over-potential φp =2.0. The corresponding dimensionless parameters are =0.05, Λ=0.77, A =1and Pe=0.36...... 60
xiii 3.6 Numerical results for potential, mole fractions, electric field lines and stream- lines. The imposed electric field is 106 V/m. The height of the channel is 20 nm and the wall surface potential is −12 mV . The electrolyte is 0.1 M NaCl solution far upstream and the over-potential φp =2.0. The corre- sponding dimensionless parameters are =0.05, Λ=0.77, A =1and Pe=0.36...... 61
3.7 Numerical results for potential, mole fractions and electric field lines calcu- lated from the Poisson-Nernst-Plank equations. The imposed electric field is 107 V/m. The height of the channel is 20nm and the dimensional surface potential of the walls is −12 mV . The electrolyte is 0.1 MNaClsolution far upstream and the dimensional patch potential is 40 mV (φp =2.0). The corresponding dimensionless parameters are =0.05, Λ=7.7, A =1and Pe=3.6...... 63
3.8 Cation mole fraction contours and streamwise velocity contours calculated from the Poisson-Nernst-Plank equations for Pe =3.6 and Pe =0. The filled colors show the results for Pe =3.6 and the lines represent the results for Pe =0. The imposed electric field is 107 V/m. The height of the chan- nel is 20 nm and the dimensional surface potential of the walls is −12 mV . The electrolyte is 0.1 MNaClsolution far upstream and the dimensional patch potential is 40 mV (φp =2.0). The corresponding dimensionless parameters are =0.05, Λ=7.7, A =1...... 65
3.9 The dimensionless distance of the vortex center away from the patch as a function of the dimensionless over-potential of the patch. The center of the vortex is the point where u =0,v =0in the bulk flow. The imposed electric field is 107 V/m. The height of the channel is 20 nm and the dimensional surface potential of the walls is −12 mV . The electrolyte is 0.1 MNaClsolution far upstream and the corresponding dimensionless parameters are =0.05, Λ=7.7, A =1...... 66
3.10 Streamlines for the EOF in a nanochannel with a patch having negative over-potential. The height of the channel is 20nm and the aspect ratio of the patch A =1. The over-potential φe = −0.4 and the imposed electric field is 107V/m. The solution far upstream is 0.1MNaClelectrolyte buffer. 67
3.11 Streamlines for EOF in nanochannels with a patch having over-potential in the form of φp = φacosnπx. The zeta potential of the wall is −1mV and the imposed electric field is 106V/m. The solution far upstream is 0.1M NaCl...... 68
xiv 3.12 The geometry of the reservoir and nanopore system (a) and the mesh used in the simulations (b). The height of the nanopores is 8 nm and the length of the pores is 1 µm. The size of the feed and the permeate reservoirs is 10 µm by 10 µm. The length of the patch is 500 nm...... 69
3.13 Results for the electric potential, streamlines and concentration of albumin: (a) electric potential (b) streamlines and concentration contours of charged species. The ionic strength of the electrolyte solution is 0.14 M and the zeta potential of the charged patch is −46 mV . The concentration of the charged species (albumin) is 0.6 mM. The valence of albumin is −17 and the diffusion coefficient is 10−10 m2/s. The height of the nanopores is 8nm and the length of the pores is 1 µm. The size of the feed and the permeate reservoirs is 10 µm by 10 µm. The length of the patch is 500 nm...... 71
3.14 Results for albumin concentration at the centerline of the channel (a) with charged patch (b) without charged patch. The ionic strength of the elec- trolyte solution is 0.14 M and the zeta potential of the charged patch is −46 mV . The concentration of the charged species (albumin) is 0.6 mM. The valence of albumin is −17 and the diffusion coefficient is 10−10 m2/s. The height of the nanopores is 8 nm and the length of the pores is 1 µm. The size of the feed and the permeate reservoirs is 10 µm by 10 µm. The length of the patch is 500 nm...... 73
4.1 The geometry of the micro-diffuser used in the experiments. The particle donor reservoir is on the right and the particle receiver reservoir is on the left of the diffuser. The walls are negatively charged and the positive elec- trode is placed in the particle receiver (left). The flow is from left to right and thus it is a diffuser. The polystyrene beads are shown and their motion is from right to left in a direction opposite to the bulk flow...... 77
4.2 The geometry of the reservoir-nozzle system (a) the system used in MD simulations (b) the cooresponding dimensions. Here only three stripe of water molecules are shown (in grey and red) and Cl− (green) and Na+ (blue) ions dispersed in water are shown...... 78
4.3 The geometry of the nozzle/diffuser used for modeling. The walls are neg- atively charged and for the polarity shown here, the flow is from left to right in both the nozzle and diffuser...... 79
xv 4.4 The surface charge density boundary condition for the potential equation is ∂φ = ± ± ( ) ∂y σy, h x , where σy is the component of the dimensionless surface charge density in the y direction ...... 87
4.5 The dimensionless imposed electric field (EI,x) for a nozzle (Figure 4.3 a) and a diffuser (Figure 4.3 b). The length of the nozzle is 65nm; the height at the inlet is 13nm and the height at outlet is 8nm. The length of the diffuser is 650 µm; the height at the inlet is 20 µm and the height at the outlet is 130 µm. The cross-section of the nozzle is rectangular, with the width in the direction into the paper much larger than the height. The potential drop over the length of the nozzle/diffuser corresponds to 80V/cm. The solid line shows the electric field for a channel (80V/cm)...... 89
4.6 Results for electroosmotic flow in the converging nanonozzle. The height of the nozzle is 13 nm at the inlet and 8 nm at the outlet and the length of the nozzle is 65 nm; =0.19 and the EDLs are overlapped. The imposed electric field is assumed to be 8000 V/mand the ζ-potential of the walls is −5mV . The pressure is zero both at inlet and the outlet and so the pressure shown here is induced by the presence of the EDLs...... 92
4.7 Results for electroosmotic flow in the experimental micro-diffuser. The height of the diffuser is 20 µm at the inlet and 130 µm at the outlet and the length of the diffuser is 650 µm; =5× 10−5 and the EDLs are thin compared to the diffuser. The imposed electric field is 8000V/m and the ζ-potential is −15mV . The pressure is zero both at the inlet and the outlet and so the pressure shown here is induced by the presence of the EDLs. . . 93
4.8 Results for pressure-driven flow in the experimental micro-nozzle as shown in Figure 4.3 (a). The height of the small end is 20 µm and 130 µm at the large end and the length of the nozzle is 650 µm. For (a), (b), (c), the dimensionless pressure is 1 at the large end and zero at the small end. For (d), (e), (f), the dimensionless pressure is 0 at the large end and 1 at the small end. The electroosmotic flow component is zero...... 95
4.9 Streamwise velocity and the streamlines for the electroosmotic flow in a micro-nozzle with adverse pressure gradient. The height of the small end is 65µm and 130µm at the large end and the length of the nozzle is 650µm. The applied voltage drop is 80V/cm and the zeta potential of the wall is −5mV . The dimensionless pressure is pi =0and po =10...... 96
xvi 4.10 Results for electroosmotic flow in the nano-nozzle. The height of the noz- zle is 14.6 nm at the inlet and 4.6 µm at the outlet and the length of the diffuser is 15 nm; =0.065. The imposed electric field is 0.5V/nm and the surface charge density is −0.01C/m2. The pressure is zero both at the inlet and the outlet and so the pressure shown here is induced by the presence of the EDLs...... 98
4.11 Results for electroosmotic flow in the nano-nozzle. The height of the noz- zle is 4.6 nm at the inlet and 14.6 µm at the outlet and the length of the diffuser is 15 nm; =0.065. The imposed electric field is 0.5V/nm and the surface charge density is −0.01C/m2. The pressure is zero both at the inlet and the outlet and so the pressure shown here is induced by the presence of the EDLs...... 99
4.12 The surface potential as a function of x. The large end of the nozzle/diffuser is 14.6 nm height; the small end is 4.6 nm height and the length is 15 nm. For (a), the surface charge density is −0.07C/m2 and the electrolyte solu- tion is 0.1 MNaCl. For (b), the surface charge density is −0.01C/m2 and the electrolyte solution is 0.001 M...... 100
4.13 The geometry of the slowly varying channel with wall defined as h(x)= 1+δsin(wx). The walls are negatively charged and for the polarity shown here, the flow is from left to right...... 107
4.14 The streamlines for electroosmotic flow in the channel shown in Figure 4.13...... 108
4.15 Mapping from the original (x, y) coordinates to the (ξ,η) coordinates, which have a rectangular mesh, for easier numerical discretization...... 110
4.16 Results for electroosmotic flow in the nano-nozzle. The height of the noz- zle is 14.6 nm at the inlet and 4.6 nm at the outlet and the length of the diffuser is 15 nm; =0.065. The imposed electric field is 0.5V/nm and the surface charge density is −0.07C/m2. The pressure is zero both at the inlet and the outlet and so the pressure shown here is induced by the presence of the EDLs...... 113
xvii 4.17 Results for electroosmotic flow in the nano-diffuser. The height of the dif- fuser is 4.6 nm at the inlet and 14.6 nm at the outlet and the length of the diffuser is 15 nm; =0.065. The imposed electric field is 0.5V/nm and the surface charge density is −0.07C/m2. The pressure is zero both at the inlet and the outlet and so the pressure shown here is induced by the presence of the EDLs...... 114
4.18 The comparison of DH (based on the analysis in Section 4.3.2) and non- linear results for streamwise velocity component ueof and pressure p. The height of the nozzle is 4.6 nm at the inlet and 14.6 µm at the outlet and the length of the diffuser is 15 nm; =0.065. The imposed electric field is 0.5V/nm and the surface charge density is −0.07C/m2. The pressure is zero both at the inlet and the outlet and so the pressure shown here is induced by the presence of the EDLs...... 116
4.19 The comparison of DH and nonlinear results for streamwise velocity com- ponent ueof and pressure p. The height of the nozzle is 4.6 nm at the inlet and 14.6µm at the outlet and the length of the diffuser is 15nm; =0.065. The imposed electric field is 0.5V/nm and the surface charge density is −0.01C/m2. The pressure is zero both at the inlet and the outlet and so the pressure shown here is induced by the presence of the EDLs...... 117
4.20 The Mesh used for the simulations. Close to the walls, extra fine mesh is used and in the center region, the mesh is coarser...... 118
4.21 Results for ionic concentration and flow field for the EOF in the nozzle- reservoir system calculated using COMSOL. The parameters are listed in table 4.2...... 119
4.22 The results for EOF in the nozzle-reservoir system based on Molecular Dynamics simulations(Shin & Singer, 2009) and continuum theory. The height of the nozzle is 14.6 nm at the inlet and 4.6 µm at the outlet and the length of the diffuser is 15 nm; =0.065. The imposed electric field is 0.5V/nmand the surface charge density is −0.07C/m2...... 121
4.23 The reservoir-nozzle system used to analyze particle transport. The reser- voir on the right is the particle donor reservoir and the one on the left is the particle receiver reservoir. Negatively charged particles always tend to move from left to right due to electromigration (uEM < 0)...... 121
xviii 4.24 Comparison of the analytical results and the experimental data (Wang & Hu, 2007). The length of the micro-diffuser is 650 µm; the inlet height is 20 µm and the outlet height is 130 µm. The electric field is 80 V/cm. The ζ-potential of the PMMA walls are −15 mV (Kirby & Hasselbrink, 2004). Both the particles and the walls are negatively charged and so uF > 0,uEM < 0...... 124
4.25 Process schematic of EOF based dynamic assembly of silica(Wang et al., 2005) ...... 131
4.26 The schematic diagram for the simplified models. (a) one-dimensional model: A homogeneous chemical reaction is happening in the bulk and a heterogeneous chemical reaction is happening at the liquid-solid interface at a different rate. The interface is moving with the deposition of the reac- tion product C. (b) two-dimensional model: EOF in rectangular channels with chemical reactions ...... 132
5.1 Schematic diagram of DNA translocation in a negatively charged nanopore. 137
5.2 The geometry of the nanopore and the DNA...... 140
5.3 The geometry of a DNA placed in a cylindrical nanopore. The computa- tional domain is the annular region between the DNA and the nanopore. . . 142
5.4 The results for potential and velocity calculated from equations 5.8 and 5.14. The radius of the nanopore is 5nm and the radius of the DNA is 1nm. The surface charge density is −0.015C/m2 on the DNA surface and −0.006C/m2 on the wall...... 144
= ∂u = 5.5 The dimensionless shear stress τ ∂r at r a as a function of electrolyte concentration. The radius of the nanopore is 5nm and the surface charge density is −0.06C/m2 while the radius of the DNA is 1nm and the surface charge density is −0.15C/m2. The electrolyte concentration varies from 0.01M ( =0.6)to0.1M ( =0.06)...... 144
5.6 The schematic diagram of the DNA-nanopore system used to investigate the uncoiling-recoiling force due to the DNA conformational change during the DNA translocation process (shown in (a)) and the resulting uncoiling- recoiling force (shown in (b)) for a 16.5kbps double stranded DNA. . . . . 149
xix 5.7 Schematic drawings of a DNA with one end dragged into a cylindrical nanopore. The radius of the nanopore is h∗ and the force acting on the DNA is f...... 152
5.8 The schematic diagram for (a) the experimental setup(Keyser et al., 2006) and (b) the force balance. The DNA is attached to a polystyrene bead and a tethering force Ft is used to immobilize the DNA. The electrical driving force Fe and the viscous drag force Fd are also shown...... 158
5.9 The geometry of the nanopore used in Keyser’s experiments. The nanopore is a 20nm long cylindrical nanopore with 5nm in radius connected to a 20nm long converging nozzle on one end and a 20nm long diverging dif- fuser on the other end. The slope of the converging or diverging wall is about 23◦...... 159
5.10 Numerical results for the potential, concentration of cations, streamwise velocity and the streamlines calculated based on the lubrication approxi- mation. The computational domain is the region between the DNA surface (left boundary) and the nanopore wall (right boundary). The nanopore is a 20nm long cylindrical nanopore with 5nm in radius connected to a 20nm long converging nozzle on one end and a 20nm long diverging diffuser on the other end. The surface charge of the nanopore is −0.06C/m2 and the applied voltage drop is 120mV . The electrolyte solution is 0.1MKCl. . . 161
5.11 Comparison of tethering forces between the experimental data, the numeri- cal results based on the lubrication approximation and the numerical results using COMSOL. The electrolyte solution is 0.1MKCl. Experimental data are recorded for one molecule at three different distances from the nanopore (filled squares 2.1µm, filled circles 2.4µm, filled triangles 2.9µm(Keyser et al., 2006). As shown in Figure 5.8, Lb is the distance between the bead and the nanopore...... 163
5.12 The geometry of the nanopore, the mesh and the numerical results for the concentration of anions, streamwise velcity. The nanopore is a 20nm long cylindrical nanopore with 5nm in radius connected to a 20nm long con- verging nozzle on one end and a 20nm long diverging diffuser on the other end. The surface charge of the nanopore is −0.06C/m2 and the applied voltage drop is 120mV . The electrolyte solution is 0.1MKCl...... 165
6.1 Illustration of DNA characterization system based on the translocation of DNA through a nanopore...... 169
xx 6.2 The geometry of the nanopore and the DNA...... 171
6.3 The algorithm used to calculate DNA velocity...... 174
6.4 The inner region and outer region used for the asymptotic analysis...... 175
6.5 The geometry of the nanopore used for calculation. A 300 nm long conical silica nanopore is connected to a 50 nm long cylindrical pore. The radii of the large and small ends of the conical pore are 50nm and 5nm respectively (Storm et al., 2003)...... 178
6.6 Velocity contours in m/s calculated based on the lubrication solutions. The length of the conical pore is 300nm and the length of the cylindrical pore is 40 nm. The inlet radius of the nanopore is 50nm and the outlet radius is 5nm. The surface charge density of the nanopore is −0.06C/m2 and the electrolyte solution is 0.1 MKCL. The radius of DNA is 1 nm...... 180
∗ 6.7 The viscous drag force due to the electroosmotic flow component Feof and ∗ the pressure driven component Feof as a function of KCl concentration. The length of the conical pore is 300nm and the length of the cylindrical pore is 40 nm. The inlet radius of the nanopore is 50nm and the outlet radius is 5nm. The surface charge density of the nanopore is −0.06C/m2. 181
6.8 DNA velocity as a function of electrolyte concentration (a) and pore surface charge (b). The length of the conical pore is 300nm and the length of the cylindrical pore is 40 nm. The inlet radius of the nanopore is 50nm and the outlet radius is 5nm. The surface charge density of the nanopore is −0.1 C/m2 for (a) and the electrolyte concentration is 1M for (b)...... 182
6.9 The results for current through the nanopore as a function of time. For the numerical results, the baseline current Ibaseline is Ibaseline = 7085 pA and the amplitude is ∆I = 160 pA. The corresponding experimental data is I = 7100pA and ∆I = 140 pA(Storm et al., 2005a)...... 185
6.10 Schematic of the partial entry case. The motion of the DNA is from left to right against the bulk electroosmotic flow which is from right to left. . . . . 186
xxi 6.11 DNA velocity and ionic current through the nanopore as DNA enters the pore. The length of the conical pore is 300nm and the length of the cylin- drical pore is 40 nm. The pore is shown at the top. The inlet radius of the nanopore is 50nm and the outlet radius is 5nm. The surface charge density of the nanopore is −0.06 C/m2 and the electrolyte concentration is 1M. . 191
6.12 DNA velocity, ionic current and pressure through the nanopore as DNA enters the pore. The length of the conical pore is 300nm and the length of the cylindrical pore is 40nm. The pore is shown at the top. The inlet radius of the nanopore is 50nm and the outlet radius is 5nm. The surface charge density of the nanopore is −0.06 C/m2 and the electrolyte concentration is 1M...... 192
6.13 Velocity and ionic current change as a function of time during the entire DNA transloction process. The length of the DNA is 16.5kbps. The length of the conical pore is 300nm and the length of the cylindrical pore is 40nm. The inlet radius of the nanopore is 50nm and the outlet radius is 5nm. The surface charge density of the nanopore is −0.06 C/m2 and the electrolyte concentration is 1M...... 193
6.14 DNA velocity and ionic current through the nanopore as DNA enters a diverging nanopore. The length of the diverging nanopore is 300nm and the length of the straight cylindrical pore is 40 nm. The pore is shown at the top. The inlet radius of the nanopore is 5nm and the outlet radius is 223nm. The surface charge density of the nanopore is −0.06 C/m2 and the electrolyte concentration is 1M...... 195
6.15 DNA velocity and ionic current through the nanopore as DNA leaves a diverging nanopore. The length of the diverging nanopore is 300nm and the length of the straight cylindrical pore is 40 nm. The pore is shown at the top. The inlet radius of the nanopore is 5nm and the outlet radius is 223nm. The surface charge density of the nanopore is −0.06 C/m2 and the electrolyte concentration is 1M...... 196
6.16 Velocity and ionic current change as a function of time during the entire DNA transloction process. The length of the diverging nanopore is 300nm and the length of the straight cylindrical pore is 40 nm. The pore is shown at the top. The inlet radius of the nanopore is 5nm and the outlet radius is 223nm. The surface charge density of the nanopore is −0.06 C/m2 and the electrolyte concentration is 1M...... 197
xxii 6.17 DNA velocity, ionic current and pressure through the nanopore as DNA enters a diverging pore. The length of the conical nanopore is 300nm and the length of the straight cylindrical pore is 40 nm. The inlet radius of the nanopore is 223nm and the outlet radius is 5nm. The surface charge density of the nanopore is −0.2 C/m2 and the electrolyte concentration is 1M...... 197
6.18 The schematic for the gene delivery device...... 199
A.1 The upwind difference method and the central difference method. Here qe is the velocity of east face using upwind difference method and qe is the face velocity approximated by the central difference...... 214
A.2 An inner node on the mesh...... 215
A.3 A boundary node at the inlet...... 228
A.4 A boundary node at the outlet...... 229
A.5 A boundary node at the lower wall...... 230
A.6 The residuals of the numerical calculation...... 233
xxiii CHAPTER 1
Introduction
1.1 Background
The subject of micro- and nanofluidics deals with controlling and manipulating of fluid
flows having length scale on the order of microns and nanometers and it is a multidis- ciplinary field comprised of physics, chemistry, engineering, and biotechnology. Micro and nanofluidics has applications to drug delivery and its control, DNA and biomolecular sensing, manipulation, the manufacture of laboratories on a microchip (lab-on-a-chip) and many other areas. The new ideas developed in these areas demand a better understanding of micro/nano scale fluid flow and particle transport phenomena.
As dimensions shrink, the effective driving and dominating forces change radically.
Conventional forces resulting from pressure, inertia, viscosity or gravity that usually plays the dominant role in macroscopic flows may not be practical in micro/nanofluidic systems while forces at interfaces such as surface tension become more important due to the in- creasing ratio of interfacial area and volume (Gad-el Hak, 1999). Electrokinetic transport is a convenient method to move materials, such as water, ions and particles, in miniature systems for fast, high-resolution and low-cost analysis and synthesis. For example, a rel- atively low imposed electric field can generate significant volume flowrates which can not
1 3 Q=1µ l/min,Pressure−driven flow Q=1µ l/min, EOF 2.5
2
1.5
1
0.5
Pressure Drop(Atm) and Applied Potential(Volts) 0 10 20 30 40 50 60 70 80 Channel Height (nm)
Figure 1.1: The required pressure drop and voltage drop for nanochannels with different channel height. The required flowrate is 1µl/min; the length of the channel is 3.5 µm and the width of the channel is 2.3µm. The zeta potential of the channel wall is −11 mV and the working fluid is 1M NaCl (Conlisk, 2005).
be achieved using the conventional pressure-driven flow. As shown in Figure 1.1, the re-
quired pressure drop for a flowrate of Q =1µl/min increases from 0.006 atm to 3 atm as the height of the channel shrinks to 10 nm while the corresponding voltage drop changes
from 0.05 Voltsto 0.33 Volts. Here the volume flowrate Q ∝ h3∆p (Conlisk et al., 2002) for the pressure-driven flow in a channel with rectangular cross section and Q ∝ h∆V for the corresponding electroosmotic flow (EOF), where Q is the flowrate; ∆p and ∆V are the required pressure and voltage drop to move the fluid. It is shown that the electroosmotic
flow is more practical than the pressure-driven flow, and this makes electroosmotic flow widely used in drug delivery devices for lower requirements in cost and power.
However, fluid flow and particle transport in microfluidic devices can be difficult to control due to many complex factors such as surface composition and buffer characteris- tics. Understanding of fundamental fluid transport and particle transport is valuable for device design and optimization. Due to the difficulties of measuring fluid motion in such a
2 small scale without disturbing the flow field, theoretical analysis and numerical simulation
become essential tools.
(a) PMMA micronozzle (b) PMMA nanonozzle
Figure 1.2: Micro and nanonozzles: (a) The Polymethyl Methacrylate (PMMA) micronoz- zles used to characterize fluid flow and particle transport. (b)The Polymethyl Methacry- late (PMMA) nanonozzles fabricated using sacrificial template imprinting (STI) methods (Wang et al., 2005).
This work is supported by the National Science Foundation Nanoscale Science and En- gineering Center (NSEC) for Affordable Nanoengineering of Polymer Biomedical Devices
(CANPBD). The primary goal of NSEC is to develop polymer-based, low-cost nanoengi- neering technology and then use it to produce nanofluidic devices. Polymethyl Methacry- late (PMMA) micro/nanonozzles (shown in Figure 1.2) have been fabricated using poly- meric materials due to their biocompatibility and/or biodegradability, low cost, and recy- clability. Efforts are being made to integrate these nanonozzles and into nanofluidic de- vices for both biomedical and non biomedical applications. The theory behind transport of species in these nanonozzles and nanopores is essential to design and implement any nano/micro fluidic device. In the present work, the fluid flow and transport of ionic species
3 and charged particles in nanochannels is addressed based on theoretical analysis and nu- merical simulations. In the following section, the basic concept of electrokinetic transport is discussed.
1.2 Electrokinetic Phenomena
Figure 1.3: An illustration of the ion distribution near a glass surface and the electric double layer (EDL). The dash line denotes the approximate edge of the electric double layer.
Electrokinetic phenomena refer to electroosmosis, electrophoresis, streaming potential and sedimentation potential, which are phenomena due to the interaction of the diffuse double layer and an applied electric field generally observed in porous medium or colloidal systems. Electroosmosis is the motion of fluid flow produced by the action of an applied electric field on a fluid having a net charge. Electrophoresis is defined as the relative motion of charged particles under an electric field. Streaming potential is the potential induced by the movement of the fluid. Sedimentation potential is the potential induced by the move- ment of charged particles (Williams et al., 1995). In this thesis, fluid flow (electroosmosis) and particle motion (electrophoresis) will be investigated.
4 When an electrolyte solution characterized by charged ions comes in contact with a charged wall, a predominantly charged layer near the wall will develop. For example, as shown in Figure 1.3, when a glass surface is immersed in aqueous solution, it undergoes chemical reactions resulting in negative surface charges. These negative surface charges, and thus negative surface potentials, attract the positive ions in the solution to the wall, while the negative charges in the solution are repelled away from the wall. In immediate contact with the surface, there is a layer of cations strongly bound to the wall, called the
Stern layer. Outside the Stern layer, there is another layer where the cations are mobile called the diffuse layer. In these two layers, the surface potential significantly influences on the ion distribution and these two layers are usually termed the electric double layers
(EDL). Away from the wall, the bulk of the solution remains electrically neutral. The thickness of EDL is characterized by the Debye length λd, which is typically in the order of nanometers. Zeta potential (ζ) is widely used for quantification of the magnitude of the electrical charge at the double layer. It is defined as the electrical potential difference between the edge of the Stern layer and a location far away from the wall where the con- centration of ionic species is the so called ’bulk’ concentration. If the finite size effect of the ions is neglected, the ζ potential becomes the potential difference between the wall and the far field.
When electrodes are placed at the ends of a channel the cations are attracted to the cathode and the anions are attracted to the anode. In the bulk flow away from the walls, the concentration of cations and anions are the same and so the electric body force will balance.
But in the EDLs, they are not the same and so there is a body force acting on the flow in the EDLs. The flow near the wall will move under such a body force and thus generate a
5 Figure 1.4: Electric double layer around a charge particle.
bulk fluid movement through the channel. This flow is usually called electroosmotic flow
(EOF) or electroosmosis.
Electrophoresis is a well-known electrokinetic phenomena and it is a complement to electroosmosis. In 1807, Reuss observed that clay particles dispersed in water migrate under influence of an applied electric field (Zwolak, 1809). Electrophoresis deals with the migration of charged particles in a liquid medium under an electric field and the particles are usually rigid and non-conducting. For a charged solid particle placed in an electrolyte solution, an electric double layer of opposite charge will develop around the particle to maintain the electric neutrality of the system as shown in Figure 1.4. Under the external applied electric field, the particle will move towards the anode or the cathode due to the body force qE where q is the net charge of the particle. The electric field also induces a body force on the ions in the electrical double layer which generate viscous drag force to retard the motion of the particle. There is another force associated with the deviation of the double layer from spherical symmetry and surface conductivity due to the excess ions
6 in the diffuse layer, which is usually called the electrophoretic relaxation force (Hubbard,
2002). The viscous drag force balances all these forces which determine the velocity of
the particle. Electrophoresis has been widely used to characterize and separate colloidal
particles and macromolecules since different particles move at different electrophoretic
velocities due to different charge to radius ratio. In microfluidic systems, it is also used to
propel functionalized particles such as DNA molecules and proteins. In the present work,
the DNA electrophoretic transport through a solid state nanopore is investigated in chapter
5 and 6. Some of the basic molecular biology of DNA is discussed in the following section.
1.3 The Molecular Biology of DNA
DNA, or deoxyribonucleic acid, is a nucleic acid that contains the genetic instructions used in the development and functioning of all known living organisms and some viruses.
It carries the information that characterizes a cell and the instructions for making a new cell
(Silverstein & Nunn, 1002). DNA is a long thread-like macromolecule made of chemical units called nucleotides, each of which includes a phosphate, a sugar and a nitrogen base as shown in Figure 1.5 (a). Two polynucleotide chains are held together by hydrogen bonds and entwine like vines, in the shape of a double helix as shown in Figure 1.5 (b) forming a
DNA molecule.
There are four types of nitrogen bases found in DNA: A(adenine), T(thymine), C(cytosine) and G(guanine). Each type of base on one strand can form a bond with just one type of base on the other strand, which is called complementary base pairing. The four bases can only form two kinds of base pairing, either A-T or C-G. For example, if the sequence of one of the double strands is known as ATTCGG, the sequence of the other strand will be
TAAGCC. This characteristic is very significant and it has been applied in some of the
7 (a) DNA chemical structure
(b) DNA double helix structure
Figure 1.5: DNA double-helix structure showing base pairs as vertical lines. The diameter of the spiral of the helix is about 2nm
8 primary DNA sequencing methods (Chang, 2005). The sequence of these four bases along
the DNA backbone encodes the genetic information and it is translated to determine the
sequence of the amino acids within proteins.
DNA chain is 22 to 26 angstroms wide (2.2 to 2.6 nanometers) and one nucleotide unit is 3.4 angstroms (0.34nm) long(M et al., 1981). Based on the number of base pairs, one can easily obtain the contour length of the DNA. For example, for a 16.5kbps λ − DNA, the contour length is 0.34 × 16.5=5.6 µm. The DNA from a human cell has a contour length of 2m however it can be compacted like a ball of string in cells, the diameter of which is only about 2 µm (Calladine et al., 2004). This gives DNA the ability to not only record huge amount of information but also to be packed into a micro sized cell. DNA
+ molecules are negatively charged at physiological pH values (ph = −log10[H ], where
[H+] is the concentration of hydrogen ions in the solution) because the phosphate groups
+ in the backbone tends to donate their protons (H3O ) to the solution to reach their acid- base equilibrium. Generally DNA has a line charge of −2e per base pair which is equivalent to a surface charge density of −0.15C/m2. Electrophoresis has become a very important mechanism that biologists use to manipulate DNA.
In liquid solution, an isolated DNA chain is subjected to Brownian motion induced by the surrounding fluid molecules and behaves like a three dimensional random walk whose nodes moves randomly while keeping the connectivity of the chain. Usually for a long polymer, the persistence length (lp), which is the maximum length of the uninterrupted polymer chain persisting in a particular direction(Calladine & Drew, 1992), determines how a DNA segment should be described. For a double-stranded DNA, the persistence length is about 50 nm and for a contour length much smaller than 50 nm, the DNA be- haves like a flexible elastic rod, while for the contour length much longer than 50 nm, the
9 properties of the DNA, such as DNA diffusivity and electrophoretic mobility, can be only
described statistically (Salamone, 1996).
Figure 1.6: Conformation of DNA confined in a pore. Here h∗ is the half height of the channel of the radius of the pore; Rg is the radius of gyration of the DNA and lp is the persistence length of DNA.
A DNA molecule in a nanochannel may extend along the channel axis to a substan- tial fraction of its full contour length. The conformations of DNA in a nanochannel can be broadly classified into three groups based on the ratio of height or radius of the channel (h∗) to the DNA radius of gyration (Rg) or the DNA persistence length (lp) (Muthukumar, 2007) as shown in Figure 1.6. Here the radius of gyration (Rg) is defined as the average square distance of the chain segments from the centre of the mass of the chain and it is a statistical measure of the volume in which the DNA molecule is contained. The specific conforma- tions that DNA adopts play an important role as the DNA interacts with its environment. If
∗ the pore or channel size is larger than the DNA radius of gyration (h >Rg), the DNA is free to coil in the nanochannel and the confinement does not alter the statistical mechanical properties of DNA (shown in Figure 1.6 (a)). If the height of the channel is smaller than
10 ∗ the radius of gyration of the DNA but larger than the persistence length (lp (Bonthuis et al., 2009) as shown in Figure 1.6 (b). If the height of the channel is smaller ∗ than the persistence length of the DNA (h ∗ restricted and in the strong confinement limit (h lp), the DNA can be considered as the link of many rigid rods deflected by the boundary of the channel as shown in Figure 1.6 (c), which is usually called the Odijk regime (Odijk, 1983). In the present work, the DNA transport through a nanopore problem falls in the Odijk regime. Besides the restriction due to the nanopore, DNA is subjected to the extension due to the applied electrical driving force and the drag force, so that it translocates like a rigid rod. The details will be discussed in chapters 5 and 6. 1.4 Applications In this thesis, electroosmotic flow and the electrokinetic transport of ionic species and particles are investigated in two kinds of channels: a slit pore and a tapering nozzle with rectangular or circular cross-section. The electroosmotic flow in a slit pore with hetero- geneous surface potential is discussed in chapter 3. The electroosmotic flow and particle transport in a tapering micro/nano nozzle is discussed in chapter 4. DNA transport through a converging nanopore is discussed in chapters 5 and 6. The potential applications of these micro/nanofluidic channel/pores are discussed in this section. Heterogeneous Surface Charge 11 The heterogeneous surface charge on the microchannels has been proposed to enhance species mixing in Lab-on-a-chip devices (Erickson & Li, 2002a). The effect of heteroge- neous surface potential can also be used to reduce the band broadening. It is known that an analyte band will become distorted, or skewed, as it travels around a turn due to differences in the path length and differences in the magnitude of the electric field between the inside and the outside of the turn. Johnson et al. (2001) have introduced surface heterogeneity to the side wall of a polymeric microchannel and demonstrated how this could reduce sample band broadening around turns. Micro/nano Nozzle/Diffuser Transdermal Drug Delivery (a) (b) Figure 1.7: Images of microneedles used for transdermal drug delivery. (a) Solid micronee- dles (150µm tall) etched from a silicon wafer were used in the first study to demonstrate microneedles for transdermal delivery (Henry et al., 1998). (b) Microneedle arrays (up to 300µm tall) in standard silicon wafer using potassium hydroxide (KOH) wet etching (Wilke et al., 2005). 12 Transdermal drug delivery is an approach used to deliver drugs across the skin and into systemic circulation. Compared to other drug delivery approaches, it delivers drug directly to the blood stream, allowing the active ingredients to bypass the acid digestive environ- ment of the stomach and filtration by the liver. Moreover, it provides a large surface area and ease of accessibility for drug administration. Transdermal drug delivery systems have been in use for over 20 years and this method of delivery has become widely recognized since the introduction of nicotine patches for smoking cessation in 1991. Transdermal delivery is severely limited by the inability of the large majority of drugs to cross skin at therapeutic rates due to the great barrier imposed by skin’s outer stratum corneum layer. To increase the permeability of skin, microneedles are used. Micronee- dles are micron-dimensioned needles that can pierce the skin in a minimally invasive man- ner without causing pain. Upon piercing skin they create micro-conduits across stratum corneum and provide a direct route for transport of drugs and vaccines into the skin. Mi- croneedles can be fabricated to be long enough to penetrate the stratum corneum, but short enough not to puncture nerve endings. This reduces the chances of pain, infection, or in- jury. Usually these microneedles are made of metal or silicon and as shown in Figure 1.7, the structures of these microneedles are similar to the micronozzles shown in Figure 1.2. In a transdermal drug delivery system, drug is generally loaded in a drug reservoir integrated into the patch and transport can occur by diffusion or electromigration or elec- trophoresis called iontophoresis if an electric field is applied. In the presence of an electric field transport can be greatly improved due to the presence of electroosmosis and electro- migration. A transdermal drug delivery patch is shown in Figure 1.8 to increase antisense oligodeoxynucleotide (ODN) across the skin for cancer treatment. There is a power source 13 Figure 1.8: Schematic drawing of a transdermal drug delivery system. There is a power source in the back. The foam ring with the drug reservoir is placed under the cathode and a microprojection array is integrated into the patch(Lin et al., 2001). in the back and the drug reservoir is placed under the cathode. A microneedle array is inte- grated in the patch. Delivery of ODN is greatly improved by the microneedles and it can be controlled by duration of the delivery, donor drug concentration, current density, and active patch area(Lin et al., 2001). The theory discussed in chapter 4 can be used to characterize the drug delivery through the microneddles in this kind of devices. 14 Figure 1.9: Gene delivery using nanotip array (Fei, 2007) under an applied electric field. Gene Delivery into Cells A gene is a portion of DNA that determines a particular characteristic in an organism. Gene delivery is the process of introducing foreign DNA into host cells and it is one of the biggest challenges in the field of gene therapy. Viruses and liposomes have been widely used as carriers in animal testing, i.e. in vivo, but safety issues, such as immune response and cytotoxicity, have limited their clinical applications. Nanonozzles can be used to di- rectly transfer repaired DNA into cells without using any carrier and thus can avoid the risks associated with introducing a secondary agent. The converging channels as shown in Figure 1.2 can pre-stretch and accelerate DNA molecules before sending them to cells. Moreover, the aperture of nanotips can carry a specific dosage of genes and leave it inside cells. A short penetration of these tiny nanotips would not cause permanent damage to the cells, the idea of which is shown in Figure 1.9. DNA is usually coiled in solution and the super-coiled DNA cannot pass through channels with dimensions smaller than its gyration diameter (usually ∼ 0.1µm) in a short time. Currently the nanopores fabricated by Wang et al. (2005) (shown in Figure 1.2) have small ends of 80nm and large ends of ∼ 1µm in 15 diameter. In the proposed gene delivery device, DNA can easily enter the large end. Inside the nanonozzle, the DNA is stretched and becomes a long worm-like shape with a much smaller radial size, and will eventually migrate through the channel. The stretching of the DNA is caused by the electrical field gradient due to the tapering geometry and the shear stress acting on the DNA due to the fluid flow. The transport of DNA through the nanopores in this problem falls in the de Gennes regime as discussed in the previous section. Boimolecular Sensing and DNA Sequencing The basic idea behind biomolecular sensing using a nanopore membrane is simple and already rather old. In 1950’s, Wallace Coulter (Coulter, 1953) first introduced a method to count particles or cells in solution. The principle of this system is to place the aperture in an electrolyte solution, containing the particles only on one side, and then constantly measure the conductance of the aperture. As the particles move through the aperture either by diffusion, by an external applied electric field or by a pressure gradient, they will cause a change in conductance due to the difference in conductance between particles and the solution. Initially this system was limited to the detection of micron-sized particles due to difficulty of fabricating small apertures in membrane. In 1970’s, DeBlois & Bean (1970); DeBlois et al. (1977) extended this method to detect submicron sized particles by using a single pore of a few hundred nanometers in diameter. The nanopore used is produced by the track-etching technique and the particles to detect are as small as 90nm in diameter. Based on the same concept, Kasianowicz et al. (1996) proposed to characterize linear polymers like DNA using the protein, α-hemolysin, as the nanopore. DNA sequencing is the process of determining the nucleotide order of a given DNA fragment. It is useful in fundamental research into why and how organisms live, as well as in applied subjects such as the treatment of genetic diseases and forensic testing. Because 16 of the key nature of DNA to living things, knowledge of a DNA sequence may come in useful in practically any biological research. For example, in medicine it can be used to identify, diagnose and potentially develop treatments for genetic diseases. Due to the development in DNA sequencing, genome-based medicine has come ever closer to reality (Zwolak, 2008). Current sequencing processes are based on the chain termination method developed by Sanger et al. (1977). It be divided into four overall steps (Chan, 2005): i) DNA isolation, ii) sample preparation, iii) sequence production, and iv) assembly and analysis. It relies on very complex sample preparation and post processing of the data and the process is very expensive and time consuming. Based on current techniques, sequencing a single human genome involves costs of about 10 million USD and several months time (Fredlake et al., 2006). Although there has been some efforts to reduce the cost of the process, it is still hard to bring the overall cost lower than $10, 000 (Bayley, 2006). The NIH has set a remarkable challenge called ”the $1000 Genome” which is to sequence the complete genome of an individual human quickly and at an accessible price. Nanopore sequencing, as shown in Figure 1.10, is one possible solution to bring the cost of DNA sequencing down and it has attracted many attentions recently. The sequenc- ing method involves a membrane containing a single nano-scale pore between two halves of an electrochemical cell filled with an electrolyte solution. An external electric field is applied and the resulting ion current flowing through the electrolyte filled nanopore is recorded versus time. As the analyte, usually biomolecules and biopolymers with dimen- sions comparable to the nanopore diameter, is driven through the channel, the change in ion current is observed. The concentration of the analyte can be determined from the frequency of these current-changing translocation events and the identity of the analyte is encoded in 17 Figure 1.10: Concept for sequencing DNA by using a single protein pore. A single-stranded DNA (or RNA) molecule moves through the pore in the transmembrane electric field. As it passes a ”contact site” each base produces a characteristic modulation of the amplitude in the single channel current (Alper, 1999). the magnitude and duration of the current change (Sexton et al., 2007). Although this tech- nique is still in its early stage, the approach can already reveal limited information about base composition of DNA or RNA (Bayley & Martin, 2000). Both biological and artificial nanopores have been used for molecular sensing. Biologi- cal nanopore resistive-pulse sensors consist of a single transmembrane protein embedded in a planar lipid bilayer support. The most commonly used protein nanopore is α-hemolysin, shown in Figure 1.11 (a) and it is usually used to detect single strand DNA (ssDNA) (Bay- ley & Cremer, 2001). Synthetic nanopores (shown in Figure 1.11 (b)), usually embedded in silicon nitride and silicon oxide nanopore membranes, are used to detect double strand DNA (dsDNA) and proteins (Sexton et al., 2007; Storm et al., 2005a). For the design of the nanopore sensor, the key mechanical and electrical interactions between the DNA and the nanopore need to be well characterized (Chang & Yang, 2004). Water Purification A non-biomedical application of nanopore is the water purification. The world’s water consumption rate is doubling every 20 years, outpacing by two times the rate of population 18 (b) Sythetic nanopore (Storm et al., 2003) (a) α−hemolysin nanopore (Bayley & Cremer, 2001) Figure 1.11: Nanopores used in the experiments to analyze DNA: (a) natural α-hemolysin nanopore and (b) synthetic nanopore. growth. The supply of fresh water is on the decrease but water demand for food, industry and people is on the rise. Because of the potentially unlimited availability of seawater, peo- ple have made great efforts to try to develop feasible and cheap desalting technologies for converting salty water to fresh water. The development of membrane technology indicated a filter could be fabricated capable of separating liquids from colloidal particles or even low molecular weight dissolved solutes (ions). Surface properties of the membrane are important to the performance and the energy cost of the purification process. In chapter 3, the electrostatic interaction between a large charged molecule and the nanopore membrane is investigated to develop an artificial kidney. Extension of this work will provide founda- tional knowledge to the optimization of nanofiltration and reverse osmosis technologies for the purpose of water purification. 19 1.5 Literature Survey 1.5.1 Electroosmotic Flow When in contact with electrolytes, many solid substrates acquire a surface charge and attracts opposite ions, creating thin layers of charges next to it, called electrical double layers. As shown in Figure 1.3, under an external electric field, the fluid in the electrical double layers acquires a momentum and drags the fluid in the bulk by viscosity. The re- sulting fluid motion is called electroosmotic flow (Probstein, 2003). The earliest model of the electrical double layer is usually attributed to Helmholtz (1879), who treated the double layer based on a physical model in which a single layer of ions is adsorbed at the surface. Later Guoy (1910) and Chapman (1913) made significant improvements by introducing a diffuse model of the electrical double layer, in which the potential at a surface decreases exponentially due to adsorbed counter-ions from the solution. The current classical electri- cal double layer model was introduced by Stern (1924) by combining the Helmholtz single adsorbed layer with the Gouy-Chapman diffuse layer. Von Smoluchowski (Smoluchowski, 1903, 1912) made several contributions to our understanding of electrokinetically driven flows, especially for conditions where the electric double layer thickness is much smaller than the channel height. Later the electroosmotic flow for an infinite capillary was stud- ied by Rice & Whitehead (1965). Solutions in a narrow slit were obtained by Burgeen & Nakache (1964) in the context of the Debye-H¬uckel approximation and Levine et al. (1975) examined electroosmotic flow in a slit channel for both thin and overlapping EDLs. Con- lisk et al. (2002) solved the problem for arbitrary ionic mole fractions and the velocity and potential for strong electrolyte solutions and have considered the case where there is a po- tential difference in the direction normal to the channel walls corresponding in some cases 20 to oppositely charged walls. The Debye-H¬uckel approximation has also been investigated by Conlisk (2005). Grahame (1953) investigated the electric double layer at a plane interface for an asym- metric electrolyte and calculated electric potential and ionic concentration in the double layers. Friedl et al. (1995) investigated the electroosmotic mobility for multivalent mix- tures for the application in capillary electrophoresis and developed an empirical formula of mobility for a number of different acidic mixtures. Zheng et al. (2003) investigated the effect of multivalent ions on electroosmotic flow for multiple electrolyte components. In micro- and nanochannels having fixed surface charges, they found that adding multivalent counter-ions into a background electrolyte solution, even in very small amounts, caused significant reduction in electroosmotic flowrate in comparison with monovalent ions, while the multivalent co-ions have little effect on the electroosmotic flowrate. Datta et al. (2009) developed a site binding model to investigate the relation between the zeta potential and the divalent cations Ca+2 and Mg+2 in a background electrolyte solution. They found that the adsorption behaviors of the two divalent cations Ca+2 and Mg+2 are different, which can be explained by the stronger Ca+2 association to the silica surface, resulting in dif- ferent zeta potential when same amount of Ca+2 and Mg+2 are added in to a background electrolyte solution. The importance of the heterogeneous surface charge on electrokinetic phenomena has been recognized in theoretical studies. The electro-osmotic flow in inhomogeneously charged pores was investigated initially by Anderson & Idol (1985). They have developed an infinite-series solution for the periodically varying ζ-potential in the flow direction and found recirculation regions. The electroosmotic flow between two slabs with periodic zeta potential has been investigated later by Ajdari (1995, 1996) to generate complex flows. In 21 that work, the flow is modeled under the Debye-H¬uckel approximation and the solution under the thin Debye layer (TDL) limit is also discussed. Chang & Yang (2004) have cal- culated the case with rectangular blocks located within the microchannel with the primary focus of enhancing a suitably defined mixing efficiency. The reversed flow pattern was also observed by Stroock et al. (2000) experimentally, and they also develop a model for the flow based on the Debye-H¬uckel approximation. The potential and velocity distribution close to a step jump in potential has also been investigated analytically by Yariv (2004). It is noted that one limitation of these studies including Ajdari (1996); Erickson & Li (2002b); Chang & Yang (2004), is the assumption that the net charge density field con- forms to an equilibrium Boltzmann distribution. The Boltzmann distribution is based on the condition that the electrochemical potential must be constant everywhere at equilibrium and there is no imposed electric field and thus no bulk flow. However, axial non-uniformity in the surface potential will cause a deviation from the equilibrium that will necessitate the complete solution of Poisson-Nernst-Planck model augmented with convection of the flow field instead of a Poisson-Boltzmann equation based approach. This can be compared to corresponding classical developments in calculating the electrophoretic mobility of parti- cles under the charge polarization and relaxation effects. EOF in micro/nano channels with varying cross-sectional area has also been investi- gated for different applications. Cervera et al. (2005) have developed a model for electroos- motic flow in conical nanopores. Their model focuses on the ionic transport and the flow field is not considered. Park et al. (2006) have investigated electroosmotic flow through a suddenly constricted cylinder and they have found eddies both in the center of the chan- nel and along the perimeter due to the induced pressure gradient. Ramirez & Conlisk (2006) have examined the effects of sudden changes in channel cross-section area on the 22 EOF numerically and observed the formation of vortices or recirculation regions near the step face. Ghosal (2002) has investigated a slowly varying wall charge and slowly varying cross-section based on lubrication theory using asymptotic methods. The problem studied in the present work is similar to Ghosal’s work but he used the Helmholtz-Smoluchowski (HS) slip approximation based on the assumption of infinitely thin EDLs, which leads to an EOF with uniform streamwise flow profile. However in view of application to nano-scale nozzle/diffusers, finite thickness EDLs are considered in the present model, which leads to an EOF with a nonuniform transverse profile. It also resolves induced pressure gradients in the problem that can not be resolved by analysis based on the HS approximation. 1.5.2 DNA Transport One recent application of nanopores is to use them as detectors and even analyzers for bio-polymers (e.g. DNA). When driven through a voltage-biased nanopore, bio-polymers cause modulation of ionic current through the pore, revealing the diameter, length and con- formation of the polymers. Recently, the translocation of bio-polymers attracts increasing attention since it offers a possible solution to rapid DNA sequencing (Bayley, 2006). Both natural and synthetic nanopores are used to study the translocation process exper- imentally. Kasianowicz et al. (1996); Bayley & Cremer (2001); Meller (2001); Butler et al. (2007) experimentally investigated the single-stranded DNA (ssDNA) molecules through a natural α-hemolysin pore (shown in Figure 1.11) to characterize the DNA translocation properties by analyzing the ionic current signals. The ssDNA translocation velocity is found to be independent of DNA length if the length of the DNA is much longer than the nanopore length, and the experimental DNA translocation rate is about 0.1 ∼ 10A/µsû . 23 Compared to the natural nanopores, synthetic solid-state nanopores are more stable, flexible in pore shape and surface properties, and they are easier to be integrated into mi- cro/nanofluidic devices. To understand the relation between the DNA translocation rate and various experimental parameters such as the applied voltage bias and the concentration of the working solution, Chang et al. (2004); Fan et al. (2005); Li et al. (2003) and Storm et al. (2005a) investigated the translocation of double stranded DNA (dsDNA) molecules through a synthetic nanopore experimentally. Li et al. (2003) demonstrated the capability of observing individual molecules of dsDNA translocation, and the experimental transloca- tion rate is ∼ 0.01m/s for 3kbps and 10kbps dsDNAs. Chen et al. (2004) showed that the DNA translocation velocity increases linearly with the applied voltage drop, and it is inde- pendent of the DNA length based on the experimental data for 3kbps, 10kbps and 48.5kbps dsDNAs. On the other hand, the experimental data of Storm et al. (2005a)on11.5kbps and 48.5kbps dsDNA showed longer DNA transports slower through the nanopore, in conflict with Li et al. (2003). Modeling and simulations based on DNA dynamics have been conducted to understand the DNA translocation process, including Brownian dynamics simulations (Muthukumar, 2007; Murphy & Muthukumara, 2007; Forrey & Muthukumara, 2007), Monte Carlo sim- ulations (Muthukumar, 1999; Kantor & Kardar, 2004; Kim et al., 2004; Chen et al., 2007) and Molecular Dynamic Simulations (O’Keeffe et al., 2003; Fyta et al., 2006). These simulations generally focus on the DNA conformational change during the translocation process, ignoring the hydrodynamic interaction between the DNA and the nanopore. The DNA translocation is considered to be controlled by the entropic barriers related to the DNA conformation, and the predicted translocational velocity decreases nonlinearly with increasing DNA length. 24 Other hydrodynamic models are developed to investigate the hydrodynamic interactions locally inside the nanopore (Ashoke et al., 1997; Ghosal, 2006, 2007a,b). These models focus on the flow field and forces locally inside the nanopore, and the predicted DNA translocation velocity is independent of DNA length, which is consistent with Li et al. (2003) and Chen et al. (2004)’s experimental findings. Ashoke et al. (1997) developed an analytical model for the translocation of a slender rod in a cylindrical channel based on the Debye-H¬uckel approximation. Ghosal (2006) developed a hydrodynamic model to in- vestigate the electrically driven DNA across a converging nanopore, and showed that the proximity of the pore walls plays an important role on the magnitude of the viscous drag on DNA. In his model, the nanopore was assumed to be neutral and the Smoluchowski slip velocity was used to approximate the electroosmotic flow (EOF). Ghosal (2007a) also investigated the dependence of DNA transport on salt concentration and the predicted the force acting on a DNA inside a nanopore (Ghosal, 2007b). In these two papers, the cylin- drical pores with constant cross-section are considered and the theoretical analysis is based on the Debye-H¬uckel approximation. For the application of nanopore sequencing using a synthetic nanopore, the DNA translo- cation rates reported in the experiments (∼ 1cm/s) are generally too high to measure the ionic current modulation resulting from the difference in shape and charge of the nitrogen bases on the DNA. To slow down the DNA transport, it is very important to understand the forces that control the translocation process. It is well accepted that the DNA is elec- trophoretically driven into the nanopore but the nature of the force resisting its transport is controversial in the literature. Some of the previous work including Muthukumar (1999, 2007); Kantor & Kardar (2004); Forrey & Muthukumara (2007) claim that the uncoiling force due to the DNA conformational change is the main resisting force while Storm et al. 25 (2005b) and Fyta et al. (2006) claim that the force due to the viscous drag on the blob- like DNA configuration outside the nanopore determines the DNA translocation rate. In the present work, all these forces are investigated and compared with the electrical driving force as well as the viscous drag force acting on the DNA inside the nanopore due to the electroosmotic flow. The results have shown that the uncoiling forces are usually negligible and the drag force acting on the blob-like DNA outside the nanopore is small compared to the electrical driving force. The main resisting force comes from the viscous drag due to the electroosmotic flow inside the nanopore. Future experimental design can be focused on the optimization of the parameters that can induce larger viscous drag on the DNA for the purpose of slowing down the DNA translocation, such as the surface charge of the nanopore and the geometry of the nanopore, Figure 1.12: The geometry of the nanopore and the DNA. An ideal model for DNA translocation would involve atomic-scale simulation of poly- mer, pore and ion dynamics in the presence of an applied electric field. However such a full MD approach is computationally very expensive, particularly when the simulation must be 26 conducted over microsecond timescales to capture blockade current events (O’Keeffe et al., 2003). In this thesis, a hydrodynamic model is developed to investigate the fluid flow and ionic transport locally inside the nanopore. The nanopores considered here generally ac- quire a conical region, connected to a small pore with constant radius, as shown in Figure 1.12. This configuration is widely employed in experimental nanopores since it is easier to fabricate and produces better current signals compared to cylindrical track-etched pores (Schiedt, 2007). Most of the previous work neglects the conical region and only consid- ered the part with constant radius for simplicity. On the other hand, previous work on this problem is generally based on the Debye-H¬uckel approximation (valid for absolute value