Shear-induced migration in flowing polymer solutions: Simulation of long-chain DNA in microchannels

Cross-streamline migration and near-wall depletion of elastic fibers in micro-channel flows

The complex dynamics of elastic fibers in viscous fluids are central to many biological and industrial systems. Fluid-structure interactions underlying these dynamics govern the shape and transport of flexible fibers, and understanding these interactions can help tune flow properties in applications such as microfluidic separation, printing and clogging. In this work, we use slender-body theory to study micromechanical dynamics that arise from the coupling between the elastic backbone of a fiber and the local straining flow that contributes to filament flipping and cross-streamline migration. The resulting transverse drift is unbiased in either direction in simple shear flow. However, a non-uniform shear rate results in bias towards regions of high shear, which we connect to the shape transitions during flips. We discover a depletion layer that forms near the boundaries of pressure-driven channel flow due to the competition between such a cross-streamline drift and steric exclusion from the walls. Finally, we develop scaling laws for the curvature of filaments during flip events, demonstrating the origin of the drift bias in non-uniform flows, and confirm this behavior from our simulations. Put together, these results shed light on the role of a local and dominant coupling between elasticity and viscous resistance in dictating long-term dynamics and transport of elastic fibers in confined flows.

Brownian dynamics simulations of bottlebrush polymers in dilute solution under simple shear and uniaxial extensional flows

We study the dynamics of bottlebrush polymer molecules in dilute solutions subjected to shear and uniaxial extensional flows using Brownian dynamics simulations with hydrodynamic interaction (HI). Bottlebrush polymers are modeled using a coarse-grained representation, consisting of a set of beads interacting pairwise via a purely repulsive potential and connected by finitely extensible nonlinear springs. We present the results for molecular stretching, stress, and solution viscosity during the startup of flow as well as under steady state as a function of side chain length while keeping the backbone length fixed. In extensional flow, the backbone fractional extension and the first normal stress difference decrease with an increase in side chain length at a fixed Weissenberg number (Wi). Using simulation results both in the presence of and in the absence of HI, we show that this is primarily a consequence of steric interaction resulting from the dense grafting of side chains. In shear flow, we observe a shear-thinning behavior in all cases, although it becomes less pronounced with increasing side chain length. Furthermore, nonmonotonicity in the backbone fractional extension is observed under shear, particularly at high Wi. We contextualize our simulation results for bottlebrush polymers with respect to existing studies in the literature for linear polymers and show that the unique dynamical features characterizing bottlebrush polymers arise on account of their additional molecular thickness due to the presence of densely grafted side chains.

Predicting reaction behavior of tethered polymers in shear flow

Kinetics of force-mediated chemical reactions of end-tethered polymers with varying chain length N in varying shear rate flow γ̇ are explored via coarse-grained Brownian dynamics simulations. At fixed γ̇, force F along a polymer increases linearly with N as previously predicted; however, contrary to existing theory, the F(N) slope increases for N above a transition length that exhibits minimal dependence on γ̇. Force profiles are used in a stochastic model of a force-mediated reaction to compute the time for x percent of a polymer population to experience a reaction, tx. Observations are insensitive to the selected value of x in that tx data for varying N and γ̇ can be consistently collapsed onto a single curve via appropriate scaling, with one master curve for systems below the transition N (small N) and another for those above (large N). Different force scaling for small and large N results in orders of magnitude difference in force-mediated reaction kinetics as represented by the population response time. Data presented illustrate the possibility of designing mechano-reactive polymer populations with highly controlled response to flow across a range in γ̇.

Compression of a confined semiflexible polymer under direct and oscillating fields

