<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">By comparing MP2 all-electron calculations performed using localized basis sets and free boundary conditions, with DFT pseudopotential calcs, done with plane waves and periodic boundary conditions you have succeeded in making the two calculations as heterogeneous as possible. Well done (converged) calculations should not depend on the basis set, pseudopotential, or boundary conditions, but they certainly will on the level of theory you use. Getting very different results with MP2 and DFT should make you suspect that some kind of inter-molecular correlation effects may play an important role. The first thing I would do is to repeat the calculation with localized basis sets at the DFT level, and see if the structure you obtain matches the one you get with QE. SB<div> <br><div><div>On Jun 20, 2011, at 6:54 AM, <a href="mailto:psavita@crlindia.com">psavita@crlindia.com</a> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><font face="Default Sans Serif,Verdana,Arial,Helvetica,sans-serif" size="2"><div>Dear Prof. Ghosh,<br><br>As per your suggestion, I performed the structure relaxation of LiH--H2O system using a bigger cubic box of dimension 50 a.u.<br>and I still find H2 abstraction instead of an intramolecular dihydrogen bond.<br>Do you think the box is still inadequate?<br>For such systems, would it be better to put some constraints initially?<br>How do we do this in QE?<br><br>Thanks for your suggestion and best regards,<br><br></div><div><font face="Default Sans Serif,Verdana,Arial,Helvetica,sans-serif" size="2"><div><font style="font-family: Georgia,Default Serif,serif; font-weight: bold;" color="#0000ff"><font color="#999900">Savita Pundlik</font></font><br style="font-family: Georgia,Default Serif,serif;"><span style="font-family: Georgia,Default Serif,serif;"><font><font color="#999900">Computational Materials </font></font></span><font color="#999900">Applied Research Group</font><br style="font-family: Georgia,Default Serif,serif;"><span style="font-family: Georgia,Default Serif,serif; font-weight: bold;"><font><font color="#999900">Computational Research Laboratories Ltd.,</font></font></span><br style="font-family: Georgia,Default Serif,serif;"><span style="font-style: italic; font-family: Georgia,Default Serif,serif;"><font><font color="#999900">Taco House, Damle Path, Off Law College Road</font></font></span><br style="font-style: italic; font-family: Georgia,Default Serif,serif;"><span style="font-style: italic; font-family: Georgia,Default Serif,serif;"><font><font color="#999900">Pune - 411004, India.</font></font></span><br><br></div></font><br><div><br></div><font color="#990099"><a href="mailto:-----pw_forum-bounces@pwscf.org">-----pw_forum-bounces@pwscf.org</a> wrote: -----<br><br></font><blockquote style="padding-right: 0px; padding-left: 5px; margin-left: 5px; border-left: 2px solid #000000; margin-right: 0px;"><font face="Courier New,Courier,monospace" size="3"><br><br><br><br><br>Dear Savita,<br><br>I am not an expert on these type of systems.......however just a small piece<br>of advice.....a single molecule of water has a huge dipole moment.....so a<br>box size of 25 bohr may not be sufficient, until and unless binding to LiH<br>quenches the dipole moment......so you may like to check whether your box is<br>sufficiently big.....usually for a single molecule of water one uses about<br>35-40 bohr box or can use dipole correction.<br><br>Prasenjit<br><br>On 17 June 2011 11:42, <a class="moz-txt-link-rfc2396E" href="mailto:psavita@crlindia.com"><psavita@crlindia.com></a> wrote:<br><br>> Hello,<br>><br>> I have been trying to get the optimum structure of LiH..H2O that has an<br>> intramolecular dihydrogen bond, using pwscf.<br>> I could get the proper structure earlier using mp2 theory and 6-31++g<br>> basis.<br>><br>> For pwscf, the species is enclosed in a cubic box of 25 bohr for gamma<br>> point calculation.<br>> I tried the pz lda, pbe gga as well we blyp PS for the atoms in this<br>> system, however, the first and the last<br>> lead to H2O..LiH, the latter LiH being linear while pbe gives LiH + H2, but<br>> no dihydrogen bond. Is there any other PS<br>> that may lead to proper structure? Or is it that plane wave approach does<br>> not suit systems of this kind?<br>><br>> Thanks in advance and best regards,<br>><br>> Savita Pundlik<br>> Computational Materials Applied Research Group<br>> Computational Research Laboratories Ltd.,<br>> Taco House, Damle Path, Off Law College Road<br>> Pune - 411004, India.<br>><br>><br>><br>> _______________________________________________<br>> Pw_forum mailing list<br>> <a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>> <a href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>><br>><br><br><br>-- <br>PRASENJIT GHOSH,<br>Assistant Professor,<br>IISER Pune,<br>First floor, Central Tower, Sai Trinity Building<br>Garware Circle, Sutarwadi, Pashan<br>Pune, Maharashtra 411021, India<br><br>Phone: +91 (20) 2590 8203<br>Fax: +91 (20) 2589 9790<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://www.democritos.it/pipermail/pw_forum/attachments/20110617/17d0a776/attachment-0001.htm">http://www.democritos.it/pipermail/pw_forum/attachments/20110617/17d0a776/attachment-0001.htm</a> <br><br>------------------------------<br><br>_______________________________________________<br>Pw_forum mailing list<br><a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a><br><br><br>End of Pw_forum Digest, Vol 48, Issue 45<br>****************************************<br></font>
</blockquote><br></div></font>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-size: 14px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">---</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-size: 14px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Stefano Baroni - SISSA</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">&</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">DEMOCRITOS National Simulation Center - Trieste</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-size: 14px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><a href="http://stefano.baroni.me/">http://stefano.baroni.me</a> [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; font-size: 12px; "><br style="font-size: 12px; "></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; font-size: 12px; "><span class="Apple-style-span" style="font-size: 14px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">La morale est une logique de l'action comme la logique est une morale de la pensée - Jean Piaget</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br></span></font></div></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-size: 14px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Please, if possible, don't</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">send me MS Word or PowerPoint attachments</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-size: 14px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Why? See:</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font><font class="Apple-style-span" color="#0023E9" size="3"><span class="Apple-style-span" style="font-size: 12px; "><a href="http://www.gnu.org/philosophy/no-word-attachments.html">http://www.gnu.org/philosophy/no-word-attachments.html</a></span></font></div></div><div><font class="Apple-style-span" color="#0023E9" size="3"><br></font></div></span><br class="Apple-interchange-newline">
</div>
<br></div></body></html>