<pre>Dear Molpro users,<br><br> I am trying to calculate some valence and the low-lying Rydberg states of<br>formaldehyde. This is part of my input:<br><br>-----------------------------------<br>multi;<br>occ,10,4,4,1;closed,2,0,0,0;<br>
config;<br>restrict,0,1,7.1,8.1,9.1,10.1,3.2,4.2,3.3,4.3,1.4;<br>wf,16,1,0;state,5;<br>wf,16,2,0;state,3;<br>wf,16,3,0;state,4;<br>wf,16,4,0;state,4;<br>orbital,2140.2<br>natorb, ci, print=10;<br>------------------------------------<br>
<br>The Rydberg states are described as single excitations to the virtual<br>orbitals given in the "restrict" part but, with this input, I m also<br>getting high-lying Rydberg states that I am not interested in.<br>
<br>I just want to calculate the low-lying Rydberg states coming from single<br>excitations from three of the occupied orbitals (1b1, 1b2 and 2b2), that<br>is,<br>the following configurations:<br><br>(3a1, 4a1, 5a1)^6 (1b1)^1 (1b2, 2b2)^4 (7a1-10a1, 3b1-4b1, 3b2-4b2, 1a2)^1<br>
<br>and<br><br>(3a1, 4a1, 5a1)^6 (1b1)^2 (1b2, 2b2)^3 (7a1-10a1, 3b1-4b1, 3b2-4b2, 1a2)^1<br><br>So, my question is: How can I further restrict the calculation to include<br>
just excitations from the 1b1, 1b2 and 2b2 occupied orbitals to the above-mentioned virtual orbitals?<br><br>If someone could help me, I would be very thankful.<br><br>Thanks a lot and best regards,<br><font color="#888888">Susana</font></pre>