Dynamics of a driven confined polyelectrolyte solution Debarshee Bagchi1, a) and Monica Olvera de la Cruz1, 2, b) 1)Department of Materials Science and Engineering, Northwestern University, Evanston, IL 60208, United States 2)Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208, United States (Dated: September 15, 2020) The transport of polyelectrolytes confined by oppositely charged surfaces and driven by a constant electric field is of interest in studies of DNA separation according to size. Using molecular dynamics simulations that include surface po- larization effect, we find that the mobilities of the polyelectrolytes and their counterions change non-monotonically with the confinement surface charge density. For an optimum value of the confinement charge density, efficient separation of polyelectrolytes can be achieved over a wide range of polyelectrolyte charge due to the differential friction imparted by the oppositely charged confinement on the polyelectrolyte chains. Furthermore, by altering the placement of the charged confinement counterions, enhanced polyelectrolyte separation can be achieved by utilizing surface polarization effect due to dielectric mismatch between the media inside and outside the confinement. I. INTRODUCTION out sieving matrices have been proposed, such as by using nanopattered structures, but these techniques are still in their Electrophoretic techniques for separating and sequencing infancy and have low separation speed and resolution. As DNA, proteins, other biomolecules and synthetic polymers ac- such, newer and more efficient techniques for separating DNA cording to size or charge, have been extensively used in recent and other bio-macromolecules are being continuously pursued times for genome analysis1,2, clinical diagnostics3, forensic (for reviews on theoretical, computational, and experimental investigations4, industrial processes5 and many other applica- works on conventional and more recent electrophoretic tech- 21–28 tions. After it was realized that the mobility of DNA in free so- niques, see ). lution is independent of its length6,7, besides the conventional In this article, we show that a separation technique without gel electrophoresis technique8, which is insensitive to sepa- involving sieving matrices or pulsed electric fields, and read- rate long molecular weight polyelectrolytes9,10, several other ily integrable with nanofluidic systems, can be conceptualized separation techniques have been subsequently developed, using a charged confinement. To optimize the resolution of such as the capillary electrophoresis11,12, pulsed-field gel the separation technique, we study the dynamics of negatively electrophoresis13, pulsed-field gradient gel electrophoresis14, charged polymer solutions under confinement by a positively pulsed-field capillary gel electrophoresis15, microchip capil- charged cylindrical nanochannel, and driven by a constant ex- lary electrophoresis16, dielectrophoresis17, entropic traps18, ternal electric field. We also examine the effect of surface nanopatterned surface19 and flat surface electrophoresis20, to polarization due to mismatch of dielectric constants between name a few. Each of these techniques has their advantages the media inside and outside the charged confinement. Pre- and shortcomings, based on factors such as separation reso- vious studies have shown that polarization effect often leads lution, speed, reproducibility, cost effectiveness, sample re- to quantitative as well as qualitative changes in the physi- quirements, automation, ease of integration etc. For exam- cal properties of charged systems29,30. Here, we demonstrate ple, pulsed-field gel electrophoresis, the most widely used when and why the properties of the confined polyelectrolyte method for large DNA molecules, suffers from problems of solution are affected by dielectric mismatch, and use the re- slow speed and difficulty of automation. On the other hand, sults of our analysis to optimize the separation technique. It is conventional capillary electrophoresis has the advantages of clear that the molecular weight dependent friction of the neg- being fast and efficient with low sample requirements and is atively charged polyelectrolytes, adsorbed on to an oppositely easily integrable, but on the downside, separation sizes are charged surface, plays a key role in separating the polyelec- limited to the kbp range for DNA electrophoresis, along with trolyte chains, a critical step for reading the genome of an the difficulty of loading a high viscosity polymer solution in- organism. side a narrow capillary. The problem of loading a high vis- Motivated by the generation of strong electrostatic corre- arXiv:2009.05656v1 [cond-mat.soft] 11 Sep 2020 cosity polymer solution inside a narrow capillary can be cir- lations when charges are present in a low dielectric environ- cumvented by using an ultradilute polymer solution, as in the ment, such as when charged particles dissolve in organic sol- transient entanglement coupling mechanism11, but then the vents in contact with aqueous solution31, or when polyelec- resolution for short DNA fragments becomes poor. Moreover, trolytes are confined in cavities with ions in the low dielec- for these methods, the separation resolution is restricted by tric media29, we consider two models, referred to as Model-I the properties of the sieving matrix. Recently, methods with- and Model-II. In Model-I, the counterions of the positively charged surface are in the high dielectric (aqueous) medium and in Model-II, the counterions are in the low dielectric medium (this happens when, say, intercalated ions in graphene a)E-mail address:[email protected] are in contact with an aqueous solution or when the charged b)E-mail address: [email protected] particles in oil release their protons at the interface with water 2 (see Fig. 2 in31). ing the software package Large-scale Atomic/Molecular Mas- Using coarse-grained molecular simulations, we demon- sively Parallel Simulator (LAMMPS)33, and for computing strate that the confined polyelectrolyte solution has rich trans- the induced charges due to dielectric mismatch, we employ port properties that can be conveniently controlled by tuning an efficient algorithm known as the Induced Charge Compu- the surface charge density (SCD) on the confinement in both tation (ICC) using the boundary element method (see34 and Model-I and Model-II. In particular, when polarization effect references therein for more details). is taken into account in Model-II, polyelectrolytes can be sep- arated with high resolution due to image charge effects and strong ionic correlations, if the SCD is of the opposite sign to that of the polyelectrolyte charge. Since the key factor for separation is the differential friction experienced by polyelec- trolyte chains of different charges (or equivalently, different Figure 1. Model-I: Typical configuration of polyelectrolyte chains sizes, as in our case) due to the charged confinement, we con- (blue), counterions (red), and surface counterions (green) inside a centrate on the case for which the confinement surface and cylindrical confinement. the polyelectrolyte are oppositely charged (although for com- pleteness, in the SI we present data of the mobility when the To characterize the dynamics of the confined polyelec- confinement surface charge is of the same sign as that of the trolyte solution, the main quantity of interest is the mobil- polyelectrolyte). Besides electrophoresis, our study is also vz relevant for understanding confinement effects on the dynam- ity m = , where vz is the average velocity of the charged qEz ics of charged soft-matter systems and in the design of smart species parallel to the direction of the applied external field. 32 nanofluidic devices . We denote the mobility of the polyelectrolyte chains and their counterions as mP and mC respectively. The mobility of the polyelectrolyte chains depends on the complex interplay of II. THE CONFINED POLYELECTROLYTE MODEL (a) the applied electric field ~E, (b) the charge on each chain −Nme, (c) the charge density on the confinement S, and (d) For our course-grained molecular dynamics simulation, the the counterion screening (cs), which we define as the fraction polyelectrolyte is modeled as a linear bead-spring chain with fcs of the total counterions that are present within a distance explicit counterions in an implicit water-like solvent. This of 1:5s from any monomer of a polyelectrolyte chain. The polyelectrolyte solution is confined inside an impenetrable counterion screening fraction fcs is an important quantity for cylindrical confinement that is comprised of discrete Lennard- understanding the transport features of the system since it is Jones (LJ) beads. Each confinement surface bead carry a related to the effective charge of the polyelectrolyte chains. partial charge and the confinement has a charge density de- noted by S. All the beads on each polyelectrolyte chain are negatively charged with one unit of elementary charge III. SIMULATION RESULTS FOR MODEL-I −e. There are typically Nc = 10 polyelectrolyte chains in- A. Mobilities without Polarization effect side the confinement, each with Nm = 40 monomers, for most of our simulations. The consecutive beads in the polyelec- trolyte chain are connected by finite extensible nonlinear elas- We begin our analysis with Model-I, for which we first tic (FENE) bonds. For overall electroneutrality, there are benchmark the dynamics of the confined polyelectrolyte solu- Ncions = NcNm positively