3
$\begingroup$

To extract the force as in the title I used the gmx trjconv command like below.

gmx trjconv -f alad.trr -s alad.tpr -n alad.ndx -pbc mol -novel -force -o alad.gro

However, the force was not extracted in alad.gro file. Of course, I got the coordinates and the force at equal intervals when getting the alad.trr file.

Another way, I used gmx traj command to get force.

gmx traj -f alad.trr -s alad.tpr -n alad.ndx -pbc -mol -ox coordi.xvg -of force.xvg

However, the force obtained from force.xvg didn't look like the value I wanted.

So, how do I get the force of molecules put in a single periodic box?

$\endgroup$
  • $\begingroup$ Could you include a figure showing force.xvg, and a better description of what it is exactly that you would like to measure, and why you think the output is not correct. $\endgroup$ – Bdrs Mar 2 at 15:45
  • $\begingroup$ However, I'm still wondering why coordinates of alad gro and coori.xvg are not same. $\endgroup$ – Yeonju Mar 5 at 11:49
  • $\begingroup$ I solved the problem. As gmx trjconv -pbc mol option only make molecules translate, the cooridinates obtained from alad.gro can correspond to forces obtained from gmx traj -f alad.trr -s alad.tpr -n alad.ndx -of force-dump.xvg(same results using gmx dump command). However, coordinates of coordi.xvg doesn't match with alad.gro . Also I couldn't trust force.xvg. whether the data is correct is tested by my own work ouput. So it's hard to say simple.. $\endgroup$ – Yeonju Mar 5 at 11:52
  • 2
    $\begingroup$ There is a proposal for a Stack Exchange site for Materials Modeling which is currently in the commitment phase. There are a lot of experts currently there and questions like these would get more activity there. It would be great if you could help building the site by committing to the proposal. $\endgroup$ – rashid Mar 5 at 11:59
0
$\begingroup$

I solved the problem. As gmx trjconv -pbc mol option only make molecules translate, the cooridinates obtained from alad.gro can correspond to forces obtained from gmx traj -f alad.trr -s alad.tpr -n alad.ndx -of force-dump.xvg.

|improve this answer|||||
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.