of freedom. Further reductions in dimensionality are possible by imposing constraints on the relative positions of the atoms in a molecule. Thus we might specify that the bond lengths and bond angles in a molecule are fixed and only the torsion angles are allowed to vary. It is such a scenario that we examine in this paper. We address two aspects of this problem: how to move a molecule subject to constraints in order to allow equilibrium averages to be computed using the canonical-ensemble Monte Carlo method [@metropolis53] and how to evaluate the energy efficiently.

The imposition of constraints in molecular modeling has been extensively studied [@frenkel02 §3.3.2, §15.1]. Let us start by elucidating the difference in the treatment of hard constraints in molecular dynamics and Monte Carlo simulations. We treat hard constraints by taking the limit where the “spring constant” for the hard degrees of freedom is infinite. Molecular dynamics simulations then consider the evolution of the resulting system over a finite time. On the other hand, if we wish to determine the equilibrium properties of a system using the Monte Carlo method, then we need to consider averaging over sufficiently long times to allow equipartition of energy among 