Dynamics of driven translocation of semiflexible polymers.

We study translocation of semiflexible polymers driven by force f_{d} inside a nanometer-scale pore using our three-dimensional Langevin dynamics model. We show that the translocation time τ increases with increasing bending rigidity κ. Similarly, the exponent ß for the scaling of τ with polymer length N,τ∼N^{ß}, increases with increasing κ as well as with increasing f_{d}. By comparing waiting times between semiflexible and fully flexible polymers we show that for realistic f_{d} translocation dynamics is to a large extent, but not completely, determined by the polymer's elastic length measured in number of Kuhn segments N_{Kuhn}. Unlike in driven translocation of flexible polymers, friction related to the polymer segment on the trans side has a considerable effect on the resulting dynamics. This friction is intermittently reduced by buckling of the polymer segment in the vicinity of the pore opening on the trans side. We show that in the experimentally relevant regime where viscosity is higher than in computer simulation models, the probability for this buckling increases with increasing f_{d}, giving rise to a larger contribution to the trans side friction at small f_{d}. Similarly to flexible polymers, we find significant center-of-mass diffusion of the cis side polymer segment which speeds up translocation. This effect is larger for smaller f_{d}. However, this speedup is smaller than the slowing down due to the trans side friction. At large enough N_{Kuhn}, the roles can be seen to be reversed, and the dynamics of flexible polymers can be reached. However, for example, polymers used in translocation experiments of DNA are elastically so short that the finite-length dynamics outlined here applies.

Uniform description of polymer ejection dynamics from capsid with and without hydrodynamics.

We use stochastic rotation dynamics (SRD) to examine the dynamics of the ejection of an initially strongly confined flexible polymer from a spherical capsid with and without hydrodynamics. The results obtained using stochastic rotation dynamics (SRD) are compared to similar Langevin simulations. Inclusion of hydrodynamic modes speeds up the ejection but also allows the part of the polymer outside the capsid to expand closer to equilibrium. This shows as higher values of radius of gyration when hydrodynamics are enabled. By examining the waiting times of individual polymer beads, we find that the waiting time t_{w} grows with the number of ejected monomers s as a sum of two exponents. When ≈63% of the polymer has ejected, the ejection enters the regime of slower dynamics. The functional form of t_{w} versus s is universal for all ejection processes starting from the same initial monomer densities. Inclusion of hydrodynamics only reduces its magnitude. Consequently, we define a universal scaling function h such that the cumulative waiting time t=N_{0}h(s/N_{0}) for large N_{0}. Our unprecedentedly precise measurements of force indicate that this form for t_{w}(s) originates from the corresponding force toward the pore decreasing superexponentially at the end of the ejection. Our measured t_{w}(s) explains the apparent superlinear scaling of the ejection time with the polymer length for short polymers. However, for asymptotically long polymers, t_{w}(s) predicts linear scaling.

Quantification of tension to explain bias dependence of driven polymer translocation dynamics.

Motivated by identifying the origin of the bias dependence of tension propagation, we investigate methods for measuring tension propagation quantitatively in computer simulations of driven polymer translocation. Here, the motion of flexible polymer chains through a narrow pore is simulated using Langevin dynamics. We measure tension forces, bead velocities, bead distances, and bond angles along the polymer at all stages of translocation with unprecedented precision. Measurements are done at a standard temperature used in simulations and at zero temperature to pin down the effect of fluctuations. The measured quantities were found to give qualitatively similar characteristics, but the bias dependence could be determined only using tension force. We find that in the scaling relation τ∼N^{ß}f_{d}^{α} for translocation time τ, the polymer length N, and the bias force f_{d}, the increase of the exponent ß with bias is caused by center-of-mass diffusion of the polymer toward the pore on the cis side. We find that this diffusion also causes the exponent α to deviate from the ideal value -1. The bias dependence of ß was found to result from combination of diffusion and pore friction and so be relevant for polymers that are too short to be considered asymptotically long. The effect is relevant in experiments all of which are made using polymers whose lengths are far below the asymptotic limit. Thereby, our results also corroborate the theoretical prediction by Sakaue's theory [Polymers 8, 424 (2016)2073-436010.3390/polym8120424] that there should not be bias dependence of ß for asymptotically long polymers. By excluding fluctuations we also show that monomer crowding at the pore exit cannot have a measurable effect on translocation dynamics under realistic conditions.

