Shear-induced migration in flowing polymer solutions: simulation of long-chain DNA in microchannels [corrected].

We simulate dilute solution dynamics of long flexible polymer molecules in pressure driven flow in channels with widths of roughly 0.1-10 times the polymer bulk radius of gyration. This is done using a self-consistent coarse-grained Langevin description of the polymer dynamics and a numerical simulation of the flow in the confined geometry that is generated by the motions of polymer segments. Results are presented for a model of DNA molecules of approximately 10-100 microm contour length in micron-scale channels. During flow, the chains migrate toward the channel centerline, in agreement with well-known experimental observations. The thickness of the resulting hydrodynamic depletion layer increases with molecular weight at constant flow strength; higher molecular weight chains therefore move with a higher average axial velocity than lower molecular weight chains. In contrast, if the hydrodynamic effects of the confining geometry are neglected, depletion of concentration is observed in the center of the channel rather than at the walls, contradicting experimental observations. The mechanisms for migration are illustrated using a simple kinetic theory dumbbell model of a confined flexible polymer. The simple theory correctly predicts the trends observed in the detailed simulations. We also examine the steady-state stretch of DNA chains as a function of channel width and flow strength. The flow strength needed to stretch a highly confined chain away from its equilibrium length is shown to increase with decreasing channel width, independent of molecular weight; this is fairly well explained using a simple blob picture.

DNA dynamics in a microchannel.

An extended Brownian dynamics simulation method is used to characterize the dynamics of long DNA molecules flowing in microchannels. The relaxation time increases due to confinement in agreement with scaling predictions. During flow the molecules migrate toward the channel center line, and thereby segregate according to molecular weight. Capturing these effects requires the detailed incorporation of solvent flow in the simulation method, demonstrating the importance of hydrodynamic effects in the dynamics of confined macromolecules.