The invisible neighbor: choosing a safer GROMACS simulation box

A simulation box can look perfectly reasonable on screen and still create a problem you cannot see: under periodic boundary conditions, the system is surrounded by copies of itself. If the box is too tight, a solute may interact with its own periodic image. That can lead to forces that do not represent the isolated system you intended to model.

Choosing a box is therefore not just a matter of fitting the molecule into a container. You need a shape that suits the solute and enough separation between periodic images. Here is a practical way to think through both decisions when preparing a system in GROMACS Wizard.

First, picture the repeating system

Periodic boundary conditions make the simulation cell repeat in space. For short-range non-bonded interactions, GROMACS uses the minimum image convention: only the nearest image of a particle is considered. The practical concern is that the solute should not interact with its own periodic copy. Long-range electrostatics need additional treatment, such as lattice-sum methods, so the minimum-image rule is not a substitute for choosing appropriate electrostatics settings.

A useful starting rule is to leave at least 1.0 nm between the solute and the box boundary. This gives at least 2.0 nm between periodic images of the solute. Treat this as a rule of thumb, not a universal guarantee: the cut-off restrictions and simulation settings for your system still matter.

Choose a shape that fits the solute

GROMACS Wizard supports cubic, orthorhombic, triclinic, rhombic dodecahedron, and truncated octahedron unit cells. Cubic boxes are familiar, but a cube is not always the most economical fit. For an approximately spherical macromolecule in solution, a rhombic dodecahedron or truncated octahedron is closer to a sphere. With the same minimum distance between macromolecular images, these shapes require fewer solvent molecules than a cube.

The rhombic dodecahedron is the smallest and most regular space-filling unit cell. The documentation reports that its volume is 71% of a cube’s volume at the same image distance, which can save about 29% of CPU time when simulating a spherical or flexible molecule in solvent. That comparison is about the specified geometry and conditions; it is not a promise of the same speed-up for every system.

For an elongated or irregular solute, a different cell shape may fit the system more naturally. The key is not to choose the most compact-looking shape automatically, but to consider the solute’s overall shape and ensure the required separation in every relevant direction.

Pick the fitting method that matches your workflow

When setting up the box, GROMACS Wizard offers two useful ways to define its initial size:

• Box lengths: Set the box dimensions directly. The initial fit is tight, so increase the dimensions as needed to maintain sufficient distance from periodic images. This option also lets you move the box to position the system. For a batch project, the initial size is based on all conformations or the whole path, making it useful when you want one consistent box size across the batch.
• Solute-box distance: Specify the distance from the solute to the box. At least 1 nm is recommended as a practical way to support the minimum image convention. In a batch project, box sizes can differ between conformations or frames.

In other words, use box lengths when consistent dimensions or manual positioning are important. Use solute-box distance when you want the padding around the system to guide the initial fit. In either case, check the resulting dimensions rather than assuming that selecting a shape alone has solved the separation problem.

A quick preparation check

1. Look at the solute’s shape: is it approximately spherical, or does it extend much farther in one direction?
2. Choose a suitable supported unit-cell shape, considering a rhombic dodecahedron or truncated octahedron for an approximately spherical solute in solvent.
3. Set the box using lengths or solute-box distance, depending on whether you need consistent batch dimensions or padding-based fitting.
4. Check that the solute has adequate clearance from the boundary. Use the 1.0 nm guideline as a starting point, then consider the GROMACS cut-off restrictions and settings for your simulation.

One detail that can be confusing later: GROMACS keeps particles in a brick-shaped volume for efficiency. When importing GROMACS results, SAMSON attempts to detect the unit-cell type, and the type can also be modified in the trajectory importer dialog.

For the complete explanation and links to the related preparation steps, see the GROMACS Wizard guide to periodic boundary conditions.

Note: SAMSON and all SAMSON Extensions are free for non-commercial use. Get SAMSON.