I am checking the interaction between two molecules, but when I am trying to optimise, the bonds are breaking, and the original structures are deforming. I want to optimise while retaining their initial structure. There must be a command like geom=connectivity to do such things. Please help.