A molecular dynamics system can look ready for production after its temperature has stabilized, yet still be changing in volume. That is the job of NPT equilibration: give the system time to settle at the chosen temperature and pressure, then check whether its density has reached a stable regime. One common source of confusion is that pressure can fluctuate widely even while density is becoming steady. Treating every pressure swing as a sign of failure can lead to unnecessary changes—or a premature stop.
This guide focuses on a practical question: How can you tell whether NPT equilibration has done enough?
Start from the equilibrated system
NPT follows energy minimization and NVT equilibration. Begin with the GRO file or batch project from the previous step; GROMACS Wizard can fill in the path to the previous successful run using its auto-fill button. Keeping the same temperature-coupling choices used during NVT helps maintain continuity. If velocity generation is enabled, check that its temperature is consistent too.
Judge progress by the density trend
The aim is for density to settle into a stable regime. Look at the running average of density over time: does it approach a plateau, or does it continue to drift? Density can still fluctuate around a target value, so a perfectly flat trace is not the expectation. The useful signal is whether the overall trend has become steady enough for your system.
Pressure is different. It may fluctuate widely during equilibration and simulation. A noisy pressure trace alone does not tell you that the system has failed to equilibrate. Consider the density behavior and the pressure-related terms together, and give the system more time if density has not stabilized.
The tutorial’s example reaches a density of about 1030 kg/m3. That is close to the experimental value of 1000 kg/m3; for the SPC/E model, the expected density is about 1008 kg/m3. These figures illustrate how to interpret that particular example—they are not a universal target for every molecular system. Use a value appropriate to your system and model.
Give it a reasonable first run, then reassess
The time needed depends on the system’s size and contents. About 100 ps is described as a reasonable starting point, not a guarantee that equilibration is complete. After the run, inspect the plots. If density is still drifting or the relevant terms are slow to converge, continue NPT equilibration from the results of the previous NPT run. You do not have to treat the first run as an all-or-nothing test.
In GROMACS Wizard, the NPT tab provides options to generate input files, run locally, or launch a cloud computation. Once the job finishes, its plots are generated and saved with the results. You can import the full trajectory, just the final frame, or a selected range of frames. If you choose not to import the trajectory, cancelling the import does not delete the results or prevent plot generation.
Check the setup before changing parameters
Default NPT parameters are intended for typical equilibration runs, and the advanced settings let you adjust items such as pressure coupling. The tutorial suggests exponential-relaxation pressure coupling with a stochastic term (c-rescale) and a time constant of about 5 ps for most projects. These are guidance values, not a substitute for system-specific judgment. Position-restraint options are filled automatically based on the input system.
If you need to continue a run, first confirm that you are using the intended previous results and that temperature coupling remains consistent with NVT. Then use the density plot to decide whether more time is needed. Once density is stable enough for your system, proceed to production molecular dynamics.
For the complete workflow and interface details, see the GROMACS Wizard NPT Equilibration tutorial.
Note: SAMSON and all SAMSON Extensions are free for non-commercial use. Get SAMSON.