Computation of pressure components due to Class II force fields.
It has been shown that purely angle dependent terms in a force field do not contribute to the total pressure in a molecular simulation. However, this is not the case for the individual components of the pressure tensor, and is also untrue for crossterms in the force field including bond stretch. In this article, we show that virial contributions to the pressure tensor are easily computed in terms of bond distance vectors and atom forces for the bond topologies present in a Class II force field. Results from a simulation of a phospholipid biomembrane using the cff97 force field show that angle and torsional crossterms make a significant contribution to the pressure tensor.
Carpenter
John E
JE
SGI, Eagan, Minnesota 55121, USA. jec@mayo.edu