Phase transitions of single semistiff polymer chains

被引:98
作者
Bastolla, U [1 ]
Grassberger, P [1 ]
机构
[1] Forschungszentrum Julich, HLRZ, D-52425 Julich, Germany
关键词
polymers; protein folding; phase transitions;
D O I
10.1007/BF02764222
中图分类号
O4 [物理学];
学科分类号
0702 ;
摘要
We study numerically a lattice model of semiflexible homopolymers with nearest neighbor (nn) attraction and energetic preference for straight joints between bonded monomers. For this we use a new Monte Carlo algorithm, the "pruned-enriched Rosenbluth Method" (PERM). It is very efficient both for relatively open configurations at high temperatures and for compact and frozen-in low-T states. This allows us to study in detail the phase diagram as a function of nn attraction epsilon and stiffness x. It shows a theta-collapse line with a transition from open coils (small epsilon) to molten compact globules (large epsilon) and a freezing transition toward a state with orientational global order (large stiffness x). Qualitatively this is similar to a recently studied mean-field theory [S. Doniach, T. Garel, and H. Orland (1996), J. Chem. Phys. 105(4), 1601], but there are important differences in details. In contrast to the mean-field theory and to naive expectations, the theta-temperature ino eases with stiffness x. The freezing temperature increases even faster, and reaches the theta-line at a finite value of x. For even stiffer chains, the freezing transition takes place directly, without the formation Of an intermediate globular state. Although being in conflict with mean-field theory, the latter had been conjectured already by Doniach ct ai. on the basis of heuristic arguments and of low-statistics Monte Carlo simulations. Finally, we discuss the relevance of the present model as a very crude model for protein folding.
引用
收藏
页码:1061 / 1078
页数:18
相关论文
共 22 条