<HTML><BODY>Dear Kirk Peterson.<br>Is it possible to combine MDF ecps and def2 basis sets?<br>Of course, we can write such a program in molpro, but is it acceptable?<br>Or mixing of ecps and basis sets is a rather controversial idea, which we should throw away immediately?<br> Sasha Medvedev.<br><br><br>Пятница, 25 октября 2013, 21:02 UTC от "Peterson, Kirk" <kipeters@wsu.edu>:<br>
<blockquote style="border-left:1px solid #0857A6; margin:10px; padding:0 0 0 10px;" class="mailru-blockquote">
<div id="">
<div class="js-helper js-readmsg-msg">
<style type="text/css"></style>
<div id="style_13827349440000000828" class="mr_read__body">
<base target="_self" href="https://e.mail.ru/">
<div id="style_13827349440000000828_BODY">
<div><br>
</div>
I can’t say that the dhf- versions are any more reliable, but certainly the new ECP’s on which they are based are more accurate than the older MWB ones that the def2 sets are matched with.
<div><br>
</div>
<div>best,</div>
<div><br>
</div>
<div>-Kirk</div>
<div><br>
<div>
<div>On Oct 25, 2013, at 1:51 PM, sasha medvedev <<a href="sentmsg?mailto=mailto%3avinsanity305@mail.ru" target="_blank">vinsanity305@mail.ru</a>> wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<div>Dear Kirk Peterson.<br>
I decided that the cardinal number n (i.e. d,t,q,5,6) of a basis set should be equal to lmax. I thought, that lmax and n are identical.<br>
Thank you for explanation of my misunderstanding, concerning the construction of correlation consistent basis sets.<br>
----------------------<br>
> <em>yes, there was an update to this file by Gerald Knizia in August of 2012. Presumably this is included in the latest production version of Molpro. In this new file most of the previous def2- sets that used the older MWB ECPs have been modified to use
the new MDF ones. These are called dhf-</em><br>
<br>
So, dhf basis sets are more reliable? (in comparison with def2 sets)<br>
Sasha Medvedev<br>
<br>
Пятница, 25 октября 2013, 15:36 UTC от "Peterson, Kirk" <<a href="sentmsg?mailto=mailto%3akipeters@wsu.edu" target="_blank">kipeters@wsu.edu</a>>:<br>
<blockquote style="border-left:1px solid #0857A6;margin:10px;padding:0 0 0 10px;" class="mailru-blockquote">
<div id="">
<div class="js-helper js-readmsg-msg">
<div id="style_13827154150000000610" class="mr_read__body">
<div id="style_13827154150000000610_BODY">Dear Sasha,
<div><br>
</div>
<div>unfortunately what you describe below about increasing the lmax of light atoms so that they would match heavier ones would spoil the very thing that makes them systematically convergent towards the CBS limit, i.e., with each “n”, functions are added that
contribute similar amounts of correlation energy. The higher angular momentum functions (f for DZ, g for TZ, etc.) are not included since they do not “belong” energetically at that level. One could of course construct sets in this way, but it’s not clear
what their convergence properties would be. And then if you wanted to do molecules containing f-block elements, you have to increase them again….</div>
<div><br>
</div>
<div>best regards,</div>
<div><br>
</div>
<div>-Kirk</div>
<div><br>
<div>
<div>On Oct 24, 2013, at 2:24 PM, sasha medvedev <<a target="_blank">vinsanity305@mail.ru</a>> wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<div>Dear Mr. Grant Hill.<br>
Thank you for a reply. <br>
<br>
>The def2 series of basis sets was not designed to systematically approach the basis set limit...<br>
<br>
These piece of information is important for us.<br>
----------------------------------------------<br>
>Think of a minimal basis for a first row element (C) and one for a transition metal (Ti),
<br>
lmax must be different simply to construct the occupied orbitals.<br>
<br>
I get your idea about Carbon and Titanium minimal atomic basis sets. But I think, that it is possible<br>
to make the correlation consistent basis sets in such a way, that all elements have the same lmax for a particular basis set.<br>
For instance, one should augment the minimal Carbon basis set (lmax=p=1) with d-functions. As the result, Titanium and Carbon will<br>
have the same value of lmax (for a particular basis set). This approach is fruitful until we take into consideration lanthanide and actinide elements
<br>
(in such case we need to go further and add f-functions to C and Ti basis sets).<br>
----------------------------------------------<br>
>One approach when using correlation consistent sets is to keep lmax as the largest for the system (usually the heaviest element) - meaning cardinal number + 1 in this case.<br>
<br>
So, in case of AB molecule we should compare lmax for A and B and then choose the bigger one for extrapolation formula?
<br>
Hm, it sounds sensible. Thanks.<br>
<br>
Alexander Medvedev<br>
<br>
On, 23 Oct 2013, 16:32 +01:00 от Grant Hill <<a target="_blank">Grant.Hill@glasgow.ac.uk</a>>:<br>
<blockquote style="border-left:1px solid #0857A6;margin:10px;padding:0 0 0 10px;" class="mailru-blockquote">
<div id="">
<div class="js-helper js-readmsg-msg">
<div id="style_13825423310000000902" class="mr_read__body">
<div id="style_13825423310000000902_BODY">Dear Sasha,<br>
<br>
The def2 series of basis sets was not designed to systematically approach the basis set limit, so the standard extrapolation formulae are not really expected to work here (although I see no reason why Molpro should stop you trying - GIGO). There may be some
more appropriate formulae in the literature, but I am not aware of them.<br>
<br>
<br>
On 23 Oct 2013, at 13:17, sasha medvedev <<a target="_blank">vinsanity305@mail.ru</a>> wrote:<br>
<br>
> Besides that, extrapolation formulae depend on the cardinal number n (which is the part of the basis set name - i.e. 2 for asvp, 3 for atzvp, 4 for aqzvp, and should be equal to the lmax). When we looked through output file, we noticed, that the same basis
set has different lmax for different atoms:<br>
<br>
This is entirely expected. Think of a minimal basis for a first row element (C) and one for a transition metal (Ti), lmax must be different simply to construct the occupied orbitals. One approach when using correlation consistent sets is to keep lmax as the
largest for the system (usually the heaviest element) - meaning cardinal number + 1 in this case. A lack of high accuracy reference data means this isn't particularly well tested, but it _seems_ to work.<br>
<br>
Best regards,<br>
<br>
Grant<br>
<br>
</div>
</div>
</div>
</div>
</blockquote>
<br>
<br>
-- <br>
s m<br>
</div>
_______________________________________________<br>
Molpro-user mailing list<br>
<a target="_blank">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 data-mce-bogus="1">
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<br>
<br>
-- <br>
s m<br>
</div>
</blockquote>
</div>
<br>
</div>
</div>
<base target="_self" href="https://e.mail.ru/">
</div>
</div>
</div>
</blockquote>
<br>
<br>-- <br>s m<br></BODY></HTML>