Researcher profile

T. A. Vilgis

T. A. Vilgis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2014arXiv

Driven translocation of a polymer: role of pore friction and crowding

Force-driven translocation of a macromolecule through a nanopore is investigated by taking into account the monomer-pore friction as well as the "crowding" of monomers on the {\it trans} - side of the membrane which counterbalance the driving force acting in the pore. The set of governing differential-algebraic equations for the translocation dynamics is derived and solved numerically. The analysis of this solution shows that the crowding of monomers on the trans side hardly affects the dynamics, but the monomer-pore friction can substantially slow down the translocation process. Moreover, the translocation exponent $α$ in the translocation time - vs. - chain length scaling law, $τ\propto N^α$, becomes smaller when monomer-pore friction coefficient increases. This is most noticeable for relatively strong forces. Our findings may explain the variety of $α$ values which were found in experiments and computer simulations.

preprint2013arXiv

Force spectroscopy of polymer desorption: Theory and Molecular Dynamics simulation

Forced detachment of a single polymer chain, strongly-adsorbed on a solid substrate, is investigated by two complementary methods: a coarse-grained analytical dynamical model, based on the Onsager stochastic equation, and Molecular Dynamics (MD) simulations with Langevin thermostat. The suggested approach makes it possible to go beyond the limitations of the conventional Bell-Evans model. We observe a series of characteristic force spikes when the pulling force is measured against the cantilever displacement during detachment at constant velocity $v_c$ (displacement control mode) and find that the average magnitude of this force increases as $v_c$ grows. The probability distributions of the pulling force and the end-monomer distance from the surface at the moment of final detachment are investigated for different adsorption energy $ε$ and pulling velocity $v_c$. Our extensive MD-simulations validate and support the main theoretical findings. Moreover, the simulation reveals a novel behavior: for a strong-friction and massive cantilever the force spikes pattern is smeared out at large $v_c$. As a challenging task for experimental bio-polymers sequencing in future we suggest the fabrication of stiff, super-light, nanometer-sized AFM probe.

preprint2012arXiv

Driven translocation of a polymer: fluctuations at work

The impact of thermal fluctuations on the translocation dynamics of a polymer chain driven through a narrow pore has been investigated theoretically and by means of extensive Molecular-Dynamics (MD) simulation. The theoretical consideration is based on the so-called velocity Langevin (V-Langevin) equation which determines the progress of the translocation in terms of the number of polymer segments, $s(t)$, that have passed through the pore at time $t$ due to a driving force $f$. The formalism is based only on the assumption that, due to thermal fluctuations, the translocation velocity $v=\dot{s}(t)$ is a Gaussian random process as suggested by our MD data. With this in mind we have derived the corresponding Fokker-Planck equation (FPE) which has a nonlinear drift term and diffusion term with a {\em time-dependent} diffusion coefficient $D(t)$. Our MD simulation reveals that the driven translocation process follows a {\em super}diffusive law with a running diffusion coefficient $D(t) \propto t^γ$ where $γ< 1$. This finding is then used in the numerical solution of the FPE which yields an important result: for comparatively small driving forces fluctuations facilitate the translocation dynamics. As a consequence, the exponent $α$ which describes the scaling of the mean translocation time $<τ>$ with the length $N$ of the polymer, $<τ> \propto N^α$ is found to diminish. Thus, taking thermal fluctuations into account, one can explain the systematic discrepancy between theoretically predicted duration of a driven translocation process, considered usually as a deterministic event, and measurements in computer simulations. In the non-driven case, $f=0$, the translocation is slightly subdiffusive and can be treated within the framework of fractional Brownian motion (fBm).

preprint2012arXiv

Force-induced breakdown of flexible polymerized membrane

