<div dir="ltr"><div>Hans,</div><br>There's failing output<div><br></div><div>Andrey</div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-01-26 12:30 GMT+04:00 Hans-Joachim Werner <span dir="ltr"><<a href="mailto:werner@theochem.uni-stuttgart.de" target="_blank">werner@theochem.uni-stuttgart.de</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Andrey,<br>
I don’t see that problem when running your input with the current distribution (2015.1). Which program version did you use? There must be something wrong, since certainly DK should not be used with ECPs. Please send the failing output.<br>
Best regards<br>
Joachim<br>
<br>
---<br>
Prof. Dr. Hans-Joachim Werner<br>
Institut für Theoretische Chemie<br>
Universität Stuttgart<br>
Pfaffenwaldring 55<br>
70569 Stuttgart, Germany<br>
e-mail: <a href="mailto:werner@theochem.uni-stuttgart.de">werner@theochem.uni-stuttgart.<wbr>de</a><br>
<div><div class="h5"><br>
<br>
<br>
<br>
> Am 25.01.2018 um 15:03 schrieb Andrey Pershin <<a href="mailto:anchizh93@gmail.com">anchizh93@gmail.com</a>>:<br>
><br>
> Dear colleagues,<br>
><br>
> I try to calculate potential curves for Ar+He with ECP for Ar. I used this code:<br>
> ***,ArHe<br>
> memory,100,M<br>
> gprint,civector,orbital=2;<br>
> geometry={angstrom<br>
> Ar<br>
> He,Ar,R(i)}<br>
> R=[20.,3.3,3.1,3.0,2.9,2.6,2.<wbr>4,2.]<br>
><br>
> basis={<br>
> ECP, ar, 10, 4 ;<br>
> 1; ! g-ul potential<br>
> 2,1.000000000,0.000000000;<br>
> 2; ! s-ul potential<br>
> 2,10.261721000,68.667788010;<br>
> 2,3.952725000,24.042766290;<br>
> 2; ! p-ul potential<br>
> 2,5.392714000,27.730763310;<br>
> 2,2.699967000,4.045459040;<br>
> 2; ! d-ul potential<br>
> 2,8.086235000,-8.137476960;<br>
> 2,4.016632000,-1.664528080;<br>
> 1; ! f-ul potential<br>
> 2,5.208459000,-3.400098450;<br>
> spdfgh,1,aug-cc-pCV5Z;C;<br>
> SP,1,EVEN,NPRIM=3,RATIO=2.5;<br>
> spdfg,2,aug-cc-pV5Z;C;<br>
> }<br>
> do i=1,#R<br>
> {RHF<br>
> wf,20,1,2<br>
> }<br>
> {casscf<br>
> closed,3,1,1,0<br>
> occ,8,3,3,0;<br>
> DYNW,2<br>
> wf,20,1,2;state,4<br>
> wf,20,2,2;state,3<br>
> wf,20,3,2;state,3<br>
> wf,20,4,2;state,2<br>
> wf,20,1,0;state,5<br>
> wf,20,2,0;state,3<br>
> wf,20,3,0;state,3<br>
> wf,20,4,0;state,2<br>
> }<br>
> {ci;closed,3,1,1,0;occ,8,3,3,<wbr>0;orbital,ignore_error;wf,sym=<wbr>1,SPIN=2;state,4;save,5101.2}<br>
> hlsdiag(1)=energd(1)<br>
> hlsdiag(2)=energd(2)<br>
> hlsdiag(3)=energd(3)<br>
> hlsdiag(4)=energd(4)<br>
> {ci;closed,3,1,1,0;occ,8,3,3,<wbr>0;wf,sym=2,SPIN=2;state,3;<wbr>save,5102.2}<br>
> hlsdiag(5)=energd(1)<br>
> hlsdiag(6)=energd(2)<br>
> hlsdiag(7)=energd(3)<br>
> {ci;closed,3,1,1,0;occ,8,3,3,<wbr>0;wf,sym=3,SPIN=2;state,3;<wbr>save,5103.2}<br>
> hlsdiag(8)=energd(1)<br>
> hlsdiag(9)=energd(2)<br>
> hlsdiag(10)=energd(3)<br>
> {ci; closed,3,1,1,0;occ,8,3,3,0;wf,<wbr>sym=4,SPIN=2;state,2;save,<wbr>5104.2}<br>
> hlsdiag(11)=energd(1)<br>
> hlsdiag(12)=energd(2)<br>
> {ci; closed,3,1,1,0;occ,8,3,3,0;wf,<wbr>sym=1,SPIN=0;State,5;save,<wbr>5105.2}<br>
> hlsdiag(13)=energd(1)<br>
> hlsdiag(14)=energd(2)<br>
> hlsdiag(15)=energd(3)<br>
> hlsdiag(16)=energd(4)<br>
> hlsdiag(17)=energd(5)<br>
> {ci; closed,3,1,1,0;occ,8,3,3,0;wf,<wbr>sym=2,SPIN=0;state,3;save,<wbr>5106.2}<br>
> hlsdiag(18)=energd(1)<br>
> hlsdiag(19)=energd(2)<br>
> hlsdiag(20)=energd(3)<br>
> {ci; closed,3,1,1,0;occ,8,3,3,0;wf,<wbr>sym=3,SPIN=0;state,3;save,<wbr>5107.2}<br>
> hlsdiag(21)=energd(1)<br>
> hlsdiag(22)=energd(2)<br>
> hlsdiag(23)=energd(3)<br>
> {ci; closed,3,1,1,0;occ,8,3,3,0;wf,<wbr>sym=4,SPIN=0;state,2;save,<wbr>5108.2}<br>
> hlsdiag(24)=energd(1)<br>
> hlsdiag(25)=energd(2)<br>
><br>
> {ci;hlsmat,ls,5101.2,5102.2,<wbr>5103.2,5104.2,5105.2,5106.2,<wbr>5107.2,5108.2}<br>
><br>
> enddo<br>
><br>
> On the RHF calculation Molpro wrote this text:<br>
><br>
> ...<br>
> ATOMIC COORDINATES<br>
><br>
> NR ATOM CHARGE X Y Z<br>
><br>
> 1 AR 8.00 0.000000000 0.000000000 -3.441963392<br>
> 2 HE 2.00 0.000000000 0.000000000 34.352559236<br>
><br>
> NUCLEAR CHARGE: 10<br>
> NUMBER OF PRIMITIVE AOS: 442<br>
> NUMBER OF SYMMETRY AOS: 327<br>
> NUMBER OF CONTRACTIONS: 309 ( 116A1 + 75B1 + 75B2 + 43A2 )<br>
> NUMBER OF CORE ORBITALS: 0 ( 0A1 + 0B1 + 0B2 + 0A2 )<br>
> NUMBER OF VALENCE ORBITALS: 5 ( 3A1 + 1B1 + 1B2 + 0A2 )<br>
><br>
><br>
> NUCLEAR REPULSION ENERGY 0.42334177<br>
><br>
> One-electron integrals computed with SEWARD<br>
><br>
> 2nd-order Douglas-Kroll-Hess method activated. Optimal DKH parametrization is used.<br>
><br>
> GLOBAL ERROR fehler on processor 0<br>
><br>
> How to check this problem?<br>
><br>
> Best Regards,<br>
> Andrey<br>
><br>
> --<br>
><br>
> Postgraduate Student Andrey Pershin<br>
> Samara University, Samara, Russia<br>
><br>
</div></div><span class="">> Без вирусов. <a href="http://www.avast.ru" rel="noreferrer" target="_blank">www.avast.ru</a><br>
> ______________________________<wbr>_________________<br>
> Molpro-user mailing list<br>
> <a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>
</span>> <a href="http://www.molpro.net/mailman/listinfo/molpro-user" rel="noreferrer" target="_blank">http://www.molpro.net/mailman/<wbr>listinfo/molpro-user</a><br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div style="font-size:small">Андрей Першин<br></div><div style="font-size:small">аспирант Самарского университета</div><div style="font-size:small">тел. +79093430429</div><div style="font-size:small">e-mail: <a href="mailto:anchizh93@gmail.com" target="_blank">anchizh93@gmail.com</a></div><div style="font-size:small"><br></div><div style="font-size:small">Postgraduate Student Andrey Pershin</div><div style="font-size:small">Samara University, Samara, Russia</div><div style="font-size:small"><div>phone. +79093430429</div><div>e-mail: <a href="mailto:anchizh93@gmail.com" target="_blank">anchizh93@gmail.com</a></div></div></div></div></div></div>
</div>