algorithms used in |Gromacs|.
Before version 4.6, |Gromacs| always used pair-lists based on groups of
-particles. These groups of particles were orginally charge-groups, which were
+particles. These groups of particles were originally charge-groups, which were
necessary with plain cut-off electrostatics. With the use of PME (or
reaction-field with a buffer), charge groups are no longer necessary
(and are ignored in the Verlet scheme). In |Gromacs| 4.6 and later, the
with bio-molecular systems dominated by water.
Without PME, the group cut-off scheme should generally be combined
-with a buffered pair-list to help avoid artefacts. However, the
+with a buffered pair-list to help avoid artifacts. However, the
group-scheme kernels that can implement this are much slower than
either the unbuffered group-scheme kernels, or the buffered
Verlet-scheme kernels. Use of the Verlet scheme is strongly encouraged
exact cut-off shift/switch always
potential-shift interactions yes yes
potential-switch interactions yes yes
-force-switch interations yes yes
+force-switch interactions yes yes
switched potential yes yes
switched forces yes yes
non-periodic systems yes Z + walls
slower. The performance of the Verlet scheme with the new non-bonded
kernels is independent of system composition and is intended to always
run with a buffered pair-list. Typically, buffer size is 0 to 10% of
-the cut-off, so you could win a bit of peformance by reducing or
+the cut-off, so you could win a bit of performance by reducing or
removing the buffer, but this might not be a good trade-off of
simulation quality.
by being set in your shell to any non-NULL value. Specific
requirements are described below if other values need to be set. You
should consult the documentation for your shell for instructions on
-how to set environment variables in the current shell, or in config
+how to set environment variables in the current shell, or in configuration
files for future shells. Note that requirements for exporting
environment variables to jobs run under batch control systems vary and
you should consult your local documentation for details.
energies as generated during your :ref:`gmx mdrun`.
The file can be transformed to a portable energy file (portable
-accross hardware platforms), the :ref:`edr` file using the program
+across hardware platforms), the :ref:`edr` file using the program
:ref:`gmx eneconv`.
See also :ref:`gmx energy`.
itp
---
-The itp file extension stands for include toplogy. These files are included in
+The itp file extension stands for include topology. These files are included in
topology files (with the :ref:`top` extension).
.. _log:
databank file format describes the positions of atoms in a molecular
structure. Coordinates are read from the ATOM and HETATM records,
until the file ends or an ENDMDL record is encountered.
-GROMACS programs can read and write a simlation box in the
+GROMACS programs can read and write a simulation box in the
CRYST1 entry.
The pdb format can also be used as a trajectory format:
-several structures, seperated by ENDMDL, can be read from
+several structures, separated by ENDMDL, can be read from
or written to one file.
Example
+++++++
-An pdb file should look like this::
+A pdb file should look like this::
ATOM 1 H1 LYS 1 14.260 6.590 34.480 1.00 0.00
ATOM 2 H2 LYS 1 13.760 5.000 34.340 1.00 0.00
rtp
---
-The rtp file extension stands for residue toplogy.
+The rtp file extension stands for residue topology.
Such a file is needed by :ref:`gmx pdb2gmx`
to make a GROMACS topology for a protein contained in a :ref:`pdb`
file. The file contains the default interaction type for the 4 bonded
:ref:`gmx pdb2gmx` sets the number exclusions to 3, which
means that interactions between atoms connected by at most 3 bonds are
excluded. Pair interactions are generated for all pairs of atoms which are
-seperated by 3 bonds (except pairs of hydrogens).
+separated by 3 bonds (except pairs of hydrogens).
When more interactions need to be excluded, or some pair interactions should
not be generated, an ``[exclusions]`` field can be added, followed by
-pairs of atom names on seperate lines. All non-bonded and pair interactions
+pairs of atom names on separate lines. All non-bonded and pair interactions
between these atoms will be excluded.
A sample is included below.
:ref:`gmx trjconv`). It is easier to do a correct "single-point" energy
evaluation with this feature than a 0-step simulation.
-Neighbour searching is normally performed for every frame in the
+Neighbor searching is normally performed for every frame in the
trajectory, since :ref:`gmx mdrun` can no longer assume anything about how the
structures were generated. If :mdp:`nstlist` is zero, then only one
-neighbour list will be constructed. Naturally, no update or constraint
+neighbor list will be constructed. Naturally, no update or constraint
algorithms are ever used.
Running a simulation in reproducible mode
more control is useful. `mdrun -nsteps 100` overrides the [.mdp] file
and executes 100 steps. `mdrun -maxh 2.5` will terminate the
simulation shortly before 2.5 hours elapse, which can be useful when
-running under cluster queues (as long as the queueing system does not
+running under cluster queues (as long as the queuing system does not
ever suspend the simulation).
Running a membrane protein embedding simulation
with a procedure similar to genbox, but since lipids are much larger
than water molecules it will lead to a large vacuum layer between the
protein and membrane if we remove all molecules where any atom is
-overlapping. Instead, this module works by first artifically shrinking
+overlapping. Instead, this module works by first artificially shrinking
the protein in the xy-plane, then it removes lipids that overlap with
a much smaller core, after which we gradually push the protein atoms
back to their initial positions, while using normal dynamics for the
To use membrane embedding, start by building a lipid bilayer that is
just-so-slightly larger in the xy-plane than what you expect to need
in the end, and make sure you have enough water outside the membrane
-to accomodate globular domains. Place the protein in the same coordinate
+to accommodate globular domains. Place the protein in the same coordinate
file (and topology) as your lipid bilayer, and make sure it is in the
orientation and position you want right in the middle of the bilayer.
and configurable intelligence to detect your hardware and make pretty
effective use of that hardware. For a lot of casual and serious use of
:ref:`gmx mdrun`, the automatic machinery works well enough. But to get the
-most from your hardware to maximise your scientific quality, read on!
+most from your hardware to maximize your scientific quality, read on!
Hardware background information
-------------------------------
inform the user about the choices made and possible consequences.
A number of command-line parameters are available to vary the default
-behaviour.
+behavior.
``-nt``
The total number of threads to use. The default, 0, will start as
``-ntomp``
The total number of OpenMP threads per rank to start. The
default, 0, will start one thread on each available core.
- Alternatively, mdrun will honour the appropriate system
+ Alternatively, mdrun will honor the appropriate system
environment variable (e.g. ``OMP_NUM_THREADS``) if set.
``-npme``
comes for free with the non-bonded communication. Particles beyond
the non-bonded cut-off are only communicated when they have
missing bonded interactions; this means that the extra cost is
- minor and nearly indepedent of the value of ``-rdd``. With dynamic
+ minor and nearly independent of the value of ``-rdd``. With dynamic
load balancing, option ``-rdd`` also sets the lower limit for the
domain decomposition cell sizes. By default ``-rdd`` is determined
by :ref:`gmx mdrun` based on the initial coordinates. The chosen value will