Free Energy Calculations by Expanded Ensemble Method for Lattice and Continuous Polymers

Abstract

Expanded ensemble Monte Carlo method suggested earlier in J. Chem. Phys. (1992, 96, 1776) was applied to calculate free energy (entropy) of a lattice polymer (self-avoiding random walk on a simple cubic lattice) with respect to the standard reference systema pure random walk with no restrictions (in the latter case the number of conformations is known exactly). An effective scheme for the initial choice of balancing factors is proposed with the following iteration procedure for their optimization. Calculations were performed for chain lengths N = 5, 8, 25, 50, and 100. Comparison of simulation results with the exact data available for short chains (N = 5, 8) indicate a high accuracy of the simulation method. For the continuous polymer model (several bead and spring chains in the periodic box of a fixed volume) the expanded ensemble Monte Carlo method was used to calculate free energy with respect to the ideal gas as a reference system. Configurations were changed with the aid of the constant temperature m...