Fast DNA translocation through a solid-state nanopore

We report translocation experiments on double-strand DNA through a silicon oxide nanopore. Samples containing DNA fragments with seven different lengths between 2000 to 96000 basepairs have been electrophoretically driven through a 10 nm pore. We find a power-law scaling of the translocation time versus length, with an exponent of 1.26 $\pm$ 0.07. This behavior is qualitatively different from the linear behavior observed in similar experiments performed with protein pores. We address the observed nonlinear scaling in a theoretical model that describes experiments where hydrodynamic drag on the section of the polymer outside the pore is the dominant force counteracting the driving. We show that this is the case in our experiments and derive a power-law scaling with an exponent of 1.18, in excellent agreement with our data.Comment: 5 pages, 2 figures. Submitted to PR

Forced Imbibition - a Tool for Determining Laplace Pressure, Drag Force and Slip Length in Capillary Filling Experiments

When a very thin capillary is inserted into a liquid, the liquid is sucked into it: this imbibition process is controlled by a balance of capillary and drag forces, which are hard to quantify experimentally, in particularly considering flow on the nanoscale. By computer experiments using a generic coarse-grained model, it is shown that an analysis of imbibition forced by a controllable external pressure quantifies relevant physical parameter such as the Laplace pressure, Darcy's permeability, effective pore radius, effective viscosity, dynamic contact angle and slip length of the fluid flowing into the pore. In determining all these parameters independently, the consistency of our analysis of such forced imbibition processes is demonstrated.Comment: 4 pages, 5 figure

Does Young's equation hold on the nanoscale? A Monte Carlo test for the binary Lennard-Jones fluid

When a phase-separated binary ($A+B$) mixture is exposed to a wall, that preferentially attracts one of the components, interfaces between A-rich and B-rich domains in general meet the wall making a contact angle $\theta$. Young's equation describes this angle in terms of a balance between the $A-B$ interfacial tension $\gamma_{AB}$ and the surface tensions $\gamma_{wA}$, $\gamma_{wB}$ between, respectively, the $A$- and $B$-rich phases and the wall, $\gamma _{AB} \cos \theta =\gamma_{wA}-\gamma_{wB}$. By Monte Carlo simulations of bridges, formed by one of the components in a binary Lennard-Jones liquid, connecting the two walls of a nanoscopic slit pore, $\theta$ is estimated from the inclination of the interfaces, as a function of the wall-fluid interaction strength. The information on the surface tensions $\gamma_{wA}$, $\gamma_{wB}$ are obtained independently from a new thermodynamic integration method, while $\gamma_{AB}$ is found from the finite-size scaling analysis of the concentration distribution function. We show that Young's equation describes the contact angles of the actual nanoscale interfaces for this model rather accurately and location of the (first order) wetting transition is estimated.Comment: 6 pages, 6 figure

Anisotropic coarse-grained statistical potentials improve the ability to identify native-like protein structures

We present a new method to extract distance and orientation dependent potentials between amino acid side chains using a database of protein structures and the standard Boltzmann device. The importance of orientation dependent interactions is first established by computing orientational order parameters for proteins with alpha-helical and beta-sheet architecture. Extraction of the anisotropic interactions requires defining local reference frames for each amino acid that uniquely determine the coordinates of the neighboring residues. Using the local reference frames and histograms of the radial and angular correlation functions for a standard set of non-homologue protein structures, we construct the anisotropic pair potentials. The performance of the orientation dependent potentials was studied using a large database of decoy proteins. The results demonstrate that the new distance and orientation dependent residue-residue potentials present a significantly improved ability to recognize native folds from a set of native and decoy protein structures.Comment: Submitted to "The Journal of Chemical Physics

Chaperone-assisted translocation of a polymer through a nanopore

Using Langevin dynamics simulations, we investigate the dynamics of chaperone-assisted translocation of a flexible polymer through a nanopore. We find that increasing the binding energy $\epsilon$ between the chaperone and the chain and the chaperone concentration $N_c$ can greatly improve the translocation probability. Particularly, with increasing the chaperone concentration a maximum translocation probability is observed for weak binding. For a fixed chaperone concentration, the histogram of translocation time $\tau$ has a transition from long-tailed distribution to Gaussian distribution with increasing $\epsilon$. $\tau$ rapidly decreases and then almost saturates with increasing binding energy for short chain, however, it has a minimum for longer chains at lower chaperone concentration. We also show that $\tau$ has a minimum as a function of the chaperone concentration. For different $\epsilon$, a nonuniversal dependence of $\tau$ on the chain length $N$ is also observed. These results can be interpreted by characteristic entropic effects for flexible polymers induced by either crowding effect from high chaperone concentration or the intersegmental binding for the high binding energy.Comment: 10 pages, to appear in J. Am. Chem. So

Driven polymer translocation through a nanopore: a manifestation of anomalous diffusion

We study the translocation dynamics of a polymer chain threaded through a nanopore by an external force. By means of diverse methods (scaling arguments, fractional calculus and Monte Carlo simulation) we show that the relevant dynamic variable, the translocated number of segments $s(t)$, displays an {\em anomalous} diffusive behavior even in the {\em presence} of an external force. The anomalous dynamics of the translocation process is governed by the same universal exponent $\alpha = 2/(2\nu +2 - \gamma_1)$, where $\nu$ is the Flory exponent and $\gamma_1$ - the surface exponent, which was established recently for the case of non-driven polymer chain threading through a nanopore. A closed analytic expression for the probability distribution function $W(s, t)$, which follows from the relevant {\em fractional} Fokker - Planck equation, is derived in terms of the polymer chain length $N$ and the applied drag force $f$. It is found that the average translocation time scales as $\tau \propto f^{-1}N^{\frac{2}{\alpha} -1}$. Also the corresponding time dependent statistical moments, $\propto t^{\alpha}$ and $\propto t^{2\alpha}$ reveal unambiguously the anomalous nature of the translocation dynamics and permit direct measurement of $\alpha$ in experiments. These findings are tested and found to be in perfect agreement with extensive Monte Carlo (MC) simulations.Comment: 6 pages, 4 figures, accepted to Europhys. Lett; some references were supplemented; typos were correcte

Representations for Three-Body T-Matrix on Unphysical Sheets: Proofs

A proof is given for the explicit representations which have been formulated in the author's previous work (nucl-th/9505028) for the Faddeev components of three-body T-matrix continued analytically on unphysical sheets of the energy Riemann surface. Also, the analogous representations for analytical continuation of the three-body scattering matrices and resolvent are proved. An algorithm to search for the three-body resonances on the base of the Faddeev differential equations is discussed.Comment: 98 Kb; LaTeX; Journal-ref was added (the title changed in the journal

Translocation of structured polynucleotides through nanopores

We investigate theoretically the translocation of structured RNA/DNA molecules through narrow pores which allow single but not double strands to pass. The unzipping of basepaired regions within the molecules presents significant kinetic barriers for the translocation process. We show that this circumstance may be exploited to determine the full basepairing pattern of polynucleotides, including RNA pseudoknots. The crucial requirement is that the translocation dynamics (i.e., the length of the translocated molecular segment) needs to be recorded as a function of time with a spatial resolution of a few nucleotides. This could be achieved, for instance, by applying a mechanical driving force for translocation and recording force-extension curves (FEC's) with a device such as an atomic force microscope or optical tweezers. Our analysis suggests that with this added spatial resolution, nanopores could be transformed into a powerful experimental tool to study the folding of nucleic acids.Comment: 9 pages, 5 figure

Anomalous Dynamics of Translocation

We study the dynamics of the passage of a polymer through a membrane pore (translocation), focusing on the scaling properties with the number of monomers $N$. The natural coordinate for translocation is the number of monomers on one side of the hole at a given time. Commonly used models which assume Brownian dynamics for this variable predict a mean (unforced) passage time $\tau$ that scales as $N^2$, even in the presence of an entropic barrier. However, the time it takes for a free polymer to diffuse a distance of the order of its radius by Rouse dynamics scales with an exponent larger than 2, and this should provide a lower bound to the translocation time. To resolve this discrepancy, we perform numerical simulations with Rouse dynamics for both phantom (in space dimensions $d=1$ and 2), and self-avoiding (in $d=2$) chains. The results indicate that for large $N$, translocation times scale in the same manner as diffusion times, but with a larger prefactor that depends on the size of the hole. Such scaling implies anomalous dynamics for the translocation process. In particular, the fluctuations in the monomer number at the hole are predicted to be non-diffusive at short times, while the average pulling velocity of the polymer in the presence of a chemical potential difference is predicted to depend on $N$.Comment: 9 pages, 9 figures. Submitted to Physical Review

Capillary filling with wall corrugations] Capillary filling in microchannels with wall corrugations: A comparative study of the Concus-Finn criterion by continuum, kinetic and atomistic approaches

We study the impact of wall corrugations in microchannels on the process of capillary filling by means of three broadly used methods - Computational Fluid Dynamics (CFD), Lattice-Boltzmann Equations (LBE) and Molecular Dynamics (MD). The numerical results of these approaches are compared and tested against the Concus-Finn (CF) criterion, which predicts pinning of the contact line at rectangular ridges perpendicular to flow for contact angles theta > 45. While for theta = 30, theta = 40 (no flow) and theta = 60 (flow) all methods are found to produce data consistent with the CF criterion, at theta = 50 the numerical experiments provide different results. Whilst pinning of the liquid front is observed both in the LB and CFD simulations, MD simulations show that molecular fluctuations allow front propagation even above the critical value predicted by the deterministic CF criterion, thereby introducing a sensitivity to the obstacle heigth.Comment: 25 pages, 8 figures, Langmuir in pres