The folding transition of biopolymers from the coil to compact structures has attracted wide research interest in the past and is well studied in polymer physics. Recent seminal works on DNA in confined devices have shown that these long biopolymers tend to collapse under an external field, which is contrary to the previously reported stretching of the chain. In this work, we capture the compression of a confined semiflexible polymer under direct and oscillating fields using a coarse-grained computer simulation model in the presence of long-range hydrodynamics. In the case of a semiflexible polymer chain, the inhomogeneous hydrodynamic drag from the center to the periphery of the coil couples with the chain bending to cause a swirling movement of the chain segments, leading to structural intertwining and compaction. Contrarily, a flexible chain of the same length lacks such structural deformation and forms a well-established tadpole structure. While bending rigidity profoundly influences the chain's folding favorability, we also found that subject to the direct field, chains in stronger confinements exhibit substantial compaction, contrary to the one in moderate confinements or bulk where such compaction is absent. However, an alternating field within an optimum frequency can effectuate this compression even in moderate or no confinement. This field-induced collapse is a quintessential hydrodynamic phenomenon, resulting in intertwined knotted structures even for shorter chains, unlike other spontaneous knotting experiments where it happens exclusively for longer chains.

Conformations and dynamic behaviors of confined wormlike chains in a pressure-driven flow

The conformations and dynamic behaviors of wormlike chains confined by a slit in a pressure-driven flow were investigated using dissipative particle dynamics method. The wormlike chains exhibit varying conformations due to the varying shear stresses across the slit. The wormlike chain solution can be well described by the power-law fluid, and the power-law index decreases with the increase in chain rigidity. We also presented that the wormlike chain undergoes tumbling motion in the vicinity of the wall in the presence of pressure-driven flow. We also found that the wormlike chains can migrate both away from the wall and slightly away from the slit center, and the migration away from the slit center increases as the chain rigidity is increased because of hydrodynamic interactions induced in a more rigid wormlike chain.

Two-fluid kinetic theory for dilute polymer solutions

We provide a Boltzmann-type kinetic description for dilute polymer solutions based on two-fluid theory. This Boltzmann-type description uses a quasiequilibrium based relaxation mechanism to model collisions between a polymer dumbbell and a solvent molecule. The model reproduces the desired macroscopic equations for the polymer-solvent mixture. The proposed kinetic scheme leads to a numerical algorithm which is along the lines of the lattice Boltzmann method. Finally, the algorithm is applied to describe the evolution of a perturbed Kolmogorov flow profile, whereby we recover the major elastic effect exhibited by a polymer solution, specifically, the suppression of the original inertial instability.

Dynamical and Structural Properties of Comb Long-Chain Branched Polymer in Shear Flow

Using hybrid multi-particle collision dynamics (MPCD) and a molecular dynamics (MD) method, we investigate the effect of arms and shear flow on dynamical and structural properties of the comb long-chain branched (LCB) polymer with dense arms. Firstly, we analyze dynamical properties of the LCB polymer by tracking the temporal changes on the end-to-end distance of both backbones and arms as well as the orientations of the backbone in the flow-gradient plane. Simultaneously, the rotation and tumbling behaviors with stable frequencies are observed. In other words, the LCB polymer undergoes a process of periodic stretched–folded–stretched state transition and rotation, whose period is obtained by fitting temporal changes on the orientation to a periodic function. In addition, the impact induced by random and fast motions of arms and the backbone will descend as the shear rate increases. By analyzing the period of rotation behavior of LCB polymers, we find that arms have a function in keeping the LCB polymer’s motion stable. Meanwhile, we find that the rotation period of the LCB polymer is mainly determined by the conformational distribution and the non-shrinkable state of the structure along the velocity-gradient direction. Secondly, structural properties are numerically characterized by the average gyration tensor of the LCB polymer. The changes in gyration are in accordance with the LCB polymer rolling when varying the shear rate. By analyzing the alignment of the LCB polymer and comparing with its linear and star counterparts, we find that the LCB polymer with very long arms, like the corresponding linear chain, has a high speed to reach its configuration expansion limit in the flow direction. However, the comb polymer with shorter arms has stronger resistance on configuration expansion against the imposed flow field. Moreover, with increasing arm length, the comb polymer in shear flow follows change from linear-polymer-like to capsule-like behavior.

Electroosmotic Flow Induced Lift Forces on Polymer Chains in Nanochannels

