Hi, I am doing non equilibrium alchemical transformations. When annihilating a ligand from the protein pocket with rlist=2.2 and couple-intramol=no everything works fine. But when creating the same ligand in a box of equilibrated water with rlist=2.2 and couple-intramol=no I get:
There are 1 perturbed non-bonded pair interactions beyond the pair-list cutoff
of 2.2 nm, which is not supported. This can happen because the system is
unstable or because intra-molecular interactions at long distances are
excluded. If the latter is the case, you can try to increase nstlist or rlist
to avoid this.The error is likely triggered by the use of couple-intramol=no
and the maximal distance in the decoupled molecule exceeding rlist.
The problem is that I am sure that the ligand perfectly fits in rlist and rlist is actually big enough to mostly fit the whole water box (4 X 4 X 4 nm circa).
I have tried to proceed the creation with different speeds but it didn’t change
And the problem doesn’t happen 100% of the time and often if I re-run the transformation it arrives to the end without problems.
Does someone have an idea/solution for the problem? (without changing couple-intramol)
ps I am creating vdw and q separately and it happens in both situations