Theory and Monte Carlo simulations for the stretching of flexible and semiflexible single polymer chains under external fields

Stretching experiments on single molecules of arbitrary length opened the way for studying the statistical mechanics of small systems. In many cases in which the thermodynamic limit is not satisfied, different macroscopic boundary conditions, corresponding to different statistical mechanics ensembles, yield different force-displacement curves. We formulate analytical expressions and develop Monte Carlo simulations to quantitatively evaluate the difference between the Helmholtz and the Gibbs ensembles for a wide range of polymer models of biological relevance. We consider generalizations of the freely jointed chain and of the worm-like chain models with extensible bonds. In all cases we show that the convergence to the thermodynamic limit upon increasing contour length is described by a suitable power law and a specific scaling exponent, characteristic of each model. (C) 2012 American Institute of Physics.