Source author record

A. Milchev

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

cond-mat.soft cond-mat.stat-mech physics.chem-ph physics.flu-dyn

Catalog footprint

Research overview

18works

4topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2014arXiv

Detachment of semiflexible polymer chains from a substrate - a Molecular Dynamics investigation

Using Molecular Dynamics simulations, we study the force-induced detachment of a coarse-grained model polymer chain from an adhesive substrate. One of the chain ends is thereby pulled at constant speed off the attractive substrate and the resulting saw-tooth profile of the measured mean force $< f >$ vs height $D$ of the end-segment over the plane is analyzed for a broad variety of parameters. It is shown that the observed characteristic oscillations in the $< f >$ - $D$ profile depend on the bending and not on the torsional stiffness of the detached chains. Allowing for the presence of hydrodynamic interactions (HI) in a setup with explicit solvent and DPD-thermostat, rather than the case of Langevin thermostat, one finds that HI have little effect on the $< f >$ - $D$ profile. Also the change of substrate affinity with respect to the solvent from solvophilic to solvophobic is found to play negligible role in the desorption process. In contrast, a changing ratio $ε_s^A / ε_s^B$ of the binding energies of $A$ - and $B$ -segments in the detachment of an $AB$ -copolymer from adhesive surface strongly changes the $< f >$ - $D$ profile whereby the $B$ -spikes vanish when $ε_s^A / ε_s^B < 0.15$ . Eventually, performing an atomistic simulation of a (bio)-polymer {\it polyglycine}, we demonstrate that the simulation results, derived from our coarse-grained model, comply favorably with those from the all-atom simulation.

preprint2014arXiv

Observation of a tricritical wedge filling transition in the 3D Ising model

In this Letter we present evidences of the occurrence of a tricritical filling transition for an Ising model in a linear wedge. We perform Monte Carlo simulations in a double wedge where antisymmetric fields act at the top and bottom wedges, decorated with specific field acting only along the wegde axes. A finite-size scaling analysis of these simulations shows a novel critical phenomenon, which is distinct from the critical filling. We adapt to tricritical filling the phenomenological theory which successfully was applied to the finite-size analysis of the critical filling in this geometry, observing good agreement between the simulations and the theoretical predictions for tricritical filling.

preprint2013arXiv

Controlling the Interactions between Soft Colloids via Surface Adsorption

By employing monomer-resolved computer simulations and analytical considerations based on polymer scaling theory, we analyze the conformations and interactions of multiarm star polymers strongly adsorbed on a smooth, two-dimensional plane. We find a stronger stretching of the arms as well as a stronger repulsive, effective interaction than in the three dimensional case. In particular, the star size scales with the number of arms $f$ as $\sim f^{1/4}$ and the effective interaction as $\sim f^{2}$ , as opposed to $\sim f^{1/5}$ and $\sim f^{3/2}$ , respectively, in three dimensions. Our results demonstrate the dramatic effect that geometric confinement can have on the effective interactions and the subsequent correlations of soft colloids in general, for which the conformation can be altered as a result of geometrical constraints imposed on them.

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.

preprint2013arXiv

Star Polymers Confined in a Nanoslit: A Simulation Test of Scaling and Self-Consistent Field Theories

The free energy cost of confining a star polymer where $f$ flexible polymer chains containing $N$ monomeric units are tethered to a central unit in a slit with two parallel repulsive walls a distance $D$ apart is considered, for good solvent conditions. Also the parallel and perpendicular components of the gyration radius of the star polymer, and the monomer density profile across the slit are obtained. Theoretical descriptions via Flory theory and scaling treatments are outlined, and compared to numerical self-consistent field calculations (applying the Scheutjens-Fleer lattice theory) and to Molecular Dynamics results for a bead-spring model. It is shown that Flory theory and self-consistent field (SCF) theory yield the correct scaling of the parallel linear dimension of the star with $N$ , $f$ and $D$ , but cannot be used for estimating the free energy cost reliably. We demonstrate that the same problem occurs already for the confinement of chains in cylindrical tubes. We also briefly discuss the problem of a free or grafted star polymer interacting with a single wall, and show that the dependence of confining force on the functionality of the star is different for a star confined in a nanoslit and a star interacting with a single wall, which is due to the absence of a symmetry plane in the latter case.

preprint2013arXiv

The escape transition of a compressed star polymer: Self-consistent field predictions tested by simulation

The escape transition of a polymer "mushroom" (a flexible chain grafted to a flat non-adsorbing substrate surface in a good solvent) occurs when the polymer is compressed by a cylindrical piston of radius $R$ , that by far exceeds the chain gyration radius. At this transition, the chain conformation abruptly changes from a two-dimensional self-avoiding walk of blobs (of diameter $H$ , the height of the piston above the substrate) to a "flower conformation", i.e. stretched almost one-dimensional string of blobs (with end-to-end distance $\approx R$ ) and an "escaped" part of the chain, the "crown", outside the piston. The extension of this problem to the case of star polymers with $f$ arms is considered, assuming that the center of the star is grafted to the substrate. The question is considered whether under compression the arms escape all together, or whether there occurs an arm by arm escape under increasing compression. Both self-consistent field calculations and Molecular Dynamics simulations are found to favor the latter scenario.

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

Dynamic Compression of in situ Grown Living Polymer Brush: Simulation and Experiment

