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.
Users
Please
log in to take part in the discussion (add own reviews or comments).