[gmx-users] box of methanol

Erik Lindahl lindahl at stanford.edu
Fri Jul 12 23:04:12 CEST 2002

nicolis at guideo.fr wrote:

>>You don't need to know the exact compressibility, just the order of 
>>magnitude. The only place where it enters is the pressure scaling, but 
>>it only affects the speed of scaling, not the equilibrium pressure.
>OK Erik, thanks for the answer, but in that case, there must be something I 
>am doing wrong because I don't get the right density. After 1000ps with the 
>methanol.itp of the $GMXDATA/top and the ffG43a1 forcefield, I end with a 
>density of 850 instead of 791.
>Maybe the 3x3x3 box (550 methanols) is too small?
Not sure; you should probably check the forcefield article and check 
their results - I
guess Wilfred and the gang included methanol in the tests.

Also make sure you use the right temperature and cut-off parameters. 
Ewald summation
will often give a too high density; Gromos-96 was parameterized with 
0.8/1.4 twin-range

You might also want to try OPLS, but there's no methanol topology there 
yet, so it has
to be handcoded (with all dihedrals).



More information about the gromacs.org_gmx-users mailing list