<div dir="ltr">1) Yes<div>2) Yes</div><div>3) Yes </div><div><br></div><div>I attach to the mail my three inputs. Thank you, again, for your help.</div><div><br></div><div><div style="font-family:arial,sans-serif;font-size:13px;color:rgb(136,136,136)">Valentina Cantatore</div><div style="font-family:arial,sans-serif;font-size:13px;color:rgb(136,136,136)">PostDoc Researcher @ Università del Piemonte Orientale </div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2014-09-18 11:26 GMT+02:00 Giovanni Cantele <span dir="ltr"><<a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank">giovanni.cantele@spin.cnr.it</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5"><br>
On 18 Sep 2014, at 10:40, Valentina Cantatore <<a href="mailto:valecantatore@gmail.com">valecantatore@gmail.com</a>> wrote:<br>
<br>
> Good morning to everyone.<br>
><br>
> I would like to calculate the interaction energy between to molecules using QE. What I set up was to fix the dimension of a cubic box at 20 A (I think it is big enough to avoid interactions between the other replicated molecules), calculate the energy of each of the two fragments (Ea and Eb) and then calculate the energy of the two fragments together (Eab). What I expect is that interaction energy (Ein) is given by Ein=Eab-(Ea+Eb).<br>
><br>
> I also used the assume_isolated option for all calculations. I sampled only gamma points.<br>
><br>
> The results obtained are not of the same order of magnitude of those obtained with a gaussian basis set on the same systems.<br>
><br>
> Where am I wrong? Is it possible run this kind of calculation with Quantum Espresso? Thank you<br>
><br>
> Valentina Cantatore<br>
> PostDoc Researcher @ Università del Piemonte Orientale<br>
><br>
<br>
<br>
</div></div>the procedure you describe is correct and this kind of calculations is of course possible with QE. Just few questions:<br>
1) are you using the same XC approximation with gaussian and plane waves runs (e.g. PBE)?<br>
2) are you sure that both calculations are converged with respect to the respective parameters (in particular basis set size)?<br>
3) did you carefully check the energy units used by the two different codes?<br>
<br>
Otherwise, you should provide more details and/or inputs to understand what is going on.<br>
<br>
Giovanni<br>
<br>
--<br>
<br>
Giovanni Cantele, PhD<br>
CNR-SPIN<br>
c/o Dipartimento di Fisica<br>
Universita' di Napoli "Federico II"<br>
Complesso Universitario M. S. Angelo - Ed. 6<br>
Via Cintia, I-80126, Napoli, Italy<br>
e-mail: <a href="mailto:giovanni.cantele@spin.cnr.it">giovanni.cantele@spin.cnr.it</a><br>
Phone: <a href="tel:%2B39%20081%20676910" value="+39081676910">+39 081 676910</a><br>
Skype contact: giocan74<br>
<br>
ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>
Web page: <a href="http://people.na.infn.it/~cantele" target="_blank">http://people.na.infn.it/~cantele</a><br>
<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br></div>