Molecular dynamics simulations of polymers: Methods for optimal Fortran programming

The structure of the general classical trajectory FORTRAN code is discussed for use in molecular dynamics simulations of polymer processes. Substantial reductions in the number of mathematical operations in the code are achieved by using statement functions, modifying the ususal chain rule derivatives, applying appropriate vector identities, and using closure relationships among the required derivatives of the potential.