Molecular dynamics simulation of a generic polymer model is applied to study melts of polymers with different types of intrinsic stiffness. Important static observables of the single chain such as gyration radius or persistence length are determined. Additionally we investigate the overall static melt structure including pair correlation function, structure function and orientational correlation function.

A recently proposed method for computer simulations in the isothermalisobaric (NPT) ensemble, based on Langevin-type equations of motion for the particle coordinates and the "piston" degree of freedom, is re-derived by straightforward application of the standard Kramers-Moyal formalism. An integration scheme is developed which reduces to a time-reversible symplectic integrator in the limit of vanishing friction. This algorithm is hence expected to be quite stable for small friction, allowing for a large time step. We discuss the optimal choice of parameters, and present some numerical test results.

No abstract

No abstract