<div dir="ltr">Dear Professor Peterson,<div><br></div><div>Thank you very much for your reply. Actually, I have acquired the spin-orbit basis for the Rb atom from Professor <span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:13px">Michael</span><span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:13px"> </span>Dolg. He gave me this link <a href="http://www.tc.uni-koeln.de/PP/clickpse.en.html" target="_blank" style="font-size:13px;font-family:arial,sans-serif">http://www.tc.uni-koeln.de/PP/<u></u>clickpse.en.html</a> to download the spin-orbit basis from the periodic table. This spin-orbit basis is the "ECP28MDF" for the Rb atom. As for the Yb atom he also gave this SO basis which is the ECP60MDF. The MDF means that the basis is fully relativistic, isn't it?! Do you think that these basis aren't spin-orbit ECP basis?</div>
<div><br></div><div>Thank you so much Professor for always answering my questions.</div><div><br></div><div>Regards,</div><div>Samir</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 27, 2014 at 6:49 PM, Peterson, Kirk <span dir="ltr"><<a href="mailto:kipeters@wsu.edu" target="_blank">kipeters@wsu.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Samir,<br>
<br>
I think the problem might be that the ECP you've chosen for Rb doesn't have any spin-orbit terms, only scalar relativistic ones. I would recommend trying the ecp28mdf ecp and basis set from the molpro library. Do you realize that the ECP you've chosen for Yb is non-relativistic? I would recommend ecp28mwb with the corresponding ecp28mwb_ano basis set.<br>
<br>
best regards,<br>
<br>
-Kirk<br>
<div><div class="h5"><br>
On Aug 26, 2014, at 10:28 AM, Samir Tohme <<a href="mailto:samstohme@gmail.com">samstohme@gmail.com</a>> wrote:<br>
<br>
> Dear MOLPRO Developers and Users,<br>
><br>
> I am applying a spin-orbit calculation for the Alkali-Ytterbium molecules particularly the RbYb compound. However, I am facing a problem in the CI program which is the following:<br>
><br>
> 1PROGRAM * CI (Multireference internally contracted CI) Authors: H.-J. Werner, P.J. Knowles, 1987<br>
><br>
><br>
> Spin-Orbit Integrals ECPLSX not found<br>
><br>
> ERROR EXIT<br>
> CURRENT STACK: CIINP CIPRO MAIN<br>
><br>
><br>
> **********************************************************************************************************************************<br>
> DATASETS * FILE NREC LENGTH (MB) RECORD NAMES<br>
> 1 22 8.64 600 500 700 960 900 950 970 1000 1100 1400<br>
> BASINP VAR GEOM ABASIS SYMINP ZMAT AOBASIS BASIS S T<br>
> 1410 1420 1200 1210 1080 1600 129 1650 1300 1700<br>
> V ECP H0 H01 AOSYM SMH P2S MOLCAS ERIS OPER<br>
> 1380 1700(1)<br>
> JKOP OPER<br>
><br>
> 2 14 19.87 700 1000 2100 7100 2101 1 5100 5200 5300 5400<br>
> GEOM BASIS RHF MCSCF RHF MCSCF MRCI MRCI MRCI MRCI<br>
> 6100 6200 6300 6400<br>
> MRCI MRCI MRCI MRCI<br>
><br>
> PROGRAMS * TOTAL CI CI CI CI CI CI CI CI CI CI<br>
> CPU TIMES * 277747.90 0.00 18.02 18.35 18.39 15.45 37.88 50.11 49.96 55.76 9302.55<br>
> REAL TIME * 278161.59 SEC<br>
> DISK USED * 673.69 MB<br>
> **********************************************************************************************************************************<br>
><br>
> It seems that the error is in the basis of the Rb and Yb atoms. I have inserted the spin-orbit basis for the two latter atoms, but the error persists! The following are the basis chosen:<br>
><br>
> basis={<br>
> ! LANL2DZ ECP EMSL Basis Set Exchange Library 7/13/13 5:02 AM<br>
> ! Elements References<br>
> ! -------- ----------<br>
> ! H - Ne: T. H. Dunning Jr. and P. J. Hay, in Methods of Electronic Structure<br>
> ! Theory, Vol. 2, H. F. Schaefer III, ed., PLENUM PRESS (1977)<br>
> ! Na - Hg: P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 270 (1985).<br>
> ! P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 284 (1985).<br>
> ! P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 299 (1985).<br>
> !<br>
> !<br>
> ! RUBIDIUM (5s,6p) -> [3s,3p]<br>
> ! RUBIDIUM (5s,6p) -> [3s,3p]<br>
> s, RB , 1.4570000, 0.9400000, 0.2706000, 0.0366000, 0.0155000<br>
> c, 1.3, -1.0617406, 1.1722125, 0.7437008<br>
> c, 1.4, 0.5052174, -0.5718404, -0.6835308, 1.1136477<br>
> c, 5.5, 1<br>
> p, RB , 3.3000000, 0.5985000, 0.2067000, 0.1947000, 0.0318000, 0.0124000<br>
> c, 1.3, -0.0729417, 0.6321782, 0.4707426<br>
> c, 4.5, -0.1250552, 1.0438206<br>
> c, 6.6, 1<br>
><br>
> ! Elements References<br>
> ! -------- ----------<br>
> ! Na - Hg: P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 270 (1985).<br>
> ! P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 284 (1985).<br>
> ! P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 299 (1985).<br>
> !<br>
> ! Effective core Potentials<br>
> ! -------------------------<br>
> ECP, rb, 28, 3 ;<br>
> 5; ! f-ul potential<br>
> 0,879.5983664,-0.0364362;<br>
> 1,142.0188695,-20.8033052;<br>
> 2,40.9346882,-115.4661910;<br>
> 2,9.9084802,-42.5061293;<br>
> 2,3.1556041,-5.5247889;<br>
> 4; ! s-ul potential<br>
> 0,70.9807215,2.9747910;<br>
> 1,38.6949807,36.1170435;<br>
> 2,15.9082035,152.2741068;<br>
> 2,3.7152659,53.7151906;<br>
> 5; ! p-ul potential<br>
> 0,93.3058684,4.9551026;<br>
> 1,47.1221370,24.2726236;<br>
> 2,48.8050640,215.2341411;<br>
> 2,14.3784949,151.6687666;<br>
> 2,3.0547807,33.8750280;<br>
> 4; ! d-ul potential<br>
> 0,33.0477430,3.0465224;<br>
> 1,9.9601880,29.4243624;<br>
> 2,3.7403577,33.4431808;<br>
> 2,0.8124208,11.1029460;<br>
><br>
> ECP,Yb,ECP60MHF<br>
><br>
> !Effective Core Potentials for spin orbit calculation<br>
> ECP,Yb,60,4,3;<br>
> 1;<br>
> 2, 1.00000000000000, 0.00000000000000;<br>
> 4;<br>
> 2, 9.29926500000000, 333.15106900000001;<br>
> 2, 4.50235700000000, -27.46097500000000;<br>
> 2, 2.24999700000000, -1.84075300000000;<br>
> 2, 1.12000000000000, 3.00925900000000;<br>
> 8;<br>
> 2, 7.44315400000000, 100.72952733333334;<br>
> 2, 7.44315400000000, 258.62780333333336;<br>
> 2, 3.74225300000000, -29.76952366666667;<br>
> 2, 3.74225300000000, -73.25967866666666;<br>
> 2, 2.37460900000000, 8.01055533333333;<br>
> 2, 2.37460900000000, 21.24672800000000;<br>
> 2, 1.18000000000000, 0.07348133333333;<br>
> 2, 1.18000000000000, -0.22957200000000;<br>
> 8;<br>
> 2, 4.75037000000000, 99.11122440000000;<br>
> 2, 4.75037000000000, 181.12201560000000;<br>
> 2, 2.37515300000000, -5.40595800000000;<br>
> 2, 2.37515300000000, -8.86254480000000;<br>
> 2, 1.18513300000000, -1.17406884000000;<br>
> 2, 1.18513300000000, -1.91775030000000;<br>
> 2, 0.58975800000000, -0.00228124000000;<br>
> 2, 0.58975800000000, -0.00457044000000;<br>
> 2;<br>
> 2, 1.38640000000000, 3.00707571428571;<br>
> 2, 1.38640000000000, 4.00943428571429;<br>
> 8;<br>
> 2, 7.44315400000000, -201.45905466666667;<br>
> 2, 7.44315400000000, 258.62780333333336;<br>
> 2, 3.74225300000000, 59.53904733333334;<br>
> 2, 3.74225300000000, -73.25967866666666;<br>
> 2, 2.37460900000000, -16.02111066666667;<br>
> 2, 2.37460900000000, 21.24672800000000;<br>
> 2, 1.18000000000000, -0.14696266666667;<br>
> 2, 1.18000000000000, -0.22957200000000;<br>
> 8;<br>
> 2, 4.75037000000000, -99.11122440000000;<br>
> 2, 4.75037000000000, 120.74801040000000;<br>
> 2, 2.37515300000000, 5.40595800000000;<br>
> 2, 2.37515300000000, -5.90836320000000;<br>
> 2, 1.18513300000000, 1.17406884000000;<br>
> 2, 1.18513300000000, -1.27850020000000;<br>
> 2, 0.58975800000000, 0.00228124000000;<br>
> 2, 0.58975800000000, -0.00304696000000;<br>
> 2;<br>
> 2, 1.38640000000000, -2.00471714285714;<br>
> 2, 1.38640000000000, 2.00471714285714;<br>
><br>
> s,Yb,12.222735,8.147737,5.078666,0.734744,0.319998,0.070091,0.027363;<br>
> p,Yb,7.086291,4.723029,1.133605,0.519137,0.227710,0.080626,0.02;<br>
> d,Yb,3.237970,1.377363,0.492827,0.168078,0.058334,0.02;<br>
> }<br>
><br>
> May anyone assist me how to solve this problem?<br>
><br>
> Thank you very much in advance.<br>
><br>
> Best regards,<br>
> Samir<br>
><br>
</div></div>> _______________________________________________<br>
> Molpro-user mailing list<br>
> <a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>
> <a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br>
<br>
</blockquote></div><br></div>