Incollection,

Osmotic equation of state for solutions of semi-flexible macromolecules: Monte Carlo computer simulation

, , , , , , and .
Abstract Book of the XXIII IUPAP International Conference on Statistical Physics, Genova, Italy, (9-13 July 2007)

Abstract

We present results of lattice Monte Carlo computer simulations of solutions of semi-flexible macromolecules using the bond fluctuation model. We have studied athermal solutions of short macromolecules (chain length N=20) with various choices for the chain stiffness in the concentration range from very dilute to concentrated solutions where the appearance of nematic phase becomes possible. Simulations were performed both in canonical and grand canonical ensemble, and we used the configurational bias Monte Carlo algorithm in the latter case. Three techniques to obtain the osmotic equation of state in Monte Carlo simulations are applied and compared in order to critically assess their efficiency and accuracy: the ``repulsive wall'' method, the thermodynamic integration method (which rests on the feasibility of simulations in the grand canonical ensemble), and the recently advocated sedimentation equilibrium method, which records the density profile in an external (e.g. gravitation-like) field and infers, via a local density approximation, the equation of state from the hydrostatic equilibrium condition. We confirm the conclusion that the latter technique is far more efficient than the repulsive wall method, but we find that the thermodynamic integration method is similarly efficient as the sedimentation equilibrium method, and seems to work better for very stiff chains (where the onset of nematic order and the formation of isotropic--nematic interfaces renders the use of the sedimentation equilibrium method problematic at large densities). The effect of a nematic wetting layer at a surface is also discussed.

Tags

Users

  • @statphys23

Comments and Reviews