A major objective of research in nanofluidics is to achieve better selectivity in manipulating the fluxes of nano-objects and in particular of biopolymers. Numerical simulations allow one to better understand the physical mechanisms at play in such situations. We performed hybrid mesoscale simulations to investigate the properties of polymers under flows in slit pores at the nanoscale. We use multiparticle collision dynamics, an algorithm that includes hydrodynamics and thermal fluctuations, to investigate the properties of fully flexible and stiff polymers under several types of flow, showing that Poiseuille flows and electroosmotic flows can lead to quantitatively and qualitatively different behaviors of the chain. In particular, a counterintuitive phenomenon occurs in the presence of an electroosmotic flow: When the monomers are attracted by the solid surfaces through van der Waals forces, shear-induced forces lead to a stronger repulsion of the polymers from these surfaces. Such focusing of the chain in the middle of the channel increases its flowing velocity, a phenomenon that may be exploited to separate different types of polymers.

Migration and Morphology of Colloidal Gel Clusters in Cylindrical Channel Flow

We report the cluster-level structural parameters of colloidal thermogelling nanoemulsions in channel flow as a function of attractive interactions and local shear stress. The spatiotemporal evolution of the gel microstructure is obtained by directly visualizing the dispersed phase near the edge of a cylindrical channel. We observe the flow of the nanoemulsion gels in a range of radial positions (r) and shear stresses between 70 and 220 Pa, finding that the r-dependent cluster sizes are due to a balance between shear forces that yield bonds and attractive interactions that rebuild the inter-colloid bonds. In addition, the largest clusters appear to be affected by confinement and accumulate toward the central axis of the channel, resulting in a volume fraction gradient. Cluster size and volume fraction variabilities are most prominent when the attractive interactions are the strongest. Specifically, a distinct transition from sparse, fluidized clusters near the walls to concentrated, large clusters toward the center is observed. These two structural states coincide with a velocity-based transition from higher shear rates near the walls to lower shear rates toward the center of the channel. We find a compounding effect where larger gel clusters, formed under strong attractions and low shear stresses, are susceptible to shear-induced migration that intensifies r-dependent heterogeneity and deviations in the flow behavior from predictive models.

Electro-hydrodynamic extraction of DNA from mixtures of DNA and bovine serum albumin

We report separation of genomic DNA (48 kbp) from bovine serum albumin (BSA) by the electro-hydrodynamic coupling between a pressure-driven flow and a parallel electric field. Electro-hydrodynamic extraction exploits this coupling to trap DNA molecules at the entrance of a microfluidic contraction channel, while allowing proteins and salts to be flushed from the device. Samples (10 μL) containing λ-DNA (1 ng) and BSA (0.3 mg) were injected directly into the device and convected to the contraction channel entrance by a flowing buffer solution. The DNA remains trapped in this region essentially indefinitely, while proteins and salts are eluted. The effectiveness of the concept has been assessed by fluorescence measurements of DNA and BSA concentrations. Electro-hydrodynamic extraction in a single-stage device was found to enhance the concentration of DNA 40-fold, while reducing the BSA concentration by four orders of magnitude. The relative concentrations of DNA to BSA at the contraction channel entrance can be as large as 1.5 : 1, corresponding to an A260/280 ratio of 1.9. The maximum yield of DNA from a salt-free solution is 50%, while salted (150 mM) solutions have a lower yield (38%).

Flexible fibers in shear flow approach attracting periodic solutions

The three-dimensional dynamics of a single non-Brownian flexible fiber in shear flow is evaluated numerically, in the absence of inertia. A wide range of ratios A of bending to hydrodynamic forces and hundreds of initial configurations are considered. We demonstrate that flexible fibers in shear flow exhibit much more complicated evolution patterns than in the case of extensional flow, where transitions to higher-order modes of characteristic shapes are observed when A exceeds consecutive threshold values. In shear flow, we identify the existence of an attracting steady configuration and different attracting periodic motions that are approached by long-lasting rolling, tumbling, and meandering dynamical modes, respectively. We demonstrate that the final stages of the rolling and tumbling modes are effective Jeffery orbits, with the constant parameter C replaced by an exponential function that either decays or increases in time, respectively, corresponding to a systematic drift of the trajectories. In the limit of C→0, the fiber aligns with the vorticity direction and in the limit of C→∞, the fiber periodically tumbles within the shear plane. For moderate values of A, a three-dimensional meandering periodic motion exists, which corresponds to intermediate values of C. Transient, close to periodic oscillations are also detected in the early stages of the modes.

Transport of polymer-modified nanoparticles in nanochannels coated with polymers

