Does minimizeEnergy work across the PBC?
Posted: Wed May 09, 2018 6:57 pm
Hi all,
I have a pre-equilibrated cell, to which I added a ligand that dips over the edge of the cell, forming bad contacts with water.
minimizeEnergy() is unable to remove the contacts - does this function work across the boundary? It can be fixed by increasing the z direction of the PBC (and recovering with a barostat if necessary), just wondering how it works.
I have a pre-equilibrated cell, to which I added a ligand that dips over the edge of the cell, forming bad contacts with water.
minimizeEnergy() is unable to remove the contacts - does this function work across the boundary? It can be fixed by increasing the z direction of the PBC (and recovering with a barostat if necessary), just wondering how it works.
Code: Select all
initial energy = 2135329600067.072 kJ/mol
after energy = 2298638365251.072 kJ/mol