I am running a variety of NPT simulations with polar, non-polar, and ionic compounds.  Although my results for density agree well with experimental values, the pressures I get from g_energy are off by 1 to 3 orders of magnitude.  In the log file, the pressure fluctuates around a lot from -400 to 400 bar, which seems to be normal according to other posts on this list, but the average (which is what g_energy gives me) is not 1.0 bar, as I specified.  Does anyone know how to correct this problem?<div>
<br></div><div>Pressure coupling parameters:</div><div><div>Pcoupl              =  berendsen</div><div>pcoupltype          =  isotropic</div><div>tau_p                   =  1.0</div><div>ref_p                   =  1.0</div>
<div>compressibility     =  4.5e-5</div></div>