Using molecular dynamics simulations based on explicit-solvent model, we study migration of polymer-modified nanoparticles through nanochannels coated with polymers. The polymers densely grafted on the spherical nanoparticle and the channel surface form spherical polymer brush (SPB) and planar polymer brush (PPB), respectively. The migration of the neutral polymer-modified nanoparticle is driven by electroosmotic flow (EOF). The effects of the electric field strength and the SPB–PPB interaction on polymer conformations and transport dynamics of the SPB are explored. The migration velocity of the SPB reduces as the interaction between the SPB and the PPB increases. For strong SPB–PPB interaction, the directional migration of the SPB can be triggered only after the electric field strength exceeds a critical value. The high EOF velocity forces the center of mass of the spherical nanoparticle to keep near the central region of the channel due to high shear rate close to the brush–fluid interface. Unlike electrophoresis of charged polymer-grafted spherical particles, the SPB adopts a more extended conformation in the plane perpendicular to the EOF direction.

Using molecular dynamics simulations based on explicit-solvent model, we study migration of polymer-modified nanoparticles through nanochannels coated with polymers.

Dumbbell kinetic theory for polymers in a combination of flow and external electric field

Combining Poiseuille flow with an external electric field is a demonstrated method to drive transverse migration in capillary electrophoresis. Despite both computational and experimental studies, a number of questions about how to best model polymers under these conditions remains. Attempts have been made to develop a kinetic theory for a bead-spring dumbbell model, but these have only been accurate at low electric field strength and have not captured the nonmonotonic relationship between migration and electric field strength. In this paper, we revisit the development of a kinetic theory for a bead-spring dumbbell in a combination of parabolic flow and an external electric field. The resultant theory yields a compact formula that predicts polymer concentration profiles that agree excellently with our Brownian dynamics simulations including the aforementioned nonmonotonic relationship. Furthermore, we compare our theoretical results to experimental data and find that our model nearly quantitatively predicts the position of the maximum in migration.

A mechano-reactive coarse-grained model of the blood-clotting agent von Willebrand factor

The von Willebrand Factor (vWF) is a large blood glycoprotein that aids in hemostasis. Within each vWF monomer, the A2 domain hosts a cleavage site for enzyme ADAMTS13, which regulates the size of vWF multimers. This cleavage site can only be exposed when an A2 domain unfolds, and the unfolding reaction energy landscape is highly sensitive to the force conditions on the domain. Based on previous optical tweezer experimental results, we advance here a new activated A2 monomer model (AA2MM) for coarse-grained modeling of vWF that accurately represents the force-based probabilistic change between the unfolded/refolded states. A system of springs is employed to mimic the complex mechanical response of vWF monomers subject to pulling forces. AA2MM was validated by comparing monomer scale simulation results to data from prior pulling experiments on vWF monomer fragments. The model was further validated by comparing multimer scale Brownian dynamics simulation results to experiments using microfluidic chamber microscopy to visualize tethered vWF proteins subject to flow. The A2 domain unfolding reaction was studied in bulk flow simulations (pure shear and elongation flow), giving evidence that elongational flow drives the vWF size regulation process in blood. The mechanoreactive, coarse-grained AA2MM accurately describes the complex mechanical coupling between human blood flow conditions and vWF protein reactivity.

Effects of Ionic Strength on Lateral Particle Migration in Shear-Thinning Xanthan Gum Solutions

Viscoelastic fluids, including particulate systems, are found in various biological and industrial systems including blood flow, food, cosmetics, and electronic materials. Particles suspended in viscoelastic fluids such as polymer solutions migrate laterally, forming spatially segregated streams in pressure-driven flow. Viscoelastic particle migration was recently applied to microfluidic technologies including particle counting and sorting and the micromechanical measurement of living cells. Understanding the effects on equilibrium particle positions of rheological properties of suspending viscoelastic fluid is essential for designing microfluidic applications. It has been considered that the shear-thinning behavior of viscoelastic fluid is a critical factor in determining the equilibrium particle positions. This work presents the lateral particle migration in two different xanthan gum-based viscoelastic fluids with similar shear-thinning viscosities and the linear viscoelastic properties. The flexibility and contour length of the xanthan gum molecules were tuned by varying the ionic strength of the solvent. Particles suspended in flexible and short xanthan gum solution, dissolved at high ionic strength, migrated toward the corners in a square channel, whereas particles in the rigid and long xanthan gum solutions in deionized water migrated toward the centerline. This work suggests that the structural properties of polymer molecules play significant roles in determining the equilibrium positions in shear-thinning fluids, despite similar bulk rheological properties. The current results are expected to be used in a wide range of applications such as cell counting and sorting.

