<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" target="_blank" style="font-size:13px;font-family:arial,sans-serif">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 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:0 0 0 .8ex;border-left:1px #ccc 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 class=""><br>
<br>
<br>
<br>
On 07/01/2014 11:56 PM, Alan Chien wrote:<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="">
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 class="h5">
<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 class=""><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 class=""><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 class=""><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 class=""><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>