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 s
aw-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 $epsilon_s^A / epsilon_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 $epsilon_s^A / epsilon_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.
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) simu
lations 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 $epsilon$ 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.
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 rad
ius. 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.
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 conditio
ns. 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.
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.
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 $sigma_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 $epsilon_s$? We examine the monomer density distribution profiles across the BB-backbone at different $epsilon_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 $epsilon_s$, too long side chains mutually block each other and freeze effectively the bottle-brush molecule.
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 invest
igate 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 $<tau>$ of bonds decays with membrane size as $<tau> propto N^{-beta}$ where $beta approx 0.50pm 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 / <tau>)$. The mean failure time $<tau_r>$ that takes to rip-off the sheet declines with growing size $N$ as a power law $<tau_r> propto N^{-phi(f)}$. We also find $<tau_r> propto exp(Delta U_0/k_BT)$ where the nucleation barrier for crack formation $Delta U_0 propto f^{-2}$, in agreement with Griffiths theory. $<tau_r>$ displays an Arrhenian dependence of $<tau_r>$ on temperature $T$. Our results indicate a rapid increase in crack spreading velocity with growing external tension $f$.
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 dist
ribution 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 $tau$ is found to decrease with the total number of network nodes $N$ by a power law $tau 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 $delta$-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.