Depletion layer dynamics of polyelectrolyte solutions under Poiseuille flow

The interfacial physics of complex liquids under flow remains a long-standing problem in fluid mechanics that is important for fields ranging from lubrication to nanofluidics. Liquids containing small amounts of high-molecular weight polymers are found to flow through channels faster than expected, a phenomenon attributed to the formation of boundary depletion layers that relaxes the no-slip boundary condition and allows the bulk of the fluid to slip past the walls. This work provides the most direct measurement to date of the dimension and composition of depletion layers of a polymer solution under flow. We anticipate extending this approach to help understand fluid dynamics in different regimes, such as flow in nanoconfinement and turbulence.

Complex liquids flow through channels faster than expected, an effect attributed to the formation of low-viscosity depletion layers at the boundaries. Characterization of depletion layer length scale, concentration, and dynamics has remained elusive due in large part to the lack of suitable real-space experimental techniques. The short length scales associated with depletion layers have traditionally prohibited direct imaging. By overcoming this limitation via adaptations of stimulated emission depletion (STED) microscopy, we directly measure the concentration profile of polymer solutions at a nonadsorbing wall under Poiseuille flow. Using this approach, we 1) confirm the theoretically predicted concentration profile governed by entropically driven depletion, 2) observe depletion layer narrowing at low to intermediate shear rates, and 3) report depletion layer composition that approaches pure solvent at unexpectedly low shear rates.

Shear-Induced Extensional Response Behaviors of Tethered von Willebrand Factor

We perform single-molecule flow experiments using confocal microscopy and a microfluidic device for shear rates up to 20,000 s^-1 and present results for the shear-induced unraveling and elongation of tethered von Willebrand factor (VWF) multimers. Further, we employ companion Brownian dynamics simulations to help explain details of our experimental observations using a parameterized coarse-grained model of VWF. We show that global conformational changes of tethered VWF can be accurately captured using a relatively simple mechanical model. Good agreement is found between experimental results and computational predictions for the threshold shear rate of extension, existence of nonhomogenous fluorescence distributions along unraveled multimer contours, and large variations in extensional response behaviors. Brownian dynamics simulations reveal the strong influence of varying chain length, tethering point location, and number of tethering locations on the underlying unraveling response. Through a complex molecule like VWF that naturally adopts a wide distribution of molecular size and has multiple binding sites within each molecule, this work demonstrates the power of tandem experiment and simulation for understanding flow-induced changes in biomechanical state and global conformation of macromolecules.

Modelling DNA extension and fragmentation in contractive microfluidic devices: a Brownian dynamics and computational fluid dynamics approach

Fragmenting DNA into short pieces is an essential manipulation in many biological studies, ranging from genome sequencing to molecular diagnosis. Among various DNA fragmentation methods, microfluidic hydrodynamic DNA fragmentation has huge advantages especially in terms of handling small-volume samples and being integrated into automatic and all-in-one DNA analysis equipment. Despite the fast progress in experimental studies and applications, a systematic understanding of how DNA molecules are distributed, stretched and fragmented in a confined microfluidic field is still lacking. In this work, we investigate the extension and fragmentation of DNA in a typical contractive microfluidic field, which consists of a shear flow-dominated area and an elongational flow-dominated area, using the Brownian dynamics-computational fluid dynamics method. Our results show that the shear flow at the straight part of the microfluidic channel and the elongational flow at the contractive bottleneck together determine the performance of DNA fragmentation. The average fragment size of DNA decreases with the increase of the strain rate of the elongational flow, and the upstream shear flow can significantly precondition the conformation of DNA to produce shorter and more uniform fragments. A systematic study of the dynamics of DNA fragmentation shows that DNA tends to break at the mid-point when the strain rate of elongational flow is small, and the breakage point largely deviates from the midpoint as the strain rate increases. Our simulation of the thorough DNA fragmentation process in a realistic microfluidic field agrees well with experimental results. We expect that our study can shed new light on the development of future microfluidic devices for DNA fragmentation and integrated DNA analysis devices.

