GMX-RMSF(1)
π GMX-RMSF(1) β the gmx rmsf page, bestie edition π
okay so. this is the manual page for gmx-rmsf, aka gmx rmsf, from the GROMACS suite (that's the manual this page lives in, mwah). it's the one that calculates atomic fluctuations. yes, the wiggly stuff. the "how much does this atom vibe-shift around its average position" stuff. iconic behavior honestly. π¦
β¨ NAME β¨
gmx-rmsf β calculate atomic fluctuations π
β¨ SYNOPSIS β¨
gmx rmsf [-f [<.xtc/.trr/...>]] [-s [<.tpr/.gro/...>]] [-n [<.ndx>]]
[-q [<.pdb>]] [-oq [<.pdb>]] [-ox [<.pdb>]] [-o [<.xvg>]]
[-od [<.xvg>]] [-oc [<.xvg>]] [-dir [<.log>]] [-b <time>]
[-e <time>] [-dt <time>] [-[no]w] [-xvg <enum>] [-[no]res]
[-[no]aniso] [-[no]fit]
β¨ DESCRIPTION β¨
so gmx rmsf computes the root mean square fluctuation β RMSF, which is literally just the standard deviation, no cap β of atomic positions in the trajectory. that trajectory comes in through -f.
and before it does all that, it optionally fits to a reference frame, which you supply with -s. cute, right?
with the -oq option, the RMSF values get converted into B-factor values, and those get written out to a .pdb file. by default, the coordinates in that output file are taken from the structure file you gave it with -s β although you can instead use coordinates read from a different .pdb file, provided with -q. π
β οΈ okay real talk for a sec: there is very little error checking here. so in this case, it is your responsibility to make sure all atoms in the structure file and the .pdb file correspond exactly to each other. exactly. like, no wiggle room. you've been warned bestie π
the -ox option writes the B-factors to a file together with the average coordinates in the trajectory.
and with the -od option, the root mean square deviation with respect to the reference structure gets calculated. different thing! similar-ish name! don't mix them up π
with the -aniso option, gmx rmsf will compute anisotropic temperature factors, and then it will also output average coordinates and a .pdb file with ANISOU records (that corresponds to the -oq or -ox option). π
β οΈ please note that the U values are orientation-dependent, so before you go comparing anything with experimental data, you should verify that you fit to the experimental coordinates. seriously. verify. πΈ
when a .pdb input file is passed to the program and the -aniso flag is set, a correlation plot of the Uij will be created β but only if any anisotropic temperature factors are present in the .pdb file. conditional queen π
with the -dir option, the average MSF (3x3) matrix gets diagonalized. this shows you the directions in which the atoms fluctuate the most and the least. main character directions, honestly π₯
β¨ OPTIONS β¨
π options to specify input files:
-f [<.xtc/.trr/...>](traj.xtc) β Trajectory:xtctrrcptgrog96pdbtngh5md-s [<.tpr/.gro/...>](topol.tpr) β Structure+mass(db):tprgrog96pdbbrkent-n [<.ndx>](index.ndx) (Optional) β Index file-q [<.pdb>](eiwit.pdb) (Optional) β Protein data bank file
π¦ options to specify output files:
-oq [<.pdb>](bfac.pdb) (Optional) β Protein data bank file-ox [<.pdb>](xaver.pdb) (Optional) β Protein data bank file-o [<.xvg>](rmsf.xvg) β xvgr/xmgr file-od [<.xvg>](rmsdev.xvg) (Optional) β xvgr/xmgr file-oc [<.xvg>](correl.xvg) (Optional) β xvgr/xmgr file-dir [<.log>](rmsf.log) (Optional) β Log file
π other options:
-b <time>(0) β time of first frame to read from the trajectory (default unit isps)-e <time>(0) β time of last frame to read from the trajectory (default unit isps)-dt <time>(0) β only use a frame whent MOD dt = first time(default unit isps)-[no]w(no) β view output.xvg,.xpm,.epsand.pdbfiles-xvg <enum>(xmgrace) β xvg plot formatting:xmgrace,xmgr,none-[no]res(no) β calculate averages for each residue-[no]aniso(no) β compute anisotropic temperature factors-[no]fit(yes) β do a least squares superposition before computing RMSF. without this, you must make sure that the reference structure and the trajectory match π€
β¨ SEE ALSO β¨
more information about GROMACS is available at <http://www.gromacs.org/> π
β¨ COPYRIGHT β¨
2026, GROMACS development team
and that's the whole page, bestie π«Ά this was GROMACS 2026.4, dated Oct 02, 2026 β go forth and fluctuate responsibly, you absolute icon πβ¨
π¦ where this page comes from
- nixpkgs unstable: gromacs 2026.4
- ubuntu stonking: gromacs-data 2026.4-1ubuntu1
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms