Modeling large Gaussian ring polymers

In many types of computer simulation, a ring polymer of length l in a particular solvent is represented as a polygon of N sides of length lp, where N≊l/lp, and lp is the persistence length. There are many established methods to make such polygons, including kink-jump, bead-spring, and dimerization. This paper introduces a very efficient and easy-to-program method, called vector shuffling, which is useful for generating such polygons in the important special case where the quality of the solvent is good, and the polymer is said to be under θ conditions.