The role of near-wall drag effects in the dynamics of tethered DNA under shear flow

We utilized single-molecule tethered particle motion (TPM) tracking, optimized for studying the behavior of short (0.922 μm) dsDNA molecules under shear flow conditions, in the proximity of a wall (surface). These experiments track the individual trajectories through a gold nanobead (40 nm in radius), attached to the loose end of the DNA molecules. Under such circumstances, local interactions with the wall become more pronounced, manifested through hydrodynamic interactions. To elucidate the mechanical mechanism that affects the statistics of the molecular trajectories of the tethered molecules, we estimate the resting diffusion coefficient of our system. Using this value and our measured data, we calculate the orthogonal distance of the extended DNA molecules from the surface. This calculation considers the hydrodynamic drag effect that emerges from the proximity of the molecule to the surface, using the Faxén correction factors. Our finding enables the construction of a scenario according to which the tension along the chain builds up with the applied shear force, driving the loose end of the DNA molecule away from the wall. With the extension from the wall, the characteristic times of the system decrease by three orders of magnitude, while the drag coefficients decay to a plateau value that indicates that the molecule still experiences hydrodynamic effects due to its proximity to the wall.

Steady state sedimentation of ultrasoft colloids

The structural and dynamical properties of ultra-soft colloids-star polymers-exposed to a uniform external force field are analyzed by applying the multiparticle collision dynamics technique, a hybrid coarse-grain mesoscale simulation approach, which captures thermal fluctuations and long-range hydrodynamic interactions. In the weak-field limit, the structure of the star polymer is nearly unchanged; however, in an intermediate regime, the radius of gyration decreases, in particular transverse to the sedimentation direction. In the limit of a strong field, the radius of gyration increases with field strength. Correspondingly, the sedimentation coefficient increases with increasing field strength, passes through a maximum, and decreases again at high field strengths. The maximum value depends on the functionality of the star polymer. High field strengths lead to symmetry breaking with trailing, strongly stretched polymer arms and a compact star-polymer body. In the weak-field-linear response regime, the sedimentation coefficient follows the scaling relation of a star polymer in terms of functionality and arm length.

Topology-Sensitive Microfluidic Filter for Polymers of Varying Stiffness

The separation of polymers based on their size, rigidity, and topology is an essential but also highly challenging task for nanoscience and engineering. Using hybrid molecular dynamics simulations that correctly take into account hydrodynamics, we have designed microfluidic channels for separating linear from ring polymers in dilute solutions. We establish that the transport velocity of the polymers is independent of their topology and rigidity when the channel walls are smooth and repulsive. However, when the walls are decorated with attractive spots arranged on lines parallel to the flow, ring polymers exhibit an order of magnitude higher transport velocity compared to linear chains. The spots induce a homeotropic-like reorientation of ring polymers close to walls leading to a tank treading motion along them, whereas linear chains are immobilized upon adsorption. This mechanism becomes more enhanced with increasing polymer rigidity. The presented technique holds thus promise for reliably separating nanoparticles based on their topology.

Anomalous Vascular Dynamics of Nanoworms within Blood Flow

Nanomaterials have been widely used in the design of drug delivery platforms. This work computationally explores the vascular dynamics of nanoworms as drug carriers within blood flow by considering the effects of nanoworm length, stiffness, and local physiological conditions such as hematocrit. We found that nanoworms with length of 8 μm and moderate stiffness are the optimal choice as drug carriers for circulating within normal vascular network due to their lower near wall margination. Compared to those of spherical rigid particles, these nanoworms demonstrate significant demargination behaviors at hematocrit 20%, induced by the local hydrodynamic interactions. Specifically, the interactions between nanoworms and red blood cells create asymmetrical local flow fields, resulting in the demargination of nanoworms. In addition, the flexibility of nanoworms enables them to conform to the deformed shape of red blood cells under shear flow, leading to their high concentration within the core region of vessels. Therefore, the long blood circulation time of nanoworms can be partially attributed to their demargination behaviors and intertwinement with red blood cells. According to these simulation results, tuning the length and stiffness of nanoworms is the key to design drug carries with reduced near wall margination within normal vascular networks and extend their blood circulation time.

