<div dir="ltr">Thank you.<br><br>On Friday, October 27, 2017 at 11:28:56 AM UTC+2, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<div><br></div><div>I did make an initial implementation. It is more or less working / passed some simple tests, but hasn't really been used, documented or extensively tested to date. </div><div><br></div><div>You are welcome to use it, but please be even more careful than normal, and report any strange results / crashes.</div><div><br></div><div>It should just be a case of adding shell parameters to atom types within the MM section to active it.</div><div><br></div><div>Matt<br><br>On Thursday, October 26, 2017 at 1:49:41 PM UTC+1, <a>liuxi...@gmail.com</a> wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Wondering if shell-core model is implemented in QM/MM in CP2K 4.1 and 5.<br><br>Regards,<br>Liu<br><br>On Tuesday, March 31, 2015 at 7:13:15 PM UTC+2, Matt W wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi,</div><div><br></div>Unfortunately (see below) the core-shell model isn't implemented for QMMM in CP2K at the moment (shell charges are ignored, if I remember correctly, and you get what you report).<div><br></div><div>I had told Teo (Teodoro Laino) that I would implement this, and do have a more-or-less working core-shell QMMM implementation somewhere. But it would take some effort to bring back into line with a current version of CP2K and debug properly. It would be plausible to work on this with a bit of help if you were interested, but would take a bit of time to get working.</div><div><br></div><div>Matt<br><br>On Tuesday, March 31, 2015 at 3:28:54 PM UTC+1, Chris Lee wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><font face="arial,sans-serif">Dear CP2K users,<br><br>I'm trying to model a hydroxylated iron oxide (Fe2O3) surface using QM/MM. The top half of the slab is modelled using QM, while the classical system is modelled using a core-shell model for the oxygen atoms only.<br><br>The issue that I'm having is that the system is supposed to be electrically neutral but I'm getting an overall negative charge as the QM oxygen atoms all have a charge of -2 but the QM iron and hydrogen atoms have zero charge. If I disable the core-shell part of the model none of the QM atoms show any charge (I think it's supposed to be this way, please correct me if I'm wrong) and the classical system is overall neutral with each atom showing the correct charge.<br><br>Any help with this would be greatly appreciated as I'd very much like to use a core-shell model for my classical system.<br><br>I've attached the input files that I'm using.<br><br>Chris<br></font></div></blockquote></div></div></blockquote></div></blockquote></div></div></blockquote></div>