<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Respectable Justin/David/Mark</DIV>
<DIV>I used ffgmx force field thinking it as suitable for bonds in the oxygen molecule. And yes it gives diffusion coefficient of SPC water to be 3.5 which matches with that given in the gromacs manual. After equilibrating my system of 1 oxygen and 255 water molecules in NPT ensemble at 1 bar and 298 K i subjected it to the NVT ensemble for about 50 ns. The msd is not a straight line. It increases linearly for first 20 ns and after some ns it decreases linearly giving negative diffusion coefficient. Could you help me out with it? The model contains L-J interaction between oxygen and water molecule. The mdp file used for production run is as follows:<BR><FONT size=2></DIV>
<DIV>;RUN CONTROL parameters.</DIV>
<DIV>integrator = md</DIV>
<DIV>dt = 0.002 ; ps !</DIV>
<DIV>nsteps = 50000000 ; total 1 ns.</DIV>
<DIV>nstcomm = 1</DIV>
<DIV>;LANGEVIN DYNAMICS parameters have not been specified as the integrator is </DIV>
<DIV>;md</DIV>
<DIV>;OUPUT CONTROL parameters.</DIV>
<DIV>nstxout = 250</DIV>
<DIV>nstvout = 1000</DIV>
<DIV>nstfout = 0</DIV>
<DIV>nstlog = 100</DIV>
<DIV>nstenergy = 100</DIV>
<DIV>energygrps = OMOL SOL</DIV>
<DIV>;NEIGHBOUR SEARCHING parameters.</DIV>
<DIV>nstlist = 10</DIV>
<DIV>ns_type = grid</DIV>
<DIV>rlist = 0.9 ;initially 1 I modified to 0.9</DIV>
<DIV>;ELECTROSTATIC and VdW parameters.</DIV>
<DIV>rcoulomb = 0.9</DIV>
<DIV>rvdw = 0.9</DIV>
<DIV>epsilon-r = 1 </DIV>
<DIV>;BERENDSEN TEMPERATURE COUPLING is on in two groups</DIV>
<DIV>Tcoupl = berendsen</DIV>
<DIV>tc-grps = OMOL SOL</DIV>
<DIV>tau_t = 0.001 0.01</DIV>
<DIV>ref_t = 298 298</DIV>
<DIV>;PRESSURE COUPLING is on</DIV>
<DIV>Pcoupl = no;berendsen; berendsen exponential relaxation coupling</DIV>
<DIV>tau_p = 0.1;modified from 0.5 </DIV>
<DIV>compressibility = 4.6e-5</DIV>
<DIV>ref_p = 1.0</DIV>
<DIV>;SIMULATED ANNEALING parameters are not specified.</DIV>
<DIV>;GENERATE VELOCITIES is on at 298 K.</DIV>
<DIV>gen_vel = no; ; generate initially</DIV>
<DIV>gen_temp = 298</DIV>
<DIV>gen_seed = 173529 ;give different values for different trials.</DIV>
<DIV>;BONDS parameters</DIV>
<DIV>constraints = all-bonds</DIV>
<DIV>constraint-algorithm = shake</DIV>
<DIV>unconstrained-start = yes</DIV></FONT>
<DIV>&nbsp;</DIV>
<DIV>and the msd curve is attached in the next mail due to insufficient space</DIV>
<DIV>Your help is necessary for me.</DIV>
<DIV>Sunil Kumar Thapa</DIV>
<DIV>Tribhuvan University, Nepal</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;</DIV></td></tr></table><br>