Lateral migration of electrospun hydrogel nanofilaments in an oscillatory flow

The recent progress in bioengineering has created great interest in the dynamics and manipulation of long, deformable macromolecules interacting with fluid flow. We report experimental data on the cross-flow migration, bending, and buckling of extremely deformable hydrogel nanofilaments conveyed by an oscillatory flow into a microchannel. The changes in migration velocity and filament orientation are related to the flow velocity and the filament's initial position, deformation, and length. The observed migration dynamics of hydrogel filaments qualitatively confirms the validity of the previously developed worm-like bead-chain hydrodynamic model. The experimental data collected may help to verify the role of hydrodynamic interactions in molecular simulations of long molecular chains dynamics.

Diffusion of active dimers in a Couette flow

We study the 3D dynamics of an elastic dimer consisting of an active swimmer bound to a passive cargo, both suspended in a Couette flow. Using numerical simulations, we determine the diffusivity of such an active dimer in the presence of long-range hydrodynamic interactions for different values of its self-propulsion speed and the Couette flow. We observe that the effect of hydrodynamic interactions is greatly enhanced under the condition that self-propulsion is strong enough to contrast the shear flow. The magnitude of the effect grows with the size of the dimer's constituents relative to their distance, which makes it appreciable under experimental conditions.

Effects of Polymer Length and Salt Concentration on the Transport of ssDNA in Nanofluidic Channels

Electrokinetic phenomena in micro/nanofluidic channels have attracted considerable attention because precise control of molecular transport in liquids is required to optically and electrically capture the behavior of single molecules. However, the detailed mechanisms of polymer transport influenced by electroosmotic flows and electric fields in micro/nanofluidic channels have not yet been elucidated. In this study, a Langevin dynamics simulation was used to investigate the electrokinetic transport of single-stranded DNA (ssDNA) in a cylindrical nanochannel, employing a coarse-grained bead-spring model that quantitatively reproduced the radius of gyration, diffusion coefficient, and electrophoretic mobility of the polymer. Using this practical scale model, transport regimes of ssDNA with respect to the ζ-potential of the channel wall, the ion concentration, and the polymer length were successfully characterized. It was found that the relationship between the radius of gyration of ssDNA and the channel radius is critical to the formation of deformation regimes in a narrow channel. We conclude that a combination of electroosmotic flow velocity gradients and electric fields due to electrically polarized channel surfaces affects the alignment of molecular conformations, such that the ssDNA is stretched/compressed at negative/positive ζ-potentials in comparatively low-concentration solutions. Furthermore, this work suggests the possibility of controlling the center-of-mass position by tuning the salt concentration. These results should be applicable to the design of molecular manipulation techniques based on liquid flows in micro/nanofluidic devices.

Dynamics of flexible fibers and vesicles in Poiseuille flow at low Reynolds number

The dynamics of flexible fibers and vesicles in unbounded planar Poiseuille flow at low Reynolds number is shown to exhibit similar basic features, when their equilibrium (moderate) aspect ratio is the same and vesicle viscosity contrast is relatively high. Tumbling, lateral migration, accumulation and shape evolution of these two types of flexible objects are analyzed numerically. The linear dependence of the accumulation position on relative bending rigidity, and other universal scalings are derived from the local shear flow approximation.

Single Molecule Hydrodynamic Separation Allows Sensitive and Quantitative Analysis of DNA Conformation and Binding Interactions in Free Solution

