Dear All,<div><br></div><div>   Is there a particular grid option that will work in general (irrespective of computational cost) for heavy elements? I want to avoid converting bases in Gaussian format to run a Gaussian calculation to check if I am getting correct energies.</div>
<div><br></div><div> Anirban:</div><div>    Did you have to try different combinations for the grid options to get it work? Thanks.</div><div><br></div><div>   Thanks.<br><br><div class="gmail_quote">On Tue, Sep 6, 2011 at 12:02 AM, Peter Knowles <span dir="ltr"><<a href="mailto:KnowlesPJ@cardiff.ac.uk">KnowlesPJ@cardiff.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Dear Anirban,<br>
<br>
I am sorry you have had this trouble with the default grid, which exists in some form for all elements heavier than krypton in the current version of Molpro.  It turned out that a code redesign has been necessary to resolve the problem, and this was completed some time ago, and will appear in the next release of the code.  I have been considering whether the change could be put back into 2010.1, but have concluded that it would not be very safe to do so, so for the moment I hope you can get by with forcing the grid in the way that you have suggested.<br>

<br>
Peter<br>
<div><div></div><div class="h5"><br>
On 1 Sep 2011, at 19:19, Kirk Peterson wrote:<br>
<br>
> Dear Anirban,<br>
><br>
> great, glad this is resolved.<br>
><br>
> -Kirk<br>
><br>
><br>
> On Sep 1, 2011, at 11:18 AM, anirban mandal wrote:<br>
><br>
>> Dear Kirk,<br>
>><br>
>> Indeed it was due to different grid. I switched the grid parameters, and now I get the same energy as obtained with Gaussian and Gamess. Use the following grid along with uks,b-lyp.<br>
>><br>
>> grid;<br>
>> radial,em,2,1,5,5,5,5,;<br>
>> angular,lebedev,0.0,0.0;<br>
>> lmin,29,47,53,53;<br>
>> lmax,53,53,53,53;<br>
>><br>
>> That should give the same energy as with gamess and gaussian.<br>
>><br>
>> Best regards,<br>
>><br>
>> Anirban.<br>
>><br>
>> Anirban Mandal<br>
>> Department of Chemistry<br>
>> Michigan State University<br>
>> East Lansing<br>
>> MI48824-1322<br>
>> From: Kirk Peterson <<a href="mailto:kipeters@wsu.edu">kipeters@wsu.edu</a>><br>
>> To: Gerald Knizia <<a href="mailto:knizia@theochem.uni-stuttgart.de">knizia@theochem.uni-stuttgart.de</a>><br>
>> Cc: anirban mandal <<a href="mailto:anirbanmandal2002@yahoo.co.in">anirbanmandal2002@yahoo.co.in</a>>; molpro-user molpro-user <<a href="mailto:molpro-user@molpro.net">molpro-user@molpro.net</a>><br>

>> Sent: Thursday, 1 September 2011 11:27 AM<br>
>> Subject: Re: [molpro-user] Problem with simple DFT calculation: Ru atom<br>
>><br>
>> The states look identical to me, at least in terms of occupation and orbital character.  Besides since the "bad" energy is much lower, that would imply that both Gaussian and Gamess get the wrong state, which doesn't seem too likely.  So it seems that maybe the changes in the radial grids is the answer.  Do you know if there is a way to reproduce the new grids with 2010.1 ?<br>

>><br>
>> -Kirk<br>
>><br>
>> On Sep 1, 2011, at 3:27 AM, Gerald Knizia wrote:<br>
>><br>
>> > On 01.09.2011 02:21, Kirk Peterson wrote:<br>
>> >> Dear Anirban,<br>
>> >><br>
>> >> actually at least on the surface it looks like a bug to me in 2010.1.  I get essentially the same value as Gaussian/Gamess with the developer's version but I get your incorrect value with the production one (2010.1).  Perhaps someone involved with the DFT code can comment?<br>

>> ><br>
>> > There were some changes to the radial grids in the core region, and this can have significant effects on total energies. Also, could it be possible that both versions converge to different states? There were some changes to the occupations of degenerate orbitals in the initial guess.<br>

>> > --<br>
>> > Gerald Knizia<br>
>><br>
>><br>
>><br>
><br>
</div></div><div class="im">> _______________________________________________<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" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br>
<br>
</div><font color="#888888">--<br>
Prof. Peter J. Knowles<br>
School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, UK<br>
Email <a href="mailto:KnowlesPJ@Cardiff.ac.uk">KnowlesPJ@Cardiff.ac.uk</a><br>
WWW <a href="http://www.cardiff.ac.uk/chemy/contactsandpeople/academicstaff/knowles.html" target="_blank">http://www.cardiff.ac.uk/chemy/contactsandpeople/academicstaff/knowles.html</a><br>
Secretary: Alison Thompson, <a href="mailto:ThompsonAF@Cardiff.ac.uk">ThompsonAF@Cardiff.ac.uk</a> Telephone <a href="tel:%2B44%2029%20208%2074805" value="+442920874805">+44 29 208 74805</a><br>
</font><div><div></div><div class="h5"><br>
<br>
<br>
<br>
_______________________________________________<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" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Regards,<br>Neeraj.<br><br><br>
</div>