<div dir="ltr">Just a follow-up to my previous post.<div><br></div><div>Based on a number of runs, and looking at the number of CSF's that come out, I've convinced myself that the double excitations are from the CLOSED-CORE orbitals to all virtuals.</div>
<div><br></div><div>The definition of "active space" in Molpro's CASPT2 (RS2) program seems to just be specific to the RS2 program itself.</div><div><br></div><div>The definition of active space for RS2 and RS2C is laid out in slightly clearer language in the "quickstart" portion of the online manual. <a href="http://www.molpro.net/info/2012.1/doc/quickstart/node24.html">http://www.molpro.net/info/2012.1/doc/quickstart/node24.html</a> </div>
<div>Quoted below in case the online manual changes in the future</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<p style="color:rgb(0,0,0);font-family:Times;font-size:medium">In case of <tt>rs2</tt> and <tt>rs3</tt>, exactly the same configuration spaces as in the MRCI are used. In this case the excitations with two electrons in the external orbital space are <em>internally contracted</em>. The total number of correlated orbitals is restricted to 16 for machines with 32-bit integers and to 32 for machines with 64-bit integers.</p>
<p style="color:rgb(0,0,0);font-family:Times;font-size:medium">In the <tt>rs2c</tt> case certain additional configuration classes involving internal and semi-internal excitations are also internally contracted (see J. Chem. Phys. <b>112</b>, 5546 (2000)). This is exactly as in the <tt>mrcic</tt> case (see above). This method is much more efficient than <tt>rs2</tt> and more suitable for large cases. In particular, in this case only the number of <em>active</em> orbitals is restricted to 16 or 32 on 32 and 64 bit machines, respectively, and any number of closed-shell (inactive) orbitals can be used (up to a maximum as defined by a program parameter).</p>
</blockquote><div>I've pretty much exhausted my questions on the limitations/capabilities of Molpro's CASPT2 program.<br></div><div><br></div><div>Final question: Why is there a limit of 32 correlated orbitals in RS2? I fiddled with the source code a bit, and it doesn't seem to be just an arbitrary number set when the program was written - but the reason is not quite transparent from a brief look at the code in src/mrci/.</div>
<div><br></div><div><br></div><div class="gmail_extra"><br>
<br><div class="gmail_quote">On Wed, Jul 2, 2014 at 12:51 PM, Alan Chien <span dir="ltr"><<a href="mailto:alandc@umich.edu" target="_blank">alandc@umich.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr">Thank you for your reply Werner, its cleared up a lot of my initial questions. <div>I had actually determined that there was a difference between CORE and FROZEN, but was still trying to figure things out a bit better.<div>
<br></div><div>I have a few clarification questions. Please forgive the length, but I'm trying to be as precise as possible so that no further clarifications will be necessary, heh.</div>
<div><br></div><div>For the CAS portion of CASPT2 (separate from the MULTI calc), the reference determinants are generated using the orbitals defined by OCC - CLOSED. </div><div>For the PT2 portion, we have to define a space from which we can do doubles excitations, i.e. correlate certain orbitals. This is defined as those orbitals in CLOSED - CORE. <u>But how do we define where those double excitations go to</u>? Do they only go to the empty orbitals that we included in OCC? Or are they automatically allowed to be excited to any/all virtuals (Which seems unlikely since virtuals usually make up a large portion of the MOs). <br>
<div><br></div><div>CASPT2 (RS2) defines "active space" as both the space used to create the CAS determinants <u>as well as</u> the orbitals which we allow to be correlated. </div><div>I find this strange as before now, I've only ever worked with CAS/RAS, and I'm used to thinking of active spaces as just the orbital set used to generate the set of reference determinants. Is this definition of active space a quirk of Molpro or is it the accepted standard when discussing CASPT2?</div>
<div><br></div></div></div><div>Finally, does RS2C treat the CORE,CLOSED,OCC cards differently when making its definition of the "active space"? This link (which I posted earlier), <a href="http://www.molpro.net/pipermail/molpro-user/2006-April/001759.html" style="font-size:13px;font-family:arial,sans-serif" target="_blank">http://www.molpro.net/pipermail/molpro-user/2006-April/001759.html</a>, seems to suggest that it does. But I am unclear as to why that may be as the manual states that the RS2 and RS2C programs are essentially the same, just one is faster or has different features implemented compared to the other. Perhaps it is just a difference in the code of each program instead of a difference having any substantive meaning?</div>
<div><br></div><div>I understand that some of these questions I can answer myself (namely the last one) with a couple of Molpro test jobs, and in lieu of provided answers, I may well determine them myself.</div><div><br>
</div>
<div>But any help is welcome,</div><div>Alan Chien</div><div><br></div></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jul 2, 2014 at 6:12 AM, Werner Győrffy <span dir="ltr"><<a href="mailto:gyorffy@theochem.uni-stuttgart.de" target="_blank">gyorffy@theochem.uni-stuttgart.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Dear Alan,<br>
<br>
CLOSED orbitals are doubly occupied in all reference configurations but excitations are allowed from them in the CASPT2 calculations if they are not frozen. If you want to define orbitals which are frozen (not correlated) in CASPT2, you should use CORE orbitals. CORE orbitals are not counted for the 32 orbital limitation in the rs2 program.<br>
<br>
Active orbitals in CASSCF are defined as OCC minus CLOSED orbitals.<br>
<br>
An (X)MS-CASPT2 calculation is still possible only with the rs2 program not with rs2c.<br>
<br>
Regards,<br>
<br>
Werner Gyorffy.<div><br>
<br>
<br>
<br>
On 07/01/2014 11:56 PM, Alan Chien wrote:<br>
</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div>
Interesting nugget of information. The example given in the multi-state<br>
CASPT2 section of the manual DOES change its active space when you<br>
change the numbers in the closed card.<br>
So it seems that for that particular example it can recognize closed input.<br>
<br>
Example: <a href="http://www.molpro.net/info/current/examples/lif_mr_mscaspt2.com" target="_blank">http://www.molpro.net/info/<u></u>current/examples/lif_mr_<u></u>mscaspt2.com</a><br>
(I edited it so that it only did r=3 ang to reduce the amt of output I<br>
had to look at. Also removed everything below RS2 section)<br>
<br>
Trying to replicate it now. In fact, been trying for a while.<br>
<br>
<br>
<br>
<br>
On Tue, Jul 1, 2014 at 3:48 PM, Alan Chien <<a href="mailto:alandc@umich.edu" target="_blank">alandc@umich.edu</a><br></div><div><div>
<mailto:<a href="mailto:alandc@umich.edu" target="_blank">alandc@umich.edu</a>>> wrote:<br>
<br>
Dear Molpro Community,<br>
<br>
I was wondering if there was something wrong with my following input<br>
for an extended multi-state CASPT2 calculation. My input is not<br>
recognizing the CLOSED card, and so is throwing the error "TOO MANY<br>
ACTIVE ORBITALS: 48 THIS VERSION ALLOWS 32." An old thread mentions<br>
that I should use RS2C -<br>
<a href="http://www.molpro.net/pipermail/molpro-user/2006-April/001759.html" target="_blank">http://www.molpro.net/<u></u>pipermail/molpro-user/2006-<u></u>April/001759.html</a> -<br>
which should recognize closed orbitals.<br>
<br>
But it's a rather old thread, and I want XMS-CASPT2 which isn't<br>
implemented for RSC2, and so I thought it worth mentioning to see if<br>
they had updated the RS2 code and perhaps my input was just off. The<br>
manual makes no mention of RS2 being unable to recognize the CLOSED<br>
card. I also tried the same input with FROZEN and that threw the<br>
same error.<br>
<br>
Any input welcome,<br>
Alan Chien<br>
<br>
<br>
***,title<br>
memory,400,m<br>
symmetry<br>
geometry={<br>
S 1.9692714536 0.9837700026 0.0000000000<br>
C 0.6745300892 -0.1941173367 0.0000000000<br>
C 1.1941024748 -1.4906632561 0.0000000000<br>
C 2.5920412491 -1.5265017306 0.0000000000<br>
C 3.1677860531 -0.2809638443 0.0000000000<br>
C -0.6745300892 0.1941173367 0.0000000000<br>
C -1.1941024748 1.4906632561 0.0000000000<br>
C -2.5920412491 1.5265017306 0.0000000000<br>
C -3.1677860531 0.2809638443 0.0000000000<br>
S -1.9692714536 -0.9837700026 0.0000000000<br>
C -4.5811338904 -0.0064991780 0.0000000000<br>
C -5.5229897499 1.0602944736 0.0000000000<br></div></div>
N -6.<a href="tel:2538051640" value="+12538051640" target="_blank">2538051640</a> <tel:<a href="tel:2538051640" value="+12538051640" target="_blank">2538051640</a>> 1.9321608642<div><br>
0.0000000000<br>
C 4.5811338904 0.0064991780 0.0000000000<br>
C 5.5229897499 -1.0602944736 0.0000000000<br></div>
N 6.<a href="tel:2538051640" value="+12538051640" target="_blank">2538051640</a> <tel:<a href="tel:2538051640" value="+12538051640" target="_blank">2538051640</a>> -1.9321608642<div><br>
0.0000000000<br>
C 5.0766433381 1.3388671764 0.0000000000<br></div>
N 5.4513787354 2.<a href="tel:4133206959" value="+14133206959" target="_blank">4133206959</a> <tel:<a href="tel:4133206959" value="+14133206959" target="_blank">4133206959</a>><div><br>
0.0000000000<br>
C -5.0766433381 -1.3388671764 0.0000000000<br></div>
N -5.4513787354 -2.<a href="tel:4133206959" value="+14133206959" target="_blank">4133206959</a> <tel:<a href="tel:4133206959" value="+14133206959" target="_blank">4133206959</a>><br>
0.0000000000<br>
H -3.1636666579 2.<a href="tel:4343293653" value="+14343293653" target="_blank">4343293653</a> <tel:<a href="tel:4343293653" value="+14343293653" target="_blank">4343293653</a>><br>
0.0000000000<br>
H -0.<a href="tel:5803308458" value="+15803308458" target="_blank">5803308458</a> <tel:<a href="tel:5803308458" value="+15803308458" target="_blank">5803308458</a>> 2.3697711304<br>
0.0000000000<br>
H 0.<a href="tel:5803308458" value="+15803308458" target="_blank">5803308458</a> <tel:<a href="tel:5803308458" value="+15803308458" target="_blank">5803308458</a>> -2.3697711304<br>
0.0000000000<br>
H 3.1636666579 -2.<a href="tel:4343293653" value="+14343293653" target="_blank">4343293653</a> <tel:<a href="tel:4343293653" value="+14343293653" target="_blank">4343293653</a>><div><br>
0.0000000000<br>
}<br>
<br>
basis=6-31G*<br>
hf<br>
{casscf;<br>
closed,31,6,30,5 !core orbs, total 72<br>
occ,31,8,30,7 !active orbs, total 76<br>
wf,148,1,0;state,2<br>
wf,148,3,0<br>
pspace,100.0<br>
}<br>
<br>
{rs2,xms=1,mix=2;<br>
state,2,2,3<br>
}<br>
<br>
<br>
<br>
<br>
<br></div>
______________________________<u></u>_________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net" 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/<u></u>listinfo/molpro-user</a><br>
<br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net" 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/<u></u>listinfo/molpro-user</a><br>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div></div>