<div>Hi everyone,<br></div><div>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?</div>
<div></div>
<div>Thank you in advance,</div><div>Alexey</div>