Though not yet implemented, the virtual site approach is applicable also to OLIVES and could be developed if there is a demand. Another way is to turn domain decomposition off, likely with a computational cost depending on hardware architecture. At our local high-performance computing (HPC) facility at CSCAA, 52 we often do not use domain decomposition given our specific hardware architecture of 1 GPU and 18 CPUs, however, this will be different on other HPC clusters where running domain decomposition could be more efficient. Ideally, GROMACS would implement an LJ bonded type that could be defined between particle pairs but included in the non-bonded list upon topology generation.