Rigidity-induced scale invariance in polymer ejection from capsid.

While the dynamics of a fully flexible polymer ejecting a capsid through a nanopore has been extensively studied, the ejection dynamics of semiflexible polymers has not been properly characterized. Here we report results from simulations of ejection dynamics of semiflexible polymers ejecting from spherical capsids. Ejections start from strongly confined polymer conformations of constant initial monomer density. We find that, unlike for fully flexible polymers, for semiflexible polymers the force measured at the pore does not show a direct relation to the instantaneous ejection velocity. The cumulative waiting time t(s), that is, the time at which a monomer s exits the capsid the last time, shows a clear change when increasing the polymer rigidity κ. The major part of an ejecting polymer is driven out of the capsid by internal pressure. At the final stage the polymer escapes the capsid by diffusion. For the driven part there is a crossover from essentially exponential growth of t with s of the fully flexible polymers to a scale-invariant form. In addition, a clear dependence of t on polymer length N_{0} was found. These findings combined give the dependence t(s)∝N_{0}^{0.55}s^{1.33} for the strongly rigid polymers. This crossover in dynamics where κ acts as a control parameter is reminiscent of a phase transition. This analogy is further enhanced by our finding a perfect data collapse of t for polymers of different N_{0} and any constant κ.

Driven polymer translocation in good and bad solvent: Effects of hydrodynamics and tension propagation.

We investigate the driven polymer translocation through a nanometer-scale pore in the presence and absence of hydrodynamics both in good and bad solvent. We present our results on tension propagating along the polymer segment on the cis side that is measured for the first time using our method that works also in the presence of hydrodynamics. For simulations we use stochastic rotation dynamics, also called multiparticle collision dynamics. We find that in the good solvent the tension propagates very similarly whether hydrodynamics is included or not. Only the tensed segment is by a constant factor shorter in the presence of hydrodynamics. The shorter tensed segment and the hydrodynamic interactions contribute to a smaller friction for the translocating polymer when hydrodynamics is included, which shows as smaller waiting times and a smaller exponent in the scaling of the translocation time with the polymer length. In the bad solvent hydrodynamics has a minimal effect on polymer translocation, in contrast to the good solvent, where it speeds up translocation. We find that under bad-solvent conditions tension does not spread appreciably along the polymer. Consequently, translocation time does not scale with the polymer length. By measuring the effective friction in a setup where a polymer in free solvent is pulled by a constant force at the end, we find that hydrodynamics does speed up collective polymer motion in the bad solvent even more effectively than in the good solvent. However, hydrodynamics has a negligible effect on the motion of individual monomers within the highly correlated globular conformation on the cis side and hence on the entire driven translocation under bad-solvent conditions.

Chaperone-assisted translocation of flexible polymers in three dimensions.

Polymer translocation through a nanometer-scale pore assisted by chaperones binding to the polymer is a process encountered in vivo for proteins. Studying the relevant models by computer simulations is computationally demanding. Accordingly, previous studies are either for stiff polymers in three dimensions or flexible polymers in two dimensions. Here, we study chaperone-assisted translocation of flexible polymers in three dimensions using Langevin dynamics. We show that differences in binding mechanisms, more specifically, whether a chaperone can bind to a single site or multiple sites on the polymer, lead to substantial differences in translocation dynamics in three dimensions. We show that the single-binding mode leads to dynamics that is very much like that in the constant-force driven translocation and accordingly mainly determined by tension propagation on the cis side. We obtain ß≈1.26 for the exponent for the scaling of the translocation time with polymer length. This fairly low value can be explained by the additional friction due to binding particles. The multiple-site binding leads to translocation the dynamics of which is mainly determined by the trans side. For this process we obtain ß≈1.36. This value can be explained by our derivation of ß=4/3 for constant-bias translocation, where translocated polymer segments form a globule on the trans side. Our results pave the way for understanding and utilizing chaperone-assisted translocation where variations in microscopic details lead to rich variations in the emerging dynamics.

Polymer ejection from strong spherical confinement.

We examine the ejection of an initially strongly confined flexible polymer from a spherical capsid through a nanoscale pore. We use molecular dynamics for unprecedentedly high initial monomer densities. We show that the time for an individual monomer to eject grows exponentially with the number of ejected monomers. By measurements of the force at the pore we show this dependence to be a consequence of the excess free energy of the polymer due to confinement growing exponentially with the number of monomers initially inside the capsid. This growth relates closely to the divergence of mixing energy in the Flory-Huggins theory at large concentration. We show that the pressure inside the capsid driving the ejection dominates the process that is characterized by the ejection time growing linearly with the lengths of different polymers. Waiting time profiles would indicate that the superlinear dependence obtained for polymers amenable to computer simulations results from a finite-size effect due to the final retraction of polymers' tails from capsids.

Criteria for minimal model of driven polymer translocation.

While the characteristics of the driven translocation for asymptotically long polymers are well understood, this is not the case for finite-sized polymers, which are relevant for real-world experiments and simulation studies. Most notably, the behavior of the exponent α, which describes the scaling of the translocation time with polymer length, when the driving force fp in the pore is changed, is under debate. By Langevin dynamics simulations of regular and modified translocation models using the freely jointed-chain polymer model we find that a previously reported incomplete model, where the trans side and fluctuations were excluded, gives rise to characteristics that are in stark contradiction with those of the complete model, for which α increases with fp. Our results suggest that contribution due to fluctuations is important. We construct a minimal model where dynamics is completely excluded to show that close alignment with a full translocation model can be achieved. Our findings set very stringent requirements for a minimal model that is supposed to describe the driven polymer translocation correctly.

Dynamics of polymer ejection from capsid.

Polymer ejection from a capsid through a nanoscale pore is an important biological process with relevance to modern biotechnology. Here, we study generic capsid ejection using Langevin dynamics. We show that even when the ejection takes place within the drift-dominated region there is a very high probability for the ejection process not to be completed. Introducing a small aligning force at the pore entrance enhances ejection dramatically. Such a pore asymmetry is a candidate for a mechanism by which viral ejection is completed. By detailed high-resolution simulations we show that such capsid ejection is an out-of-equilibrium process that shares many common features with the much studied driven polymer translocation through a pore in a wall or a membrane. We find that the ejection times scale with polymer length, τ ∼ N(α). We show that for the pore without the asymmetry the previous predictions corroborated by Monte Carlo simulations do not hold. For the pore with the asymmetry the scaling exponent varies with the initial monomer density (monomers per capsid volume) ρ inside the capsid. For very low densities ρ ≤ 0.002 the polymer is only weakly confined by the capsid, and we measure α = 1.33, which is close to α=1.4 obtained for polymer translocation. At intermediate densities the scaling exponents α = 1.25 and 1.21 for ρ = 0.01 and 0.02, respectively. These scalings are in accord with a crude derivation for the lower limit α = 1.2. For the asymmetrical pore precise scaling breaks down, when the density exceeds the value for complete confinement by the capsid, ρ ⪆ 0.25. The high-resolution data show that the capsid ejection for both pores, analogously to polymer translocation, can be characterized as a multiplicative stochastic process that is dominated by small-scale transitions.

Sedimentation of knotted polymers.

