Rotational isomeric state chains on a high coordination lattice: Dynamic Monte Carlo algorithm details

A high coordination lattice model for simulating coarse-grained rotational isomeric state (RIS) chains has been under development recently. Initially, the model was developed for chains with symmetric torsional potential ener gy functions, E(phi)=-E(phi). A single-bead move Monte Carlo algorithm was used and found to be effective in simulating polyethylene chains. A modific ation was subsequently developed to allow for the simulation of chains with an asymmetric torsional potential, E(phi)not equal-E(phi). The single-bead move Monte Carlo (MC) algorithm employed previously was found to be ineffe ctive following this modification. Similar kinetic effects have been seen p reviously with single-bead MC moves on the cubic lattice, which lead to the Hilhorst-Deutch modification (two-bead crankshaft MC move) of the Verdier- Stockmayer (single-bead move MC) algorithm. A reptation MC move applied to this model appears problematic. A multiple-bead MC move is developed using the pivot algorithm formalism in order to avoid the lattice model specific kinetic effect seen with only single-bead MC moves. This allows for the eff ective simulation of vinyl polymers with asymmetric torsional potentials su ch as polypropylene and polyvinyl chloride. Polypropylene (PP) and poly(vin yl chloride) (PVC) chains of varying stereochemical structure are simulated . The chains are found to relax with reasonable efficiency. Polypropylene a nd polyvinyl chloride chains are reverse mapped back to the fully atomistic description. The solubility parameters of the reverse-mapped atomistic str uctures are found to agree reasonably well with experimental values. (C) 20 00 American Institute of Physics. [S0021-9606(00)50722-3].