We consider the fracture of a free-standing two-dimensional (2D) elastic-brittle network to be used as protective coating subject to constant tensile stress applied on its rim. Using a Molecular Dynamics simulation with Langevin thermostat, we investigate the scission and recombination of bonds, and the formation of cracks in the 2D graphene-like hexagonal sheet for different pulling force $f$ and temperature $T$. We find that bond rupture occurs almost always at the sheet periphery and the First Mean Breakage Time $<τ>$ of bonds decays with membrane size as $<τ> \propto N^{-β}$ where $β\approx 0.50\pm 0.03$ and $N$ denotes the number of atoms in the membrane. The probability distribution of bond scission times $t$ is given by a Poisson function $W(t) \propto t^{1/3} \exp (-t / <τ>)$. The mean failure time $<τ_r>$ that takes to rip-off the sheet declines with growing size $N$ as a power law $<τ_r> \propto N^{-ϕ(f)}$. We also find $<τ_r> \propto \exp(ΔU_0/k_BT)$ where the nucleation barrier for crack formation $ΔU_0 \propto f^{-2}$, in agreement with Griffith&#39;s theory. $<τ_r>$ displays an Arrhenian dependence of $<τ_r>$ on temperature $T$. Our results indicate a rapid increase in crack spreading velocity with growing external tension $f$.

preprint2012arXiv

Forced translocation of a polymer: dynamical scaling vs. MD-simulation

