<div dir="ltr">Hello every Developer and Researcher,<div>I tried the example in https://www.cp2k.org/exercises:2017_ethz_mmm:qmmm several days ago. It is a geo_opt example of 48 KCl atoms using QMMM method. And I also look up the paper listed in the tutorial and tried to calculate the band gap of the system. But the result is not satisfactory. </div><div>If I use 1-layer QM atoms (16 atoms), the HOMO - LUMO gap is as following</div><div><div class="prettyprint" style="background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; word-wrap: break-word;"><code class="prettyprint"><div class="subprettyprint"><div class="subprettyprint">HOMO - LUMO gap [eV] :    4.926845</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.668970</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.734728</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.728571</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.733307</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.692899</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.722199</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.681768</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.678310</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.700719</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662065</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.685438</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.663563</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662155</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662057</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662044</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662065</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662064</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662034</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662058</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662063</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.662052</div></div></code></div><br>If I use 2-layer QM atoms (32 atoms), the result is </div><div><div class="prettyprint" style="background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; word-wrap: break-word;"><code class="prettyprint"><div class="subprettyprint"><div class="subprettyprint">HOMO - LUMO gap [eV] :    4.782295</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.271173</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.447092</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.470585</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.499691</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.487309</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.495278</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.486996</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481340</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481181</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481379</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481309</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481298</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481304</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481308</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481297</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.481311</div></div></code></div><div><br></div>If I use full QM,</div><div><div class="prettyprint" style="background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; word-wrap: break-word;"><code class="prettyprint"><div class="subprettyprint"><div class="subprettyprint">HOMO - LUMO gap [eV] :    4.934751</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.683587</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.923171</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.941706</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.934531</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.924194</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.916435</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.917274</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.915809</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.916783</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.913765</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.916681</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.914778</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.916037</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.918605</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.916562</div><div class="subprettyprint"> HOMO - LUMO gap [eV] :    4.916525</div></div></code></div><br>The question is the 2-layer QM atoms should be more approaching to the result of full QM than the 1-layer atoms. But the fact is the opposite. The setting about printing the HOMO-LUMO gap is as follows</div><div><div class="prettyprint" style="background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; word-wrap: break-word;"><code class="prettyprint"><div class="subprettyprint"><div class="subprettyprint">&MO_CUBES</div><div class="subprettyprint">        NHOMO 1</div><div class="subprettyprint">        NLUMO 1</div><div class="subprettyprint">        WRITE_CUBE F</div><div class="subprettyprint">&END MO_CUBES</div></div></code></div><div><br></div>The most possible reason I can think of is something wrong appear in the setting of full QM calculation.</div><div>Thanks for your reply ! <br><br></div><div><br><br></div></div>