We investigate the sedimentation of knotted polymers by means of stochastic rotation dynamics, a molecular dynamics algorithm that takes hydrodynamics fully into account. We show that the sedimentation coefficient s, related to the terminal velocity of the knotted polymers, increases linearly with the average crossing number n(c) of the corresponding ideal knot. This provides direct computational confirmation of this relation, postulated on the basis of sedimentation experiments by Rybenkov et al. [J. Mol. Biol. 267, 299 (1997)]. Such a relation was previously shown to hold with simulations for knot electrophoresis. We also show that there is an accurate linear dependence of s on the inverse of the radius of gyration R(g)(-1), more specifically with the inverse of the R(g) component that is perpendicular to the direction along which the polymer sediments. When the polymer sediments in a slab, the walls affect the results appreciably. However, R(g)(-1) remains to a good precision linearly dependent on n(c). Therefore, R(g)(-1) is a good measure of a knot's complexity.

Thermal conduction and interface effects in nanoscale Fermi-Pasta-Ulam conductors.

We perform classical nonequilibrium molecular dynamics simulations to calculate heat flow through a microscopic junction connecting two larger reservoirs. In contrast to earlier papers, we also include the reservoirs in the simulated region to study the effect of the bulk-nanostructure interfaces and the bulk conductance. The scalar Fermi-Pasta-Ulam (FPU) model is used to describe the effects of anharmonic interactions in a simple manner. The temperature profile close to the junction in the low-temperature limit is shown to exhibit strong directional features that fade out when temperature increases. Simulating both the FPU chain and the two bulk regions is also shown to eliminate the nonmonotonous temperature variations found for simpler geometries and models. We show that, with sufficiently large reservoirs, the temperature profile in the chain does not depend on the details of thermalization used at the boundaries.

Event distributions of polymer translocation.

We present event distributions for the polymer translocation obtained by extensive Langevin dynamics simulations. Such distributions have not been reported previously and they provide new understanding of the stochastic characteristics of the process. We extract at a high length scale resolution distributions of polymer segments that continuously traverse through a nanoscale pore. The obtained log-normal distributions together with the characteristics of polymer translocation suggest that it is describable as a multiplicative stochastic process. In spite of its clear out-of-equilibrium nature the forced translocation is surprisingly similar to the unforced case. We find forms for the distributions almost unaltered with a common cut-off length. We show that the individual short-segment and short-time movements inside the pore give the scaling relations τ∼Nα and τ∼f-ß for the polymer translocation.

Unforced polymer translocation compared to the forced case.

We present results for unforced polymer translocation from simulations using Langevin dynamics in two dimensions (2D) to four dimensions and stochastic rotation dynamics supporting hydrodynamic modes in three dimensions (3D). We compare our results to forced translocation and a simplified model where the polymer escapes from an infinite pore. The simple model shows that the scaling behavior of unforced translocation is independent of the dimension of the side to which the polymer is translocating. We find that, unlike its forced counterpart, unforced translocation dynamics is insensitive to pore design. Hydrodynamics is seen to markedly speed up the unforced translocation process but not to affect the scaling relations. Average mean-squared displacement shows scaling with average transition time in unforced but not in forced translocation. The waiting-time distribution in unforced translocation follows closely Poissonian distribution. Our measured transfer probabilities align well with those obtained from an equilibrium theory in 3D, but somewhat worse in 2D, where a polymer's relaxation toward equilibrium with respect to its translocation time is slower. Consequently, in stark contrast to forced translocation, unforced translocation is seen to remain close to equilibrium and shows clear universality.

Pore-polymer interaction reveals nonuniversality in forced polymer translocation.

We present a numerical study of forced polymer translocation by using two separate pore models. Both of them have been extensively used in previous forced translocation studies. We show that variations in the pore model affect the forced translocation characteristics significantly in the biologically relevant range of the pore force, i.e., the driving force. Details of the model are shown to change even the obtained scaling relations, which is a strong indication of strongly out-of-equilibrium dynamics in the computational studies which have not yet succeeded in addressing the characteristics of the forced translocation for biopolymers at realistic length scale.

