A Molecular Dynamics Technique to Extract Forces in Soft Matter Systems Under Compression With Constant Solvent Chemical Potential

Molecular dynamics simulations of opposing polymer brushes at varying surface separation distances were performed to develop a method for conducting a static compression of soft matter. As all separation distances were represented by independent simulations, the proper solvent density for every level of compression needed to be determined to acquire realistic data. This was accomplished by maintaining a constant solvent chemical potential for each separation distance. In doing so, each independent simulation is equilibrated with all others, reproducing conditions encountered experimentally in force spectroscopy measurements. Chemical potential was determined using the Widom test particle insertion method. Force information was extracted from pressure profiles, such that unphysical forces occurring within the surface layers were not accounted for in the calculation. Each individual simulation was a canonical ensemble molecular dynamics simulation, but taken together they approximate a grand canonical ensemble for the solvent particles by holding their chemical potential constant.