I am trying to run a constrained walk simulation for committor analysis of transition path ensembles. I constrain the QM atoms as well as some of the residues interacting with it using "cons harm force" and run the dynamics, and this calculation starts to become erratic within few steps and crashes. I have attached the input and the output files here. I would appreciate any help that can make me understand this crash and fix it. Thanks!