Dear Gromacs users,<div><br></div><div>Could someone please shed light on the following problem? My system consists of an alkanethiol SAM (in the xy plane) with a layer of water on top. The mdrun command works when the box vectors are 1.1, 1.3, and 2.5, but the SAM flies apart and the simulation crashes. I'm trying to lengthen the box in the z direction (and of course, resolvating the system) to see if that solves the problem. Unfortunately, this seems to give rise to domain decomposition issues... I get the following output in my md.log file:</div>
<div><br></div><div><div>Initial maximum inter charge-group distances:</div><div> two-body bonded interactions: 0.424 nm, LJ-14, atoms 36 39</div><div> multi-body bonded interactions: 0.424 nm, Ryckaert-Bell., atoms 36 39</div>
<div>Minimum cell size due to bonded interactions: 0.467 nm</div><div>Maximum distance for 5 constraints, at 120 deg. angles, all-trans: 0.794 nm</div><div>Estimated maximum distance required for P-LINCS: 0.794 nm</div><div>
This distance will limit the DD cell size, you can override this with -rcon</div><div>Scaling the initial minimum size with 1/0.8 (option -dds) = 1.25</div><div>Optimizing the DD grid for 6 cells with a minimum initial size of 0.992 nm</div>
<div>The maximum allowed number of cells is: X 1 Y 1 Z 5</div></div><div><div>Fatal error:</div><div>There is no domain decomposition for 6 nodes that is compatible with the given box and a minimum cell size of 0.9925 nm</div>
<div>Change the number of nodes or mdrun option -rcon or -dds or your LINCS settings</div></div><div><br></div><div>I've experimented a bit with -rcon and with the number of nodes, but I really have no idea what I'm doing! Any thoughts would be appreciated.</div>
<div>Olivia<br clear="all"><div><br></div><br>
</div>