<div dir="ltr"><div><div><div><div>Dear Paolo,<br><br></div>Assuming everything was set up correctly (for which I would have to see your input files), the magnitude of the adsorption energy is way too large to be accounted for by dipole corrections. What I would suggest is for you to isolate the TiO2 slab from Anastase_Surface_without_H_anisole_1.out and run an scf calculation on it. I would then see what the difference is with the final energy you obtain from Anastase_Surface_without_H.out. If the difference is very large, there's your answer. This suggestion is triggered by the observation that the reconstruction on the "bottom" surface of the slabs in the two output files appears different. This may not be the answer but worth a shot.<br><br></div>Beyond this, there are a couple things to note about your calculations. The first one is that it appears that you are not using van der Waals corrections. Since your molecules are perpendicular, it shouldn't matter that much but most commonly, such a molecule would adsorb in the parallel configuration, so presumably you will investigate that configuration as well. So it may be worthwhile investigating the effect of some vdW schemes. Secondly, your molecule appears to travel a lot on the slab during the geometry optimization. This is a little unusual, I would investigate this.<br><br></div>Hope this helps.<br></div>Hande<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 25, 2018 at 7:01 PM, Paolo Costa <span dir="ltr"><<a href="mailto:paolo.costa85@gmail.com" target="_blank">paolo.costa85@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span class=""><br clear="all"><div>Dear experts,</div><div><br></div><div>I got nonsense value of adsorption energy by computing 4-iodoanisole on top of TiO2 anastase surface (101). </div><div>By calculating : (E<font size="1">surface+anisole</font>)-E<font size="1">anisole</font>-E<font size="1">s<wbr>urface, </font>I got an adsoprtion energy of 0.66 ry (8.98 eV), which is too high. </div><div>Now I am trying to apply the dipole correction to see if this is the problem causing an unusual adsorption energy. </div></span><div>In attachment I put the output files.</div><span class=""><div>I would be glad to receive any suggestions to understand the problem. </div><div><br></div><div>Thanks.</div><div><br></div></span><span class="HOEnZb"><font color="#888888"><div>Paolo</div></font></span><span class="">-- <br><div class="m_2684118038942672275gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Paolo Costa, Ph.D.<div>Postdoctoral Researcher</div><div>Department of Chemistry and Biomolecular Sciences</div><div>University of Ottawa</div><div>10 Marie Curie, Ottawa, ON K1N 6N5, Canada</div><div>Room number: DRO 326 (D'Iorio Hall)</div></div></div>
</span></div>
<br>______________________________<wbr>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.<wbr>org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-<wbr>espresso.org/mailman/listinfo/<wbr>users</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div>Hande Toffoli<br></div><div>Department of Physics</div><div>Office 439</div><div>Middle East Technical University</div><div>Ankara 06800, Turkey</div><div>Tel : +90 312 210 3264</div><div><a href="https://www.compmatsci.com/" target="_blank">https://www.compmatsci.com/</a><br></div><div><br></div></div></div></div></div></div>
</div>