Hi everyone, I have a model for liquid benzene. In this model all bonds are treated as harmonic potentials. When I use larger constants for these potentials (to make molecules more planar), the average density of my benzene (from NPT ensemble) becomes lower. Does virial depend on intramolecular interactions? If a want to make benzene molecules completely rigid (through the virtual sites, for example), can I use the same parameters for intermolecular interactions? Thank you in advance, Alexey
_______________________________________________ gmx-users mailing list [email protected] http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [email protected]. Can't post? Read http://www.gromacs.org/mailing_lists/users.php

