<div dir="ltr">Dear colleagues,<div><br></div><div>I try to calculate potential curves for Ar+He with ECP for Ar. I used this code:</div><div><div><font face="monospace, monospace" size="1">***,ArHe</font></div><div><font face="monospace, monospace" size="1">   memory,100,M</font></div><div><font face="monospace, monospace" size="1">   gprint,civector,orbital=2;</font></div><div><font face="monospace, monospace" size="1">   geometry={angstrom</font></div><div><font face="monospace, monospace" size="1">   Ar</font></div><div><font face="monospace, monospace" size="1">   He,Ar,R(i)}</font></div><div><font face="monospace, monospace" size="1">   R=[20.,3.3,3.1,3.0,2.9,2.6,2.4,2.]</font></div><div><font face="monospace, monospace" size="1"><br></font></div><div><font face="monospace, monospace" size="1"> basis={</font></div><div><font face="monospace, monospace" size="1">ECP, ar, 10, 4 ;</font></div><div><font face="monospace, monospace" size="1">1; !  g-ul potential</font></div><div><font face="monospace, monospace" size="1">2,1.000000000,0.000000000;</font></div><div><font face="monospace, monospace" size="1">2; !  s-ul potential</font></div><div><font face="monospace, monospace" size="1">2,10.261721000,68.667788010;</font></div><div><font face="monospace, monospace" size="1">2,3.952725000,24.042766290;</font></div><div><font face="monospace, monospace" size="1">2; !  p-ul potential</font></div><div><font face="monospace, monospace" size="1">2,5.392714000,27.730763310;</font></div><div><font face="monospace, monospace" size="1">2,2.699967000,4.045459040;</font></div><div><font face="monospace, monospace" size="1">2; !  d-ul potential</font></div><div><font face="monospace, monospace" size="1">2,8.086235000,-8.137476960;</font></div><div><font face="monospace, monospace" size="1">2,4.016632000,-1.664528080;</font></div><div><font face="monospace, monospace" size="1">1; !  f-ul potential</font></div><div><font face="monospace, monospace" size="1">2,5.208459000,-3.400098450;</font></div><div><font face="monospace, monospace" size="1">spdfgh,1,aug-cc-pCV5Z;C;</font></div><div><font face="monospace, monospace" size="1">SP,1,EVEN,NPRIM=3,RATIO=2.5;</font></div><div><font face="monospace, monospace" size="1">spdfg,2,aug-cc-pV5Z;C;</font></div><div><font face="monospace, monospace" size="1"> }</font></div><div><font face="monospace, monospace" size="1">do i=1,#R<br></font></div><div><font face="monospace, monospace" size="1">{RHF</font></div><div><font face="monospace, monospace" size="1">wf,20,1,2</font></div><div><font face="monospace, monospace" size="1">}</font></div><div><font face="monospace, monospace" size="1"> {casscf</font></div><div><font face="monospace, monospace" size="1">closed,3,1,1,0</font></div><div><font face="monospace, monospace" size="1">occ,8,3,3,0;</font></div><div><font face="monospace, monospace" size="1">DYNW,2</font></div><div><font face="monospace, monospace" size="1"> wf,20,1,2;state,4</font></div><div><div><font face="monospace, monospace" size="1"> wf,20,2,2;state,3</font></div><div><font face="monospace, monospace" size="1"> wf,20,3,2;state,3</font></div><div><font face="monospace, monospace" size="1"> wf,20,4,2;state,2</font></div><div><font face="monospace, monospace" size="1"> wf,20,1,0;state,5</font></div><div><font face="monospace, monospace" size="1"> wf,20,2,0;state,3</font></div><div><font face="monospace, monospace" size="1"> wf,20,3,0;state,3</font></div><div><font face="monospace, monospace" size="1"> wf,20,4,0;state,2</font></div><div><font face="monospace, monospace" size="1">}</font></div><div><font face="monospace, monospace" size="1">  {ci;closed,3,1,1,0;occ,8,3,3,0;orbital,ignore_error;wf,sym=1,SPIN=2;state,4;save,5101.2}</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(1)=energd(1)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(2)=energd(2)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(3)=energd(3)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(4)=energd(4)</font></div><div><font face="monospace, monospace" size="1">  {ci;closed,3,1,1,0;occ,8,3,3,0;wf,sym=2,SPIN=2;state,3;save,5102.2}</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(5)=energd(1)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(6)=energd(2)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(7)=energd(3)</font></div><div><font face="monospace, monospace" size="1">  {ci;closed,3,1,1,0;occ,8,3,3,0;wf,sym=3,SPIN=2;state,3;save,5103.2}</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(8)=energd(1)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(9)=energd(2)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(10)=energd(3)</font></div><div><font face="monospace, monospace" size="1">  {ci; closed,3,1,1,0;occ,8,3,3,0;wf,sym=4,SPIN=2;state,2;save,5104.2}</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(11)=energd(1)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(12)=energd(2)</font></div><div><font face="monospace, monospace" size="1">  {ci; closed,3,1,1,0;occ,8,3,3,0;wf,sym=1,SPIN=0;State,5;save,5105.2}</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(13)=energd(1)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(14)=energd(2)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(15)=energd(3)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(16)=energd(4)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(17)=energd(5)</font></div><div><font face="monospace, monospace" size="1">  {ci; closed,3,1,1,0;occ,8,3,3,0;wf,sym=2,SPIN=0;state,3;save,5106.2}</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(18)=energd(1)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(19)=energd(2)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(20)=energd(3)</font></div><div><font face="monospace, monospace" size="1">  {ci; closed,3,1,1,0;occ,8,3,3,0;wf,sym=3,SPIN=0;state,3;save,5107.2}</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(21)=energd(1)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(22)=energd(2)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(23)=energd(3)</font></div><div><font face="monospace, monospace" size="1">  {ci; closed,3,1,1,0;occ,8,3,3,0;wf,sym=4,SPIN=0;state,2;save,5108.2}</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(24)=energd(1)</font></div><div><font face="monospace, monospace" size="1"> hlsdiag(25)=energd(2)</font></div></div><div><div><font face="monospace, monospace" size="1"><br></font></div><div><font face="monospace, monospace" size="1">  {ci;hlsmat,ls,5101.2,5102.2,5103.2,5104.2,5105.2,5106.2,5107.2,5108.2}</font></div><div><font face="monospace, monospace" size="1"><br></font></div><div><font face="monospace, monospace" size="1">enddo</font></div></div><div><br></div><div>On the RHF calculation Molpro wrote this text:<br><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">...</font></div><div><font face="monospace, monospace">ATOMIC COORDINATES</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> NR  ATOM    CHARGE       X              Y              Z</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">   1  AR      8.00    0.000000000    0.000000000   -3.441963392</font></div><div><font face="monospace, monospace">   2  HE      2.00    0.000000000    0.000000000   34.352559236</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> NUCLEAR CHARGE:                   10</font></div><div><font face="monospace, monospace"> NUMBER OF PRIMITIVE AOS:         442</font></div><div><font face="monospace, monospace"> NUMBER OF SYMMETRY AOS:          327</font></div><div><font face="monospace, monospace"> NUMBER OF CONTRACTIONS:          309   ( 116A1  +  75B1  +  75B2  +  43A2  )</font></div><div><font face="monospace, monospace"> NUMBER OF CORE ORBITALS:           0   (   0A1  +   0B1  +   0B2  +   0A2  )</font></div><div><font face="monospace, monospace"> NUMBER OF VALENCE ORBITALS:        5   (   3A1  +   1B1  +   1B2  +   0A2  )</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> NUCLEAR REPULSION ENERGY    0.42334177</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> One-electron integrals computed with SEWARD</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> 2nd-order Douglas-Kroll-Hess method activated. Optimal DKH parametrization is used.</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> GLOBAL ERROR fehler on processor   0</font></div></div><div><br></div><div>How to check this problem?</div><div><br></div><div>Best Regards,<br></div><div>Andrey<br></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div style="font-size:small"><br></div><div style="font-size:small">Postgraduate Student Andrey Pershin<br>Samara University, Samara, Russia<br></div></div></div>
</div></div><div id="DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br> <table style="border-top:1px solid #d3d4de">
        <tr>
      <td style="width:55px;padding-top:18px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
                <td style="width:470px;padding-top:17px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Без вирусов. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank" style="color:#4453ea">www.avast.ru</a>                 </td>
        </tr>
</table>
<a href="#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"></a></div>