<div dir="ltr">Dear Juerg,<div>Thank you very much for your message.</div><div>This may explain the reason my system behaves abnormally (i.e. water molecules in the QM region interact strongly with the MM region to the extend that they form clusters near the QM boundary). </div><div>I have one last question:</div><div>Is it meaningful to use the scc/a1-a2.spl files and add the SUBSYS/KIND/DFTB3_PARAM with the DIAGONAL_DFTB3 and HB_SR_PARAM=4.00?</div><div><br></div><div><br></div><div>Thanks again,</div><div>Saleh</div><div><br></div></div>