Dear everyone,<br><br>  Hello,<br><br>  Excuse me, I have two questions about the frozen core in the SAPT calculation in molpro 2009. <br><br>(i) in the manual,  <a href="http://www.molpro.net/info/current/doc/manual/node417.html">http://www.molpro.net/info/current/doc/manual/node417.html</a><br>
  it is said <strong><tt><br><br>*********<br>   SAPT_FROZENA</tt></strong><dl><dd>Number of frozen electrons in the
      response calculations for monomer A (default 0)
</dd><dt>***************<br></dt></dl>  If I understood correctly, it is the number of electrons (number of orbital *2 in closed-shell case). But when I do a SAPT for neon dimer, setting SAPT_FROZENA=6<br><br>  the output file said <br>
<br>   <br> Occupied space for monomer A:    5<br> Virtual  space for monomer A:   87<br> <br> Occupied space for monomer B:    5<br> Virtual  space for monomer B:   87<br> ? Error<br> ? frozena>noa!<br> ? The problem occurs in sapt_interface<br>
<br>  does it imply the FROZENA actually is the number of orbital to be frozen?<br><br> (ii) I set FROZENA=0,1,2; FROZENB=1, the  E1tot+E2tot  is invariant respect with the frozen core set up. Why the frozen core does not change the interaction energy?<br>
<br>Here is an input file, if you need any other information, please let me known<br><br>*************************<br>memory,80,m<br>symmetry,nosym<br>geomtyp=xyz<br>geometry={<br>2<br>title<br>Ne 0. 0. 0.<br>Ne 0. 0. 3.0}<br>
 <br>basis=avtz<br> <br>!wf records<br>ca=2101.2<br>cb=2102.2<br> <br> <br>!monomer A<br>dummy,2<br>{hf; save,$ca}<br>sapt;monomerA<br> <br>!monomer B<br>dummy,1<br>{hf; start,atdens; save,$cb}<br>sapt;monomerB<br> <br>!interaction contributions<br>
{sapt,SAPT_FROZENA=6;intermol,ca=$ca,cb=$cb}<br> <br>*************************<br><br><br>Thank you very much in advance<br>Best regards<br>Cong Wang<br>Ph. D. Student<br><br>Department of Chemistry
<br>Laboratory for Instruction in Swedish
<br>University of Helsinki
<br>A.I. Virtanens plats 1
<br>P.O. Box 55
<br>FI-00014 University of Helsinki
<br>FINLAND

<br> <br>