Dear all,<br>&nbsp; Recently I made some tests, and I want to get same potential with the following method. (gmx 4.0.5.)<br>&nbsp; Firstly, I run a complete simulation as the reference. Then I got a log file. <br>&nbsp; Secondly, I use trjconv(-sep) to get each frame to a separate *.gro file.<br>&nbsp; Thirdly, I use *.gro and run md for just one step (delta t=2fs), and another input files are same with the reference one. Then I just want to read potential in step 0 in the log file.&nbsp;&nbsp; <br>&nbsp; I compare the potential which I got from different log files, and I found that there are big difference between them. But in my opinion, as they have the same structure and topol.top they should have the same (or little difference) Energies like "G96Angle, Proper Dih., Improper Dih., LJ-14, Coulomb-14, LJ (SR), Disper. corr., Coulomb (SR), Coul. recip., and Potential".<br>&nbsp; <br>&nbsp; The following is the job script lines :<br>&nbsp; ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>&nbsp; # Reference<br>&nbsp; grompp -f md.mdp(nsteps=50000) -c npt.gro -t npt.cpt -p topol.top -o md.tpr<br>&nbsp; mpirun -np 8 mdrun_mpi -deffnm md<br><br>&nbsp; # Comparison<br>&nbsp; echo 0|trjconv -f md.xtc -s md.tpr -o ./md_gro/md.gro -sep -pbc mol -ur compact<br>&nbsp; grompp -f md1.mdp(nsteps=1) -c ./md_gro/mdi.gro(i=0,1,...) -p topol.top -o ./tpr/md_1_i.tpr<br>&nbsp; mpirun -np 8 mdrun_mpi -s ./tpr/md_1_i.tpr -g ./log/md_1_i.log<br>&nbsp; ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>&nbsp; Do you have some further ideas?<br>&nbsp; Thanks in advance!<br><br>&nbsp; best wishes,<br>&nbsp; lin.<br><br><br><br><br><a href="http://t.sina.com.cn/1324608620?s=RGFuZGVsaW9uMDIxMg==" target="_blank"><br></a>