Force from gmx traj

GROMACS version:2020.4-MODIFIED
GROMACS modification: Yes
Here post your question

  1. I used gmx traj -fo to evaluate force, but the output has a unit of “Force (kJ mol\S-1\N nm\S-1\N) .”
    In the command-line description, it wasn’t clearly explained how force is calculated and why it has a strange unit.

  2. What is the difference between -cf and -of? I read that :
    " Options -cv and -cf write the average velocities and average forces as temperature factors to a .pdb file with the average coordinates or the coordinates at -ctime. The temperature factors are scaled such that the maximum is 10."
    I don’t understand the purpose of a pdf file and what “temperature factor scaling” means.

Thank you,
Di

Hi,
with options -of the tools extracted the forces storaged in traj.trr and write it in xvg format. During molecular dynamics simulation, forces on each atom are computed as -dV/dr (see Molecular Dynamics — GROMACS 2021.3 documentation) and storaged in traj.trr (based on the setting in the mdp file)

The unit may look strange (in xmgrace format), but the unit is [kJ]/[mol nm] in line with unit used in GROMACS (Definitions and Units — GROMACS 2021.3 documentation )

the option -cf write the average forces as temperature factors to a .pdb file.
I guess, one application can be to be able to visualize this information using visualization software. Some visualization tool allows to color atoms according the values located in pdb at the position of the temperature factor. But I never used this option my self. Maybe other users know more.

Best
Alessandra

Thanks Alevilla. That’s helpful. So you are saying the -cf option makes a heat map of forces. I guess that makes sense.