gromacs - molecular dynamics simulation suite (the Groningen Machine for Chemical Simulations) is a full-featured suite of programs to perform molecular dynamics simulations - in other words, to simulate the behavior of systems with hundreds to millions of particles, using Newtonian equations of motion. It is primarily used for research on proteins, lipids, and polymers, but can be applied to a wide variety of chemical and biological research questions. The following commands make up the GROMACS suite. Please refer to their individual man pages for further details.
Trajectory analysis
gmx-gangle Calculate anglesgmx-distance Calculate distances between pairs of positionsgmx-freevolume Calculate free volumegmx-sasa Compute solvent accessible surface areagmx-select Print general information about selections
Generating topologies and coordinates
gmx-editconf Edit the box and write subgroupsgmx-protonate Protonate structuresgmx-x2top Generate a primitive topology from coordinatesgmx-solvate Solvate a systemgmx-insert-molecules Insert molecules into existing vacanciesgmx-genconf Multiply a conformation in 'random' orientationsgmx-genion Generate monoatomic ions on energetically favorable positionsgmx-genrestr Generate position restraints or distance restraints for index groupsgmx-pdb2gmx Convert coordinate files to topology and FF-compliant coordinate files
Running a simulation
gmx-grompp Make a run input filegmx-mdrun Perform a simulation, do a normal mode analysis or an energy minimizationgmx-convert-tpr Make a modifed run-input file
Viewing trajectories
gmx-nmtraj Generate a virtual oscillating trajectory from an eigenvectorgmx-view View a trajectory on an X-Windows terminal
Processing energies
gmx-enemat Extract an energy matrix from an energy filegmx-energy Writes energies to xvg files and display averagesgmx-mdrun (Re)calculate energies for trajectory frames with -rerun
Converting files
gmx-editconf Convert and manipulates structure filesgmx-eneconv Convert energy filesgmx-sigeps Convert c6/12 or c6/cn combinations to and from sigma/epsilongmx-trjcat Concatenate trajectory filesgmx-trjconv Convert and manipulates trajectory filesgmx-xpm2ps Convert XPM (XPixelMap) matrices to postscript or XPM
Tools
gmx-analyze Analyze data setsgmx-dyndom Interpolate and extrapolate structure rotationsgmx-filter Frequency filter trajectories, useful for making smooth moviesgmx-lie Estimate free energy from linear combinationsgmx-morph Interpolate linearly between conformationsgmx-pme_error Estimate the error of using PME with a given input filegmx-sham Compute free energies or other histograms from histogramsgmx-spatial Calculate the spatial distribution functiongmx-traj Plot x, v, f, box, temperature and rotational energy from trajectoriesgmx-tune_pme Time mdrun as a function of PME ranks to optimize settingsgmx-wham Perform weighted histogram analysis after umbrella samplinggmx-check Check and compare filesgmx-dump Make binary files human readablegmx-make_ndx Make index filesgmx-mk_angndx Generate index files for 'gmx angle'gmx-trjorder Order molecules according to their distance to a groupgmx-xpm2ps Convert XPM (XPixelMap) matrices to postscript or XPM
Distances between structures
gmx-cluster Cluster structuresgmx-confrms Fit two structures and calculates the RMSDgmx-rms Calculate RMSDs with a reference structure and RMSD matricesgmx-rmsf Calculate atomic fluctuations
Distances in structures over time
gmx-mindist Calculate the minimum distance between two groupsgmx-mdmat Calculate residue contact mapsgmx-polystat Calculate static properties of polymersgmx-rmsdist Calculate atom pair distances averaged with power -2, -3 or -6
Mass distribution properties over time
gmx-gyrate Calculate the radius of gyrationgmx-msd Calculates mean square displacementsgmx-polystat Calculate static properties of polymersgmx-rdf Calculate radial distribution functionsgmx-rotacf Calculate the rotational correlation function for moleculesgmx-rotmat Plot the rotation matrix for fitting to a reference structuregmx-sans Compute small angle neutron scattering spectragmx-saxs Compute small angle X-ray scattering spectragmx-traj Plot x, v, f, box, temperature and rotational energy from trajectoriesgmx-vanhove Compute Van Hove displacement and correlation functions
Analyzing bonded interactions
gmx-angle Calculate distributions and correlations for angles and dihedralsgmx-mk_angndx Generate index files for 'gmx angle'
Structural properties
gmx-anadock Cluster structures from Autodock runsgmx-bundle Analyze bundles of axes, e.g., helicesgmx-clustsize Calculate size distributions of atomic clustersgmx-disre Analyze distance restraintsgmx-hbond Compute and analyze hydrogen bondsgmx-order Compute the order parameter per atom for carbon tailsgmx-principal Calculate principal axes of inertia for a group of atomsgmx-rdf Calculate radial distribution functionsgmx-saltbr Compute salt bridgesgmx-sorient Analyze solvent orientation around solutesgmx-spol Analyze solvent dipole orientation and polarization around solutes
Kinetic properties
gmx-bar Calculate free energy difference estimates through Bennett's acceptance ratiogmx-current Calculate dielectric constants and current autocorrelation functiongmx-dos Analyze density of states and properties based on thatgmx-dyecoupl Extract dye dynamics from trajectoriesgmx-principal Calculate principal axes of inertia for a group of atomsgmx-tcaf Calculate viscosities of liquidsgmx-traj Plot x, v, f, box, temperature and rotational energy from trajectoriesgmx-vanhove Compute Van Hove displacement and correlation functionsgmx-velacc Calculate velocity autocorrelation functions
Electrostatic properties
gmx-current Calculate dielectric constants and current autocorrelation functiongmx-dielectric Calculate frequency dependent dielectric constantsgmx-dipoles Compute the total dipole plus fluctuationsgmx-potential Calculate the electrostatic potential across the boxgmx-spol Analyze solvent dipole orientation and polarization around solutesgmx-genion Generate monoatomic ions on energetically favorable positions
Protein-specific analysis
gmx-do_dssp Assign secondary structure and calculate solvent accessible surface areagmx-chi Calculate everything you want to know about chi and other dihedralsgmx-helix Calculate basic properties of alpha helicesgmx-helixorient Calculate local pitch/bending/rotation/orientation inside helicesgmx-rama Compute Ramachandran plotsgmx-wheel Plot helical wheels
Interfaces
gmx-bundle Analyze bundles of axes, e.g., helicesgmx-density Calculate the density of the systemgmx-densmap Calculate 2D planar or axial-radial density mapsgmx-densorder Calculate surface fluctuationsgmx-h2order Compute the orientation of water moleculesgmx-hydorder Compute tetrahedrality parameters around a given atomgmx-order Compute the order parameter per atom for carbon tailsgmx-potential Calculate the electrostatic potential across the box
Covariance analysis
gmx-anaeig Analyze the eigenvectorsgmx-covar Calculate and diagonalize the covariance matrixgmx-make_edi Generate input files for essential dynamics sampling
Normal modes
gmx-anaeig Analyze the normal modesgmx-nmeig Diagonalize the Hessian for normal mode analysisgmx-nmtraj Generate a virtual oscillating trajectory from an eigenvectorgmx-nmens Generate an ensemble of structures from the normal modesgmx-grompp Make a run input filegmx-mdrun Find a potential energy minimum and calculate the Hessian
Consult the manual at <http://www.gromacs.org/content/view/27/42/> for an introduction to molecular dynamics in general and GROMACS in particular, as well as an overview of the individual programs. The shorter HTML reference is available in /usr/share/doc/gromacs/html/ . The development of GROMACS is mainly funded by academic research grants. To help us fund development, the authors humbly ask that you cite the GROMACS papers: H.J.C. Berendsen, D. van der Spoel and R. van Drunen. GROMACS: A message-passing parallel molecular dynamics implementation. Comp. Phys. Comm. 91, 43-56 (1995) Erik Lindahl, Berk Hess and David van der Spoel. GROMACS 3.0: A package for molecular simulation and trajectory analysis. J. Mol. Mod. 7, 306-317 (2001) B. Hess, C. Kutzner, D. van der Spoel, and E. Lindahl. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 4, 3, 435-447 (2008), <http://dx.doi.org/10.1021/ct700301q> Current developers: David van der Spoel <[email protected]>
Berk Hess <[email protected]>
Erik Lindahl <[email protected]> A full list of present and former contributors is available at <http://www.gromacs.org> This manual page is largely based on the GROMACS online reference, and was prepared in this format by Nicholas Breen <[email protected]>. GROMACS has no major known bugs, but be warned that it stresses your CPU more than most software. Systems with slightly flaky hardware may prove unreliable while running heavy-duty simulations. If at all possible, please try to reproduce bugs on another machine before reporting them.