Analysis of DNA elasticity.

With a model that incorporates hydrodynamics directly, we show that flow experiments can be used for detecting some characteristics of the DNA elasticity which manifest themselves clearly at large length scales but cannot be observed by mechanical forcing experiments even at very small length scales. By systematic analysis, the conclusiveness of different experimental methods is evaluated. For the wormlike chain, confirmed as the correct model for DNA, we find an underlying scaling relation between its extension and flow velocity of the form L(p) approximately v(0.155), which emphasizes the significance of hydrodynamics.

Critical evaluation of the computational methods used in the forced polymer translocation.

In forced polymer translocation, the average translocation time tau scales with respect to pore force f and polymer length N as tau approximately f;{-1}N;{beta} . We demonstrate that an artifact in the Metropolis Monte Carlo method resulting in breakage of the force scaling with large f may be responsible for some of the controversies between different computationally obtained results and also between computational and experimental results. Using Langevin dynamics simulations we show that the scaling exponent beta

Dimensional effects in dynamic fragmentation of brittle materials.

It has been shown previously that dynamic fragmentation of brittle D -dimensional objects in a D -dimensional space gives rise to a power-law contribution to the fragment-size distribution with a universal scaling exponent 2-1/D . We demonstrate that in fragmentation of two-dimensional brittle objects in three-dimensional space, an additional fragmentation mechanism appears, which causes scale-invariant secondary breaking of existing fragments. Due to this mechanism, the power law in the fragment-size distribution has now a scaling exponent of approximately 1.17 .

Exponential and power-law mass distributions in brittle fragmentation.

Generic arguments, a minimal numerical model, and fragmentation experiments with gypsum disk are used to investigate the fragment-size distribution that results from dynamic brittle fragmentation. Fragmentation is initiated by random nucleation of cracks due to material inhomogeneities, and its dynamics are pictured as a process of propagating cracks that are unstable against side-branch formation. The initial cracks and side branches both merge mutually to form fragments. The side branches have a finite penetration depth as a result of inherent damping. Generic arguments imply that close to the minimum strain (or impact energy) required for fragmentation, the number of fragments of size s scales as s(-(2D-1)/D) f(1) (- (2/lambda)(D) s)+ f(2) (- s(-1 )(0 ) (lambda+ s(1/D) )(D) ), where D is the Euclidean dimension of the space, lambda is the penetration depth, and f(1) and f(2) can be approximated by exponential functions. Simulation results and experiments can both be described by this theoretical fragment-size distribution. The typical largest fragment size s(0) was found to diverge at the minimum strain required for fragmentation as it is inversely related to the density of initially formed cracks. Our results also indicate that scaling of s(0) close to this divergence depends on, e.g., loading conditions, and thus is not universal. At the same time, the density of fragment surface vanishes as L-1, L being the linear dimension of the brittle solid. The results obtained provide an explanation as to why the fragment-size distributions found in nature can have two components, an exponential as well as a power-law component, with varying relative weights.

Universal dynamic fragmentation in D dimensions.

A generic model is introduced for brittle fragmentation in D dimensions, and this model is shown to lead to a fragment-size distribution with two distinct components. In the small fragment-size limit a scale-invariant size distribution results from a crack branching-merging process. At larger sizes the distribution becomes exponential as a result of a Poisson process, which introduces a large-scale cutoff. Numerical simulations are used to demonstrate the validity of the distribution for D=2. Data from laboratory-scale experiments and large-scale quarry blastings of granitic gneiss confirm its validity for D=3. In the experiments the nonzero grain size of rock causes deviation from the ideal model distribution in the small-size limit. The size of the cutoff seems to diverge at the minimum energy sufficient for fragmentation to occur, but the scaling exponent is not universal.

Comment on "Scaling behavior in explosive fragmentation".

We discuss the data analysis and the conclusions based upon the analysis given in the paper by Diehl et al. [Phys. Rev. E 62, 4742 (2000)].