<html><body bgcolor="#FFFFFF"><div>Aleem,</div><div><br></div><div>When you use a ECP, core electrons are being replaced (that's why the total energies change so much, as Gerald has pointed out). Thus you have to drop these from your wf, occ, closed, con, etc directives. In your case you have accidentally calculated props of the doubly negative ions. </div><div><br></div><div>Regards.</div><div><br></div><div>Kirk<br><br>Sent from my iPhone</div><div><br>On Jan 3, 2011, at 11:36 PM, syed aleem <<a href="mailto:aleem3176@gmail.com">aleem3176@gmail.com</a>> wrote:<br><br></div><div></div><blockquote type="cite"><div>Dear Kirk,<br>     <br>Here am providing my complete input for both calculations.<br><br>Input with ECP<br><br>***,ionization potential<br>memory,100,m<br>basis={h=aug-cc-pvtz,<br>       c=aug-cc-pvtz,<br>       f=aug-cc-pvqz<br>
       ECP,F,2,3,0;<br>       1; 2,1.000000,0.000000;<br>       1; 2,16.762800,96.241900;<br>       1; 2,19.857600,-4.581900;<br>       1; 2,0.031600,-0.001500;}<br>geomtyp=xyz<br>geometry={5<br>equilibrium geometry<br>c     0.000000    0.000000   -0.636134<br>
h     0.000000    1.029894   -0.984650<br>h    -0.891914   -0.514947   -0.984650<br>h     0.891914   -0.514947   -0.984650<br>f     0.000000    0.000000    0.752306<br>end}<br> rhf;occ,7,2;wf,18,1,0;orbital,2130.2;orbprint,5;<br>
 multi;occ,9,3;closed,4,0;start,2130.2;wf,17,1,1;state,2;                        <br> mrci;option,maxit=100;occ,9,3;closed,4,0;wf,17,1,1;state,2;select;con,2,2,1,0,0,2,2,0;con,2,1,2,0,0,2,2,0;<br> exit;<br><br>Input without ECP<br>
***,ionization potential<br>memory,100,m<br>basis={h=aug-cc-pvtz,<br>       c=aug-cc-pvtz,<br>       f=aug-cc-pvqz;}<br>      !ECP,F,2,3,0;<br>      !1; 2,1.000000,0.000000;<br>      !1; 2,16.762800,96.241900;<br>      !1; 2,19.857600,-4.581900;<br>
      !1; 2,0.031600,-0.001500;}<br>geomtyp=xyz<br>geometry={5<br>equilibrium geometry<br>c     0.000000    0.000000   -0.636134<br>h     0.000000    1.029894   -0.984650<br>h    -0.891914   -0.514947   -0.984650<br>h     0.891914   -0.514947   -0.984650<br>
f     0.000000    0.000000    0.752306<br>end}<br>rhf;occ,7,2;wf,18,1,0;orbital,2130.2;orbprint,5;<br>multi;occ,9,3;closed,4,0;start,2130.2;wf,17,1,1;state,2;   <br>mrci;option,maxit=100;occ,9,3;closed,4,0;wf,17,1,1;state,2;select;con,2,2,1,0,0,2,2,0;con,2,1,2,0,0,2,2,0;<br>
exit;<br><br>Thanks for your reply,<br><br>Regards,<br>Aleem.<br><br><br><div class="gmail_quote">On 4 January 2011 12:26, Kirk Peterson <span dir="ltr"><<a href="mailto:kipeters@wsu.edu"><a href="mailto:kipeters@wsu.edu">kipeters@wsu.edu</a></a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div bgcolor="#FFFFFF"><div>Aleem,</div><div><br></div><div>It's difficult to see where the problem lies without the full inputs. Your ECP input looks ok. Note that your F basis should probably be used uncontracted since it's contraction coeffs are optimized for all-electron calculations. </div>
<div><br></div><div>-Kirk<br><br>Sent from my iPhone</div><div><div></div><div class="h5"><div><br>On Jan 3, 2011, at 9:54 PM, syed aleem <<a href="mailto:aleem3176@gmail.com" target="_blank"><a href="mailto:aleem3176@gmail.com">aleem3176@gmail.com</a></a>> wrote:<br>
<br></div><div></div><blockquote type="cite"><div>Dear Users,<br><br>I am trying to use relativistic effective core potential for core electrons of Fluorine atom in ch3f molecule<br><br>My input according to manual is<br>
<br>***,ionization potential<br>memory,100,m<br>basis={h=aug-cc-pvtz,<br>
       c=aug-cc-pvtz,<br>       f=aug-cc-pvqz<br>       ECP,F,2,3,0;<br>       1; 2,1.000000,0.000000;<br>       1; 2,16.762800,96.241900;<br>       1; 2,19.857600,-4.581900;<br>       1; 2,0.031600,-0.001500;}<br>geomtyp=xyz<br>

geometry={5<br>geometry parameters,wave function for multi and mrci calculation were followed<br><br>But the problem is when I submitted this job. state energies were drastically changed. I am providing both outputs (one with ecp and one without)<br>

<br><u><b>With ECP</b></u><br><br>     MRCI            MRCI            MULTI           MULTI           RHF-SCF<br>    -64.08157998    -64.10092216    -63.67925159    -63.69710084    -63.49847009<br><br><u><b>With out ECP</b></u><br>

<br>               MRCI            MRCI            MULTI           MULTI           RHF-SCF<br>   -138.93548516   -139.06809598   -138.57347287   -138.70164845   -139.10421155<br><br>Could you please check whether there was any mistake in my input. I would be grateful if you explain how to deal with  ECP (effective core potentials in molpro calculations)<br>

<br>Thanking you in advance for your help<br><br>Regards,<br><br>Aleem,<br>School of Chemistry,<br>University of Hyderabad,<br>INDIA.<br><br>
</div></blockquote></div></div><blockquote type="cite"><div><span>_______________________________________________</span><br><span>Molpro-user mailing list</span><br><span><a href="mailto:Molpro-user@molpro.net" target="_blank"><a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a></a></span><br>
<span><a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank"><a href="http://www.molpro.net/mailman/listinfo/molpro-user">http://www.molpro.net/mailman/listinfo/molpro-user</a></a></span><br></div></blockquote></div></blockquote></div><br>
</div></blockquote></body></html>