A comparative dynamic Monte Carlo simulation study of polydisperse living polymer brushes, created by surface initiated living polymerization, and conventional polymer monodisperse brush, comprising linear polymer chains, grafted to a planar substrate under good solvent conditions, is presented. The living brush is created by end-monomer (de)polymerization reaction after placing an array of initiators on a grafting plane in contact with a solution of initially non-bonded segments (monomers). At equilibrium, the monomer density profile ϕ(z) of the LPB is found to decline as ϕ(z) ~ z^{-α} with the distance from the grafting plane z, while the distribution of chain lengths in the brush scales as c(N) ~ N^{-τ}. The measured values α= 0.64 and τ= 1.70 are very close to those, predicted within the framework of the Diffusion-Limited Aggregation theory, α= 2/3 and τ= 7/4. At varying mean degree of polymerization (from L = 28 to L = 170) and effective grafting density (from σ_g = 0.0625 to σ_g = 1.0), we observe a nearly perfect agreement in the force-distance behavior of the simulated LPB with own experimental data obtained from colloidal probe AFM analysis on PNIPAAm brush and with data obtained by Plunkett et. al., [Langmuir 2006, 22, 4259] from SFA measurements on same polymer.

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'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, "trumpet", "stem-trumpet", 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) "trumpet" (weak adsorption and small pulling force), (ii) "stem-trumpet" (weak adsorption and moderate force), and (iii) "stem" (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

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

Absorption/Expulsion of Oligomers and Linear Macromolecules in a Polymer Brush

The absorption of free linear chains in a polymer brush was studied with respect to chain size $L$ and compatibility $χ$ with the brush by means of Monte Carlo (MC) simulations and Density Functional Theory (DFT) / Self-Consistent Field Theory (SCFT) at both moderate, $σ_g = 0.25$ , and high, $σ_g = 1.00$ , grafting densities using a bead-spring model. Different concentrations of the free chains $0.0625 \le ϕ_o \le 0.375$ are examined. Contrary to the case of $χ= 0$ when all species are almost completely ejected by the polymer brush irrespective of their length $L$ , for $χ< 0$ we find that the degree of absorption (absorbed amount) $Γ(L)$ undergoes a sharp crossover from weak to strong ( $\approx 100%$ ) absorption, discriminating between oligomers, $1\le L\le 8$ , and longer chains. For a moderately dense brush, $σ_g = 0.25$ , the longer species, $L > 8$ , populate predominantly the deep inner part of the brush whereas in a dense brush $σ_g = 1.00$ they penetrate into the "fluffy" tail of the dense brush only. Gyration radius $R_g$ and end-to-end distance $R_e$ of absorbed chains thereby scale with length $L$ as free polymers in the bulk. Using both MC and DFT/SCFT methods for brushes of different chain length $32 \le N \le 256$ , we demonstrate the existence of unique {\em critical} value of compatibility $χ= χ^{c}<0$ . For $χ^{c}(ϕ_o)$ the energy of free chains attains the {\em same} value, irrespective of length $L$ whereas the entropy of free chain displays a pronounced minimum. At $χ^{c}$ all density profiles of absorbing chains with different $L$ intersect at the same distance from the grafting plane. The penetration/expulsion kinetics of free chains into the polymer brush after an instantaneous change in their compatibility $χ$ displays a rather rich behavior. We find three distinct regimes of penetration kinetics of free chains regarding the length $L$ : I ( $1\le L\le 8$ ), II ( $8 \le L \le N$ ), and III ( $L > N$ ), in which the time of absorption $τ$ grows with $L$ at a different rate. During the initial stages of penetration into the brush one observes a power-law increase of $Γ\propto t^α$ with power $α\propto -\ln ϕ_o$ whereby penetration of the free chains into the brush gets {\em slower} as their concentration rises.

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.

preprint2009arXiv

Hydrokinetic simulations of nanoscopic precursor films in rough channels

We report on simulations of capillary filling of high-wetting fluids in nano-channels with and without obstacles. We use atomistic (molecular dynamics) and hydrokinetic (lattice-Boltzmann) approaches which point out clear evidence of the formation of thin precursor films, moving ahead of the main capillary front. The dynamics of the precursor films is found to obey a square-root law as the main capillary front, z^2(t) ~ t, although with a larger prefactor, which we find to take the same value for the different geometries (2D-3D) under inspection. The two methods show a quantitative agreement which indicates that the formation and propagation of thin precursors can be handled at a mesoscopic/hydrokinetic level. This can be considered as a validation of the Lattice-Boltzmann (LB) method and opens the possibility of using hydrokinetic methods to explore space-time scales and complex geometries of direct experimental relevance. Then, LB approach is used to study the fluid behaviour in a nano-channel when the precursor film encounters a square obstacle. A complete parametric analysis is performed which suggests that thin-film precursors may have an important influence on the efficiency of nanochannel-coating strategies.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

A. Milchev

Research overview

Is this your work?

18 publications

Detachment of semiflexible polymer chains from a substrate - a Molecular Dynamics investigation

Observation of a tricritical wedge filling transition in the 3D Ising model

Controlling the Interactions between Soft Colloids via Surface Adsorption

Force spectroscopy of polymer desorption: Theory and Molecular Dynamics simulation

Star Polymers Confined in a Nanoslit: A Simulation Test of Scaling and Self-Consistent Field Theories

The escape transition of a compressed star polymer: Self-consistent field predictions tested by simulation

Driven translocation of a polymer: fluctuations at work

Dynamic Compression of in situ Grown Living Polymer Brush: Simulation and Experiment

Force-induced breakdown of flexible polymerized membrane

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

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

Tension enhancement in branched macromolecules upon adhesion on a solid substrate

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

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

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

Absorption/Expulsion of Oligomers and Linear Macromolecules in a Polymer Brush

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

Hydrokinetic simulations of nanoscopic precursor films in rough channels

See the researcher in context

Building this map preview