<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Amit,<div class=""><br class=""></div><div class="">as far as I know, one doesn’t have to treat ECP-based calculations any differently compared to those without. It is a little surprising you get the same CP corrected binding energies with these two basis sets, but DZ results can be strange sometimes.</div><div class=""><br class=""></div><div class="">regards,</div><div class=""><br class=""></div><div class="">-Kirk</div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On May 6, 2016, at 5:21 AM, Amit Sharma <<a href="mailto:greifsw@gmail.com" class="">greifsw@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div dir="ltr" class="">Hi all, <div class="">I want to calculate BSSE where one of the atoms has an ECP potential basis on it. Should BSSE be computed differently for ECP? There is some discussion in the literature about switching off the core but I am not sure how to do that in molpro. </div><div class=""><br class=""></div><div class="">As a test I compute BSSE for two noble gas atoms with AVDZ basis and get a BSSE of about 47 cm-1. Then I switch the basis on one of the atoms to be def2-SVPD and I get BSSE which is twice that of the AVDZ result (about 97 cm-1). </div><div class=""><br class=""></div><div class="">However, the CP corrected interaction energy is approximately the same between the two calculations. 400 cm-1 vs 403 cm-1, so maybe nothing special needs to be done. Can anyone please help clarify this for me.</div><div class="">Thanks</div><div class="">Amit</div><div class=""><br class=""></div><div class=""><div style="margin: 0px; line-height: normal;" class=""><span style="font-family:Menlo;font-size:13px;color:rgb(206,121,36)" class="">***</span><span style="font-family:Menlo;font-size:13px" class="">, Ar </span><span style="font-family:Menlo;font-size:13px;color:rgb(206,121,36)" class="">+</span><font face="Menlo" class=""> he bsse </font></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">geometry</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">{Ar</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">He Ar Rx}</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">Rx</span><span style="color:rgb(206,121,36)" class="">=5.0</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">!basis</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">avdz,he</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">avdz</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">basis</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">def2-SVPD,he</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">avdz</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> text,calculation </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> complex</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> hf !RHF </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> total system</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> ccsd</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">t</span><span style="color:rgb(213,59,211)" class="">)</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> e_Arhe</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">energy !save energy in variable e_Arhe</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> text,cp calculation </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> HE</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> dummy,Ar !make Ar a dummy</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> hf !RHF </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> He</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> ccsd</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">t</span><span style="color:rgb(213,59,211)" class="">)</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> e_he</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">energy !save energy in variable e_he</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> text,cp calculation </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> Ar</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> dummy,he !make he dummy</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> hf</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> ccsd</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">t</span><span style="color:rgb(213,59,211)" class="">)</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> e_Ar</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">energy !save energy in variable e_Ar</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> text,separate calculation </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> he</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> symmetry,x,y</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> geometry</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">{he} !geometry </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> he alone</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> hf !RHF </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> He</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> ccsd</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">t</span><span style="color:rgb(213,59,211)" class="">)</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> e_he_inf</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">energy !save energy in variable e_he_inf</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> text,separate calculation </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> Ar</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> symmetry,x,y</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> geometry</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">{Ar} !geometry </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> Ar alone</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> hf !scf </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> Ar</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> ccsd</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">t</span><span style="color:rgb(213,59,211)" class="">)</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class=""> e_Ar_inf</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">energy !save energy in variable e_Ar_inf</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo; min-height: 15px;" class=""><span style="" class=""></span><br class=""></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">de</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">i</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">=</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">e_Arhe-e_he_inf-e_Ar_inf</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">*</span><span style="" class="">tocm !compute uncorrected interaction energy</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">de_cp</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">i</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">=</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">e_Arhe-e_he-e_Ar</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">*</span><span style="" class="">tocm !compute counter-poise corrected interaction energy</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">bsse_he</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">i</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">=</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">e_he-e_he_inf</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">*</span><span style="" class="">tocm !BSSE </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> he</span></div><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">bsse_Ar</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">i</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">=</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">e_Ar-e_Ar_inf</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">*</span><span style="" class="">tocm !BSSE </span><span style="color:rgb(206,121,36)" class="">for</span><span style="" class=""> Ar</span></div><p style="margin:0px;font-size:13px;line-height:normal;font-family:Menlo" class="">
</p><div style="margin: 0px; font-size: 13px; line-height: normal; font-family: Menlo;" class=""><span style="" class="">bsse_tot</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">i</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">=</span><span style="" class="">bsse_he</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">i</span><span style="color:rgb(213,59,211)" class="">)</span><span style="color:rgb(206,121,36)" class="">+</span><span style="" class="">bsse_Ar</span><span style="color:rgb(213,59,211)" class="">(</span><span style="" class="">i</span><span style="color:rgb(213,59,211)" class="">)</span><span style="" class=""> !total BSSE</span></div></div></div>
_______________________________________________<br class="">Molpro-user mailing list<br class=""><a href="mailto:Molpro-user@molpro.net" class="">Molpro-user@molpro.net</a><br class="">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.molpro.net_mailman_listinfo_molpro-2Duser&d=CwICAg&c=C3yme8gMkxg_ihJNXS06ZyWk4EJm8LdrrvxQb-Je7sw&r=mqaXWsZ0ghQR-pXvPLjY_g&m=B4tbS47XX15ILr4ZVAykT-Xbolq5pVnAkv_J6VroMFA&s=F11UmXf2GReqdN-sR8iHnnhZnWxhUD-cTwkh4cU8yJo&e= </div></blockquote></div><br class=""></div></body></html>