Limited tools exist that are capable of monitoring nucleic acid conformations, fluctuations, and distributions in free solution environments. Single molecule free solution hydrodynamic separation enables the unique ability to quantitatively analyze nucleic acid biophysics in free solution. Single molecule fluorescent burst data and separation chromatograms can give layered insight into global DNA conformation, binding interactions, and molecular distributions. First, we show that global conformation of individual DNA molecules can be directly visualized by examining single molecule fluorescent burst shapes and that DNA exists in a dynamic equilibrium of fluctuating conformations as it is driven by Poiseuille flow through micron-sized channels. We then show that this dynamic equilibrium of DNA conformations is reflected as shifts in hydrodynamic mobility that can be perturbed using salt and ionic strength to affect packing density. Next, we demonstrate that these shifts in hydrodynamic mobility can be used to investigate hybridization thermodynamics and binding interactions. We distinguish and classify multiple interactions within a single sample, and demonstrate quantification amidst large concentration differences for the detection of rare species. Finally, we demonstrate that these differences can resolve perfect complement, 2 bp mismatched, and 3 bp mismatched sequences. Such a system can be used to garner diverse information about DNA conformation and structure, and potentially be extended to other molecules and mixed-species interactions, such as between nucleic acids and proteins or synthetic polymers.

Nonlinear electro-osmosis of dilute non-adsorbing polymer solutions with low ionic strength

Nonlinear electro-osmotic behaviour of dilute non-adsorbing polymer solutions with low salinity is investigated using Brownian dynamics simulations and a kinetic theory. In the Brownian simulations, the hydrodynamic interaction between the polymers and a no-slip wall is considered using the Rotne-Prager approximation of the Blake tensor. In a plug flow under a sufficiently strong applied electric field, the polymer migrates toward the bulk, forming a depletion layer thicker than the equilibrium one. Consequently, the electro-osmotic mobility increases nonlinearly with increasing electric field and becomes saturated. This nonlinear mobility does not depend qualitatively on the details of the rheological properties of the polymer solution. Analytical calculations using the kinetic theory for the same system quantitatively reproduce the results of the Brownian dynamics simulation well.

Dynamics of flexible fibers in shear flow

Dynamics of flexible non-Brownian fibers in shear flow at low-Reynolds-number are analyzed numerically for a wide range of the ratios A of the fiber bending force to the viscous drag force. Initially, the fibers are aligned with the flow, and later they move in the plane perpendicular to the flow vorticity. A surprisingly rich spectrum of different modes is observed when the value of A is systematically changed, with sharp transitions between coiled and straightening out modes, period-doubling bifurcations from periodic to migrating solutions, irregular dynamics, and chaos.

Mesoscopic simulation of single DNA dynamics in rotational flows

In this numerical study, the transport and dynamics of an isolated DNA in rotational flow generated in a microchannel have been investigated using dissipative particle dynamics. Often, inertial flow through microchannels with a sudden change in surface structure facilitates a re-circulation or vortex region. The conformation and mobility of the bio-polymer under the influence of such rotating fluid inside a square cavity of the microchannel is analyzed. The flexible polymer chain is found to migrate towards the rotating region and follows the vortex streamline. The orientation, size and tumbling period of polymer strands are affected by the strength of the microvortex. At elevated flow rates, the macromolecule prefers to remain inside the vortex and a hydrodynamic trap is formed. Moreover, residence time of the single molecule in the microcavity is significantly influenced by the chain length and flow strength. Further, it has been demonstrated that, such entrapment duration can be strategically altered by modifying the hydrophobicity of the microchannel.

Visualization of large elongated DNA molecules

Long and linear DNA molecules are the mainstream single-molecule analytes for a variety of biochemical analysis within microfluidic devices, including functionalized surfaces and nanostructures. However, for biochemical analysis, large DNA molecules have to be unraveled, elongated, and visualized to obtain biochemical and genomic information. To date, elongated DNA molecules have been exploited in the development of a number of genome analysis systems as well as for the study of polymer physics due to the advantage of direct visualization of single DNA molecule. Moreover, each single DNA molecule provides individual information, which makes it useful for stochastic event analysis. Therefore, numerous studies of enzymatic random motions have been performed on a large elongated DNA molecule. In this review, we introduce mechanisms to elongate DNA molecules using microfluidics and nanostructures in the beginning. Secondly, we discuss how elongated DNA molecules have been utilized to obtain biochemical and genomic information by direct visualization of DNA molecules. Finally, we reviewed the approaches used to study the interaction of proteins and large DNA molecules. Although DNA-protein interactions have been investigated for many decades, it is noticeable that there have been significant achievements for the last five years. Therefore, we focus mainly on recent developments for monitoring enzymatic activity on large elongated DNA molecules.