<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Wafaa,<div><br></div><div>probably the lack of tight functions in your basis set via something like cc-pwCVTZ wouldn't affect your RHF convergence. That is probably something more fundamental and I'm surprised you didn't run into any problems when you optimized your geometry in the frozen core approximation. In regards to the difference between cc-pCVTZ and cc-pwCVTZ, these were both developed to recover core correlation effects, but the cc-pCVTZ describes both the core-core and core-valence correlation on equal footing while the cc-pwCVTZ is biased towards the intershell core-valence correlation. They do both converge to the same basis set limit but the latter will give slightly higher energies but generally more accurate core-valence correlation "effects" on bond lengths and relative energetics.</div><div><br></div><div>Yes, the explicitly correlated methods, e.g., UCCSD(T)-F12, are capable of very high accuracy with relatively small basis sets. I'm quite biased, but you should try the cc-pVnZ-F12 series of basis sets (n=D, T, Q). For core correlation the cc-pCVnZ-F12 sets have just become available. The latter are in Molpro 2009.2 already and for 2009.1 you can download the sets (including new auxiliary RI fitting sets) from our website, <a href="http://tyr0.chem.wsu.edu/~kipeters/basis.html">http://tyr0.chem.wsu.edu/~kipeters/basis.html</a>.</div><div><br></div><div>best regards,</div><div><br></div><div>Kirk</div><div><br><div><div>On Feb 28, 2010, at 12:34 PM, Wafaa Fawzy wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div lang="EN-US" link="blue" vlink="purple" style="word-wrap: break-word;
-webkit-nbsp-mode: space;-webkit-line-break: after-white-space">
<div class="Section1"><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Dear
Kirk,</span> <o:p></o:p></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Thank
you so much for your quick response. May be my problem I that I did not change
the basis set, I used aug-cc-pVTZ. I looked up the </span>cc-pCVTZ or
cc-pwCVTZ basis, I found that both of them include 7s,69, 4d, and 3f orbitals
for Oxygen. Do you know what is the difference between the two? Did you try
the explicit correlation methods, e.g. UCCSD(T)-F12? What do you think of this
method?<o:p></o:p></p><p class="MsoNormal">Thanks and Regards, Wafaa <span style="font-size:11.0pt;
font-family:"Calibri","sans-serif""><o:p></o:p></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D"><o:p> </o:p></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in"><p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> Kirk Peterson
[mailto:kipeters@wsu.edu] <br>
<b>Sent:</b> Sunday, February 28, 2010 11:37 AM<br>
<b>To:</b> Wafaa Fawzy<br>
<b>Cc:</b> <a href="mailto:molpro-user@molpro.net">molpro-user@molpro.net</a><br>
<b>Subject:</b> Re: [molpro-user] CCSd(T) with core<o:p></o:p></span></p>
</div>
</div><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">Dear Wafaa,<o:p></o:p></p>
<div><p class="MsoNormal"><o:p> </o:p></p>
</div>
<div><p class="MsoNormal">in principle the core directive as you've shown it only
affects the uccsd(t) calculation and not the rhf step, which is where your job
is failing. So without any other info I would say your optimization has taken a
step to a geometry where the rhf is struggling to converge. Inspection of
your log file (not the out file) may be helpful. Perhaps you could attach
your input file? Note that when correlating the core electrons you should
change your basis set to one that is designed for this purpose, e.g., cc-pCVTZ
or cc-pwCVTZ.<o:p></o:p></p>
</div>
<div><p class="MsoNormal"><o:p> </o:p></p>
</div>
<div><p class="MsoNormal">regards,<o:p></o:p></p>
</div>
<div><p class="MsoNormal"><o:p> </o:p></p>
</div>
<div><p class="MsoNormal">Kirk<o:p></o:p></p>
</div>
<div><p class="MsoNormal"><o:p> </o:p></p>
<div>
<div><p class="MsoNormal">On Feb 27, 2010, at 10:10 AM, Wafaa Fawzy wrote:<o:p></o:p></p>
</div><p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
<div>
<div>
<div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Hi
All,<o:p></o:p></span></p>
</div>
<div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">I
ran optimization calculations that converged without the core electrons being
correlated. I repeated the same job with the option core: “uccsd(t);core”
but the optimization does not converge, I get an error saying : no
convergence, the problem occurs in rhfpro”. <o:p></o:p></span></p>
</div>
<div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">In
my calculations, it is important to get accurate correlation energy.
Does anyone has an explanation or suggestions for solving the problem?<o:p></o:p></span></p>
</div>
<div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Thanks
in advance,<o:p></o:p></span></p>
</div>
<div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Wafaa
Fawzy<o:p></o:p></span></p>
</div>
<div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Dept.
of Chem.<o:p></o:p></span></p>
</div>
<div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">MSU<o:p></o:p></span></p>
</div>
<div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""> <o:p></o:p></span></p>
</div>
</div><p class="MsoNormal"><span style="font-size:13.5pt;font-family:"Helvetica","sans-serif"">_______________________________________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user">http://www.molpro.net/mailman/listinfo/molpro-user</a><o:p></o:p></span></p>
</div>
</div><p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</blockquote></div><br></div></body></html>