<div dir="ltr">Hi everyone, <div><br></div><div>@Matt: Thank you very much for the clarification, it was helpful.</div><div><br></div><div>Since full QM/MM periodicity is not possible, I now would like to get the QM/MM simulation of water working without QM/MM periodicity, but with electrostatic embedding (E_COUPL = COULOMB). The QM method is still DFTB3. </div><div>And I have tried almost everything to get it working, but my water molecules still fall apart. </div><div>I have tried...</div><div><ul><li>different GMAX values in the EWALD section</li><li>different ALPHA values in the EWALD section</li><li>different O_SPLINE values in the EWALD section</li><li>the QM part completely without electrostatics (no POISSON section)</li><li>DFTB3 with and without dispersion correction</li><li>different time step sizes</li><li>Standard BOMD and second generation CPMD</li><li>...</li></ul></div><div>Nothing seems to really solve the problem. I have also tried more things such as mechanical embedding with and without QM/MM periodicity. </div><div><br></div><div>Here is a link to a sample input package which can be started directly, all files included:</div><div><br></div><div>https://drive.google.com/open?id=0B-Ij17NhUNZoeUpETFJhS1hpbFE<br></div><div><br></div><div>Any help I would greatly appreciate. </div><div><br></div><div>Many thanks,</div><div>Jadzia</div><div><br></div><div><br><br>On Thursday, September 7, 2017 at 12:26:24 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">Hi again,<div><br></div><div>sorry I don't really have time to give long enough replies to be clear, so perhaps a bit confusing.</div><div><br></div><div>QM-QM coupling is not needed in your case - your QM box is the same size as the QMMM box. It calculates multipole moments of the QM charge distribution and subtracts artificial interactions if the QM box is smaller than QMMM and adds back in correct ones.</div><div><br></div><div>QM-MM long range interactions as you say are the about including the ewald summed potential from the classical part - I think that this is not possible except using E_COUPL = GAUSS, which in turn requires DFT. See the paper from Laio et al for details.</div><div><br></div><div>There are several different QMMM methods available, and they have different ranges of applicability.</div><div><br></div><div>Matt<br><br>On Thursday, September 7, 2017 at 10:28:58 AM UTC+1, Jadzia wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thank you for your answer Matt. <div><br></div><div>Since error messages can be very helpful if they are precise, I would like summarize my understanding so far related to the error message I have gotten:<div><span style="color:rgb(32,18,77);font-family:monospace,monospace;font-size:x-small"><br></span></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><span style="color:rgb(32,18,77);font-family:monospace,monospace;font-size:x-small">QM-QM long range correction not possible with COULOMB coupling.</span></div></blockquote><div><br></div><div>Firstly, from what you wrote I implicitly conclude that the periodic section which I have added to my input file went into the direction of what you had in mind when you wrote about including the coupling of the QM region with the periodic MM images. In this case the above error message which I got seems to be not fully suitable since I had disabled the QM-QM long range correction (i.e. the MULTIPOLE section). To my understanding (which is admittedly incomplete) the problem might be of a more general nature and summarized as follows: <span style="color:rgb(32,18,77);font-family:monospace,monospace;font-size:x-small"><i>QM/MM long range correction not possible with COULOMB coupling</i></span></div><div><font size="2" face="arial, sans-serif" color="#000000">because it seems to cover both cases QM-MM and optionally also QM-QM long term corrections (if the multipole section is not switched off) for the QM region. </font><br></div><div><div><br></div><div>Secondly, the error message seems to be too general in another aspect, because you wrote that QM/MM periodic couplings are possible with DFT (and Coulomb coupling I guess), in which case the error message does not hold in its current general form. </div><div><br></div><div><br></div><div><br></div><div><br></div><div><br><div><br>On Thursday, September 7, 2017 at 9:26:55 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">As I said, I wasn't sure if there was a periodic correction implemented using Coulomb coupling for DFTB/semi-empirical.<div><br></div><div>Looks like you just have to run without periodic coupling, or switch to DFT to use periodic coupling.</div><div><br>On Thursday, September 7, 2017 at 4:41:31 AM UTC+1, Jadzia wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Another thing I wonder about: Matt you said that I should keep the section QMMM periodic section so that the QM region feels the potential of the periodically replicated MM region. Which I would call a "MM-QM long range correction". But the error message of my last post states that a QM-QM long range correction is not possible with Coulomb coupling when I add this section. Therefore I wonder if the periodic section which I set up for the MM-QM long range correction was indeed what you had in mind...<br></div></blockquote></div></div></blockquote></div></div></div></div></div></blockquote></div></div></blockquote></div></div>