We suggest a theoretical description of the force-induced translocation dynamics of a polymer chain through a nanopore. Our consideration is based on the tensile (Pincus) blob picture of a pulled chain and the notion of propagating front of tensile force along the chain backbone, suggested recently by T. Sakaue. The driving force is associated with a chemical potential gradient that acts on each chain segment inside the pore. Depending on its strength, different regimes of polymer motion (named after the typical chain conformation, &#34;trumpet&#34;, &#34;stem-trumpet&#34;, etc.) occur. Assuming that the local driving and drag forces are equal (i.e., in a quasi-static approximation), we derive an equation of motion for the tensile front position $X(t)$. We show that the scaling law for the average translocation time $<τ>$ changes from $<τ> \sim N^{2ν}/f^{1/ν}$ to $<τ> \sim N^{1+ν}/f$ (for the free-draining case) as the dimensionless force ${\widetilde f}_{R} = a N^νf /T$ (where $a$, $N$, $ν$, $f$, $T$ are the Kuhn segment length, the chain length, the Flory exponent, the driving force, and the temperature, respectively) increases. These and other predictions are tested by Molecular Dynamics (MD) simulation. Data from our computer experiment indicates indeed that the translocation scaling exponent $α$ grows with the pulling force ${\widetilde f}_{R}$) albeit the observed exponent $α$ stays systematically smaller than the theoretically predicted value. This might be associated with fluctuations which are neglected in the quasi-static approximation.

preprint2012arXiv

Polymer Detachment Kinetics from Adsorbing Surface: Theory, Simulation and Similarity to Infiltration into Porous Medium

The force-assisted desorption kinetics of a macromolecule from adhesive surface is studied theoretically, using the notion of tensile (Pincus) blobs, as well as by means of Monte-Carlo (MC) and Molecular Dynamics (MD) simulations. We show that the change of detached monomers with time is governed by a differential equation which is equivalent to the nonlinear porous medium equation (PME), employed widely in transport modeling of hydrogeological systems. Depending on the pulling force and the strength of adsorption, three kinetic regimes can be distinguished: (i) &#34;trumpet&#34; (weak adsorption and small pulling force), (ii) &#34;stem-trumpet&#34; (weak adsorption and moderate force), and (iii) &#34;stem&#34; (strong adsorption and large force). Interestingly, in all regimes the number of desorbed beads $M(t)$, and the height of the first monomer (which experiences a pulling force) $R(t)$ above the surface follow an universal square-root-of-time law. Consequently, the total time of detachment $<τ_d>$, scales with polymer length $N$ as $<τ_d> \propto N^2$. Our main theoretical conclusions are tested and found in agreement with data from extensive MC- and MD-simulations.

preprint2012arXiv

Tension enhancement in branched macromolecules upon adhesion on a solid substrate

The effect of self-generated tension in the backbone of a bottle-brush (BB) macromolecule, adsorbed on an attractive surface, is studied by means of Molecular Dynamics simulations of a coarse-grained bead-spring model in the good solvent regime. The BB-molecule is modeled as a backbone chain of $L$ beads, connected by breakable bonds and with side chains, tethered pairwise to each monomer of the backbone. Our investigation is focused on several key questions that determine the bond scission mechanism and the ensuing degradation kinetics: how are frequency of bond scission and self-induced tension distributed along the BB-backbone at different grafting density $σ_g$ of the side chains? How does tension $f$ depend on the length of the side chains $N$, and on the strength of surface adhesion $ε_s$? We examine the monomer density distribution profiles across the BB-backbone at different $ε_s$ and relate it to adsorption-induced morphological changes of the macromolecule whereby side chains partially desorb while the remaining chains spread better on the surface. Our simulation data are found to be in qualitative agreement with experimental results and recent theoretical predictions. Yet we demonstrate that the interval of parameter values where these predictions hold is limited in $N$. Thus, at high values of $ε_s$, too long side chains mutually block each other and freeze effectively the bottle-brush molecule.

preprint2012arXiv

Thermal decomposition of a honeycomb-network sheet - A Molecular Dynamics simulation study

The thermal degradation of a graphene-like two-dimensional triangular membrane with bonds undergoing temperature-induced scission is studied by means of Molecular Dynamics simulation using Langevin thermostat. We demonstrate that the probability distribution of breaking bonds is highly peaked at the rim of the membrane sheet at lower temperature whereas at higher temperature bonds break at random anywhere in the hexagonal flake. The mean breakage time $τ$ is found to decrease with the total number of network nodes $N$ by a power law $τ\propto N^{-0.5}$ and reveals an Arrhenian dependence on temperature $T$. Scission times are themselves exponentially distributed. The fragmentation kinetics of the average number of clusters can be described by first-order chemical reactions between network nodes $n_i$ of different coordination. The distribution of fragments sizes evolves with time elapsed from a $δ$-function through a bimodal one into a single-peaked again at late times. Our simulation results are complemented by a set of $1^{st}$-order kinetic differential equations for $n_i$ which can be solved exactly and compared to data derived from the computer experiment, providing deeper insight into the thermolysis mechanism.

preprint2011arXiv

Dynamics of pulled desorption with effects of excluded volume interaction: The p-Laplacian diffusion equation and its exact solution

We analyze the dynamics of desorption of a polymer molecule which is pulled at one of its ends with force $f$, trying to desorb it. We assume a monomer to desorb when the pulling force on it exceeds a critical value $f_{c}$. We formulate an equation for the average position of the $n^{th}$ monomer, which takes into account excluded volume interaction through the blob-picture of a polymer under external constraints. The approach leads to a diffusion equation with a $p$-Laplacian for the propagation of the stretching along the chain. This has to be solved subject to a moving boundary condition. Interestingly, within this approach, the problem can be solved exactly in the trumpet, stem-flower and stem regimes. In the trumpet regime, we get $τ=τ_{0}n_d^{2}$ where $n_d$ is the number of monomers that have desorbed at the time $τ$. $τ_{0}$ is known only numerically, but for $f$ close to $f_{c}$, it is found to be $τ_{0}\sim f_c/(f^{2/3}-f_{c}^{2/3})$. If one used simple Rouse dynamics, this result changes to {\normalsize $τ\sim f_c n_d^2/(f-f_{c})$.} In the other regimes too, one can find exact solution, and interestingly, in all regimes $τ\sim n_d^2$.

preprint2011arXiv

Dynamics of two topologically entangled chains

Starting from a given topological invariant, we argue that it is possible to construct a topological field theory with a finite number of Feynman diagrams and an amplitude of gauge invariant objects that is a function of that invariant. This is for example the case of the Gauss linking number and of the abelian BF models which has been already successfully applied in the statistical mechanics of polymers. In this work it is shown that a suitable generalization of the BF model can be applied also to polymer dynamics, where the polymer trajectories are not static, but change their shape during time.

preprint2011arXiv

Polymer chain scission at constant tension - an example of force-induced collective behaviour

The breakage of a polymer chain of segments, coupled by anharmonic bonds with applied constant external tensile force is studied by means of Molecular Dynamics simulation. We show that the mean life time of the chain becomes progressively independent of the number of bonds as the pulling force grows. The latter affects also the rupture rates of individual bonds along the polymer backbone manifesting the essential role of inertial effects in the fragmentation process. The role of local defects, temperature and friction in the scission kinetics is also examined.

preprint2011arXiv

Thermal Degradation of Unstrained Single Polymer Chain: Non-linear Effects at Work

We examine the thermally-induced fracture of an unstrained polymer chain of discrete segments coupled by an anharmonic potential by means of Molecular Dynamics simulation with a Langevin thermostat. Cases of both under- and over-damped dynamics are investigated, and a comparison with recent studies of bond scission in model polymers with harmonic interactions is performed. We find that the polymer degradation changes qualitatively between the inertial regime and that of heavily damped dynamics. The role of bond healing (recombination) is also studied and probability distributions for the recombination times and overstretched bond lengths are obtained. Our extensive simulations reveal many properties of the scission dynamics in agreement with the notion of random breakdown of independent bonds, e.g., the mean time of chain rupture, $<τ>$ follows an Arrhenian behavior with temperature $T$, and depends on the number of bonds $N$ in the polymer as $<τ> \propto N^{-1}$. In contrast, the rupture rates of the individual bonds along the polymer backbone indicate clearly the presence of self-induced inhomogeneity resulting from the interplay of thermal noise and nonlinearity. Eventually we examine the fragmentation kinetics during thermolysis. We demonstrate that both the probability distribution function of fragment sizes as well as the mean length of fragments at subsequent times $t$ characterize degradation as predominantly a first order reaction.

preprint2010arXiv

Thermal Breakage and Self-Healing of a Polymer Chain under Tensile Stress

We consider the thermal breakage of a tethered polymer chain of discrete segments coupled by Morse potentials under constant tensile stress. The chain dynamics at the onset of fracture is studied analytically by Kramers-Langer multidimensional theory and by extensive Molecular Dynamics simulations in 1D- and 3D-space. Comparison with simulation data in one- and three dimensions demonstrates that the Kramers-Langer theory provides good qualitative description of the process of bond-scission as caused by a {\em collective} unstable mode. We derive distributions of the probability for scission over the successive bonds along the chain which reveal the influence of chain ends on rupture in good agreement with theory. The breakage time distribution of an individual bond is found to follow an exponential law as predicted by theory. Special attention is focused on the recombination (self-healing) of broken bonds. Theoretically derived expressions for the recombination time and distance distributions comply with MD observations and indicate that the energy barrier position crossing is not a good criterion for true rupture. It is shown that the fraction of self-healing bonds increases with rising temperature and friction.

preprint1996arXiv

Size and scaling in ideal polymer networks

The scattering function and radius of gyration of an ideal polymer network are calculated depending on the strength of the bonds that form the crosslinks. Our calculations are based on an {\it exact} theorem for the characteristic function of a polydisperse phantom network that allows for treating the crosslinks between pairs of randomly selected monomers as quenched variables without resorting to replica methods. From this new approach it is found that the scattering function of an ideal network obeys a master curve which depends on one single parameter $x= (ak)^2 N/M$, where $ak$ is the product of the persistence length times the scattering wavevector, $N$ the total number of monomers and $M$ the crosslinks in the system. By varying the crosslinking potential from infinity (hard $δ$-constraints) to zero (free chain), we have also studied the crossover of the radius of gyration from the collapsed regime where $R_{\mbox{\tiny g}}\simeq {\cal O}(1)$ to the extended regime $R_{\mbox{\tiny g}}\simeq {\cal O}(\sqrt{N})$. In the crossover regime the network size $R_{\mbox{\tiny g}}$ is found to be proportional to $(N/M)^{1/4}$.