6533b7ddfe1ef96bd127478f
RESEARCH PRODUCT
Monte Carlo simulation of crystalline polyethylene
Roman MartoňákWolfgang PaulKurt Bindersubject
PhysicsQuantum Monte CarloMonte Carlo methodDegrees of freedom (physics and chemistry)General Physics and AstronomyHybrid Monte Carlosymbols.namesakeMetropolis–Hastings algorithmHardware and ArchitectureDynamic Monte Carlo methodsymbolsStatistical physicsvan der Waals forceMonte Carlo molecular modelingdescription
Abstract We consider here the problem of constructing an efficient algorithm for a classical Monte Carlo simulation of crystalline polyethylene with unconstrained bond lengths and angles. This macromolecular crystal presents a particular example of a system with many different energy scales, ranging from soft ones represented by nonbonded van der Waals interactions, to stiff ones, represented in particular by bond stretching. A proper sampling of all the energy scales poses a problem and it is shown that a standard Metropolis algorithm employing just local moves is not very efficient at low temperatures. As a solution it is proposed to employ also global moves consisting of displacements of the center of mass of the whole chains in all three spatial directions as well as rotations of the chains around an axis parallel to the z -axis, which act on the degrees of freedom associated with the chain packing. It is shown that by properly alternating such global moves with standard local moves, the statistical inefficiency of the algorithm is considerably reduced.
| year | journal | country | edition | language |
|---|---|---|---|---|
| 1996-12-01 | Computer Physics Communications |