We present a unified scaling description for the dynamics of monomers of a semiflexible chain under good solvent condition in the free draining limit. We consider both the cases where the contour length $L$ is comparable to the persistence length $el
l_p$ and the case $Lgg ell_p$. Our theory captures the early time monomer dynamics of a stiff chain characterized by $t^{3/4}$ dependence for the mean square displacement(MSD) of the monomers, but predicts a first crossover to the Rouse regime of $t^{2 u/{1+2 u}}$ for $tau_1 sim ell_p^3$, and a second crossover to the purely diffusive dynamics for the entire chain at $tau_2 sim L^{5/2}$. We confirm the predictions of this scaling description by studying monomer dynamics of dilute solution of semi-flexible chains under good solvent conditions obtained from our Brownian dynamics (BD) simulation studies for a large choice of chain lengths with number of monomers per chain N = 16 - 2048 and persistence length $ell_p = 1 - 500$ Lennard-Jones (LJ) units. These BD simulation results further confirm the absence of Gaussian regime for a 2d swollen chain from the slope of the plot of $langle R_N^2 rangle/2L ell_p sim L/ell_p$ which around $L/ell_p sim 1$ changes suddenly from $left(L/ell_p right) rightarrow left(L/ell_p right)^{0.5} $, also manifested in the power law decay for the bond autocorrelation function disproving the validity of the WLC in 2d. We further observe that the normalized transverse fluctuations of the semiflexible chains for different stiffness $sqrt{langle l_{bot}^2rangle}/L$ as a function of renormalized contour length $L/ell_p$ collapse on the same master plot and exhibits power law scaling $sqrt{langle l_{bot}^2rangle}/L sim (L/ell_p)^eta $ at extreme limits, where $eta = 0.5$ for extremely stiff chains ($L/ell_p gg 1$), and $eta = -0.25$ for fully flexible chains.
Using a coarse-grained bead-spring model for semi-flexible macromolecules forming a polymer brush, structure and dynamics of the polymers is investigated, varying chain stiffness and grafting density. The anchoring condition for the grafted chains is
chosen such that their first bonds are oriented along the normal to the substrate plane. Compression of such a semi-flexible brush by a planar piston is observed to be a two-stage process: for small compressions the chains contract by buckling deformation whereas for larger compression the chains exhibit a collective (almost uniform) bending deformation. Thus, the stiff polymer brush undergoes a 2-nd order phase transition of collective bond reorientation. The pressure, required to keep the stiff brush at a given degree of compression, is thereby significantly smaller than for an otherwise identical brush made of entirely flexible polymer chains! While both the brush height and the chain linear dimension in the z-direction perpendicular to the substrate increase monotonically with increasing chain stiffness, lateral (xy) chain linear dimensions exhibit a maximum at intermediate chain stiffness. Increasing the grafting density leads to a strong decrease of these lateral dimensions, compatible with an exponential decay. Also the recovery kinetics after removal of the compressing piston is studied, and found to follow a power-law / exponential decay with time. A simple mean-field theoretical consideration, accounting for the buckling/bending behavior of semi-flexible polymer brushes under compression, is suggested.
Polymer translocation through a nano-pore in a thin membrane is studied using a coarse-grained bead-spring model and Langevin dynamics simulation with a particular emphasis to explore out of equilibrium characteristics of the translocating chain. We
analyze the out of equilibrium chain conformations both at the $cis$ and the $trans$ side separately either as a function of the time during the translocation process or as as function of the monomer index $m$ inside the pore. A detailed picture of translocation emerges by monitoring the center of mass of the translocating chain, longitudinal and transverse components of the gyration radii and the end to end vector. We observe that polymer configurations at the $cis$ side are distinctly different from those at the $trans$ side. During the translocation, and immediately afterwards, the chain is clearly out of equilibrium, as different parts of the chain are characterized by a series of effective Flory exponents. We further notice that immediately after the translocation the last set of beads that have just translocated take a relatively compact structure compared to the first set of beads that translocated earlier, and the chain immediately after translocation is described by an effective Flory exponent $0.45 pm 0.01$. The analysis of these results is further strengthened by looking at the conformations of chain segments of equal length as they cross from the $cis$ to the $trans$ side, We discuss implications of these results to the theoretical estimates and numerical simulation studies of the translocation exponent reported by various groups.
We propose a new coarse-grained model for the description of liquid-vapor phase separation of colloid-polymer mixtures. The hard-sphere repulsion between colloids and between colloids and polymers, which is used in the well-known Asakura-Oosawa (AO)
model, is replaced by Weeks-Chandler-Anderson potentials. Similarly, a soft potential of height comparable to thermal energy is used for the polymer-polymer interaction, rather than treating polymers as ideal gas particles. It is shown by grand-canonical Monte Carlo simulations that this model leads to a coexistence curve that almost coincides with that of the AO model and the Ising critical behavior of static quantities is reproduced. Then the main advantage of the model is exploited - its suitability for Molecular Dynamics simulations - to study the dynamics of mean square displacements of the particles, transport coefficients such as the self-diffusion and interdiffusion coefficients, and dynamic structure factors. While the self-diffusion of polymers increases slightly when the critical point is approached, the self-diffusion of colloids decreases and at criticality the colloid self-diffusion coefficient is about a factor of 10 smaller than that of the polymers. Critical slowing down of interdiffusion is observed, which is qualitatively similar to symmetric binary Lennard-Jones mixtures, for which no dynamic asymmetry of self-diffusion coefficients occurs.
A description of phase separation kinetics for solid binary (A,B) mixtures in thin film geometry based on the Kawasaki spin-exchange kinetic Ising model is presented in a discrete lattice molecular field formulation. It is shown that the model descri
bes the interplay of wetting layer formation and lateral phase separation, which leads to a characteristic domain size $ell(t)$ in the directions parallel to the confining walls that grows according to the Lifshitz-Slyozov $t^{1/3}$ law with time $t$ after the quench. Near the critical point of the model, the description is shown to be equivalent to the standard treatments based on Ginzburg-Landau models. Unlike the latter, the present treatment is reliable also at temperatures far below criticality, where the correlation length in the bulk is only of the order of a lattice spacing, and steep concentration variations may occur near the walls, invalidating the gradient square approximation. A further merit is that the relation to the interaction parameters in the bulk and at the walls is always transparent, and the correct free energy at low temperatures is consistent with the time evolution by construction.
A polymer chain containing $N$ monomers confined in a finite cylindrical tube of diameter $D$ grafted at a distance $L$ from the open end of the tube may undergo a rather abrupt transition, where part of the chain escapes from the tube to form a crow
n-like coil outside of the tube. When this problem is studied by Monte Carlo simulation of self-avoiding walks on the simple cubic lattice applying a cylindrical confinement and using the standard pruned-enriched Rosenbluth method (PERM), one obtains spurious results, however: with increasing chain length the transition gets weaker and weaker, due to insufficient sampling of the escaped states, as a detailed analysis shows. In order to solve this problem, a new variant of a biased sequential sampling algorithm with re-sampling is proposed, force-biased PERM: the difficulty of sampling both phases in the region of the first order transition with the correct weights is treated by applying a force at the free end pulling it out of the tube. Different strengths of this force need to be used and reweighting techniques are applied. Using rather long chains (up to N=18000) and wide tubes (up to D=29 lattice spacings), the free energy of the chain, its end-to-end distance, the number of imprisoned monomers can be estimated, as well as the order parameter and its distribution. It is suggested that this new algorithm should be useful for other problems involving state changes of polymers, where the different states belong to rather disjunct valleys in the phase space of the system.
When systems that can undergo phase separation between two coexisting phases in the bulk are confined in thin film geometry between parallel walls, the phase behavior can be profoundly modified. These phenomena shall be described and exemplified by c
omputer simulations of the Asakura-Oosawa model for colloid-polymer mixtures, but applications to other soft matter systems (e.g. confined polymer blends) will also be mentioned. Typically a wall will prefer one of the phases, and hence the composition of the system in the direction perpendicular to the walls will not be homogeneous. If both walls are of the same kind, this effect leads to a distortion of the phase diagram of the system in thin film geometry, in comparison with the bulk, analogous to the phenomenon of capillary condensation of simple fluids in thin capillaries. In the case of competing walls, where both walls prefer different phases of the two phases coexisting in the bulk, a state with an interface parallel to the walls gets stabilized. The transition from the disordered phase to this soft mode phase is rounded by the finite thickness of the film and not a sharp phase transition. However, a sharp transition can occur where this interface gets localized at (one of) the walls. The relation of this interface localization transition to wetting phenomena is discussed. Finally, an outlook to related phenomena is given, such as the effects of confinement in cylindrical pores on the phase behavior, and more complicated ordering phenomena (lamellar mesophases of block copolymers or nematic phases of liquid crystals under confinement).
We revisit the classical problem of a polymer confined in a slit in both of its static and dynamic aspects. We confirm a number of well known scaling predictions and analyse their range of validity by means of comprehensive Molecular Dynamics simulat
ions using a coarse-grained bead-spring model of a flexible polymer chain. The normal and parallel components of the average end-to-end distance, mean radius of gyration and their distributions, the density profile, the force exerted on the slit walls, and the local bond orientation characteristics are obtained in slits of width $D$ = $4 div 10$ (in units of the bead radius) and for chain lengths $N=50 div 300$. We demonstrate that a wide range of static chain properties in normal direction can be described {em quantitatively} by analytic model - independent expressions in perfect agreement with computer experiment. In particular, the observed profile of confinement-induced bond orientation, is shown to closely match theory predictions. The anisotropy of confinement is found to be manifested most dramatically in the dynamic behavior of the polymer chain. We examine the relation between characteristic times for translational diffusion and lateral relaxation. It is demonstrated that the scaling predictions for lateral and normal relaxation times are in good agreement with our observations. A novel feature is the observed coupling of normal and lateral modes with two vastly different relaxation times. We show that the impact of grafting on lateral relaxation is equivalent to doubling the chain length.
