A new simulation method for the determination of forces in polymer/colloid systems

We introduce a new simulation method, which we call the contact-distribution method, for the determination of the Helmholtz potential for polymer/colloid systems from lattice Monte-Carlo simulations. This method allows one to obtain forces between finite or semi-infinite objects of any arbitrary shape and dimensions in the presence of polymer chains in solution or physisorbed or chemisorbed at interfaces. We illustrate the application of the method using two examples: (i) the interaction between the tip of an atomic force microscope (AFM) and a single, end-grafted polymer chain and (ii) the interaction between an AFM tip and a polymer brush. Numerical results for the first two cases illustrate how the method can be used to confirm and extend scaling laws for forces and Helmholtz potentials, to examine the effects of the shapes and sizes of the objects and to examine conformational transitions in the polymer chains.