Molecular dynamics simulations are useful tools to unveil molecular mechanisms of polymer phase separation,self-assembly,adsorption,and so on.Due to large molecular size and slow relaxation of the polymer chains,a great amount of issues related to large-distance chain displacement cannot be tackled easily with conventional molecular dynamic simulations.Systematic coarse-graining and enhanced sampling methods are two types of improvements that can boost spatiotemporal scales in polymer simulations.We present two typical ways to obtain the coarse-graining potential either by fitting to correct liquid structures or by fitting to available thermodynamic properties of polymer systems.The newly proposed anisotropic coarse-grained particle model can be used to describe aggregation and assembly of polymeric building blocks from disk-like micelles to Janus particles.We also present a stochastic polymerization model combined with coarse-grained simulations to investigate the problems strongly influenced by the coupling of polymerization and excluded volume effects.Finally,a facile implementation of integrated tempering sampling method is illustrated to be very efficient on bypassing local energy minima and having access to true equilibrium polymer structures.
LI YanChunWANG YongLeiLI ZhanWeiLIU HongLü ZhongYuan
In this paper, we present the coil-to-globule(CG) transitions of homopolymers and multiblock copolymers with different topology and stiffness by using molecular dynamics with integrated tempering sampling method. The sampling method was a novel enhanced method that efficiently sampled the energy space with low computational costs. The method proved to be efficient and precise to study the structural transitions of polymer chains with complex topological constraint, which may not be easily done by using conventional Monte Carlo method. The topological constraint affects the globule shape of the polymer chain, thus further influencing the CG transition. We found that increasing the topological constraint generally decreased CG transition temperature for homopolymers. For semiflexible chains, an additional first-order like symmetry-broken transition emerged. For block copolymers, the topological constraint did not obviously change the transition temperature, but greatly reduced the energy signal of the CG transition.