<div>Dear cp2k users/developers,</div><div><br></div><div>I'm trying to simulate guanine DNA guanine quadruplex using QMMM approach. MM force field parameters are prepared with AmberTools20. QM part is constituted of Na+ cation coordinated by 8 guanine molecules (so called G-octet), while other nucleobases, counterions and solvating water molecules are treated with MM.</div><div><br></div><div>System was previously equlibrated with NVT ensemble and MM approach with no problem encountered. Also, when I try to run QM part of the system only using same DFT settings as in the case of QMMM calculation I have no problems.</div><div><br></div><div>However, when I try to run QMMM calculation under present setup (all input files in the attachment) I constantly get shift of the Cons Qty during my run, for both NVE and NVT runs (graphs from my last attempt are in the attachment - last 500fs).</div><div><br></div><div>I tried different calculation setup (different timestep, ewald summation parms etc.) by looking to the previous topics here in the group but with no success.</div><div><br></div><div>Maybe there is something trivial bothering me here, qmmm link or similar...</div><div>I used IMOMM method for the link and also distributed charges across neighbouring MM atoms to neutralize MM subsystem.</div><div><br></div><div>Can anybody help me with this?</div><div><br></div><div>Best,</div><div>Bane</div>