A simulation box can look like a setup detail, but choosing it poorly can affect both the reliability and the cost of a molecular dynamics run. Make it too small, and a solute may interact with its own periodic copies. Make it unnecessarily large, and you may spend computation on solvent that does not help answer your question. The useful goal is to choose a box shape and size that protect the simulation from unwanted image interactions without adding avoidable volume.
Why periodic images matter
With periodic boundary conditions, the simulation treats the box as part of a repeating array. When a particle leaves one side, its periodic copy enters from the other. For short-range non-bonded interactions, GROMACS uses the minimum image convention: it considers the nearest image of each particle. This makes periodic simulations practical, but it does not mean the solute can be placed arbitrarily close to its copies. If a solute interacts with one of its own periodic images, the resulting forces may be wrong.
A practical rule of thumb when preparing a solvated system is to leave at least 1.0 nm between the solute and the box boundary. That gives at least 2.0 nm between periodic images of the solute. Treat this as a setup guideline, then check the cutoff restrictions relevant to your simulation in the GROMACS manual. The box should be sized with the actual system and simulation settings in mind—not just made visually comfortable around the molecule.
Choose a shape that fits the solute
Cubic, orthorhombic, and triclinic cells are available, as are rhombic dodecahedron and truncated octahedron cells. For an approximately spherical macromolecule in solution, the latter two can be more efficient: they are closer to a sphere than a cube and need fewer solvent molecules for a given minimum distance between macromolecular images.
The rhombic dodecahedron is the smallest and most regular space-filling unit cell. At the same image distance, its volume is 71% of a cube’s—about 29% less volume. The documentation notes that this can save about 29% of CPU time when simulating a spherical or flexible molecule in solvent. This is a reason to consider the shape, not a guarantee of a particular speedup for every system or workload.
Two ways to size the box
In GROMACS Wizard, the initial fit can be based on Box lengths or Solute-box distance. With Box lengths, the box fits tightly to the system, so increase its size as needed to preserve the recommended clearance. 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, which is useful when you want one box size for the batch.
With Solute-box distance, you 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, each conformation or path frame can have a different box size. Consider whether consistent dimensions or a fit to each structure better suits your comparison.
A quick preparation check
- Pick a cell shape that makes sense for the solute; consider a rhombic dodecahedron or truncated octahedron for an approximately spherical molecule.
- Check the solute-to-boundary clearance rather than judging the box by appearance alone.
- Keep the minimum image convention and your simulation’s cutoff restrictions in view when selecting dimensions.
- For batches, decide whether a shared box or conformation-specific boxes are more appropriate.
One detail that can be surprising: GROMACS keeps particles in a brick-shaped volume for efficiency. When importing GROMACS results, SAMSON tries to detect the unit-cell type, and the importer dialog lets you modify it if needed.
For box setup steps and further details, see the GROMACS Wizard guide to periodic boundary conditions.
Note: SAMSON and all SAMSON Extensions are free for non-commercial use. Get SAMSON.