[molpro-user] puzzling XeF2 result (corrected)

Hoffman, Gerald GHOFFMAN at edinboro.edu
Tue Jul 25 14:44:54 CEST 2017


I made an embarrassing mistake in my previous post, which has been corrected below.  However, I'm still getting the wrong energy:

Using the input file below, I attempted an SCF calculation on the XeF2 molecule using the Stuttgart ECP and basis set for Xe and the Aug-cc-pVQZ basis set for F.  I obtained a total energy of -58.43 Hartrees.  This is clearly wrong.  Each of the fluorines by itself ought to have an energy of -99. Hartrees, and in previous versions of Molpro, the energy of this molecule is around -214. Hartrees.  The calculation of XeF gives a reasonable value (-115. Hartrees).

Help would be appreciated.

Input file:

***,XeF2
memory,50,M

angstrom

geometry={
Xe    0.0000    0.0000    0.0000
F1    0.0000    0.0000   -1.99
F2    0.0000    0.0000    1.99
}

basis={
!
! Xenon ecp
!
ecp,Xe,46,5,3;
1;                                                                              ! h
2,1.00,0.00
2;                                                                              ! s-h
2,3.940263,122.76382934;2,2.277264,8.30885115;
2;                                                                              ! p-h
2,3.028373,68.82300437;2,1.394319,3.64674233;
2;                                                                              ! d-h
2,2.122605,23.65207854;2,0.798669,3.25844113;
2;                                                                              ! f-h
2,6.164360,-47.70319876;2,1.542374,-6.54113991;
1;                                                                              ! g-h
2,1.847892,-7.10585060;
2;                                                                              ! so-ecp
2,3.028373,-1.74609113*(2/3);2,1.394319,2.19487257*(2/3);
2;                                                                              ! so-ecp
2,2.122605,-1.04567591*(2/5);2,0.798669,0.27199291*(2/5);
1;                                                                              ! so-ecp
2,6.164360,0.20772207*(2/7);
!
! Xenon basis
!
spdfg,Xe,ecp46mwb_vqz
s,Xe,0.040702,0.019340;c,1.1,1.00;c,2.2,1.00;
p,Xe,0.026108,0.012447;c,1.1,1.00;c,2.2,1.00;
d,Xe,0.089444;c,1.1,1.00;
f,Xe,0.132301;c,1.1,1.00;
!
! Fluorine-basis
!
spdfg,F,avqz
s,2,0.034376;c,1.1,1.00;
p,2,0.026272;c,1.1,1.00;
d,2,0.0828;c,1.1,1.00;
}

cpp,init,1;
Xe,1,0.8345,,,0.97391;

{rhf;wf,26,1,0;}



Gerry Hoffman


Gerald J. Hoffman, Ph. D.

Associate Professor of Chemistry

Edinboro University of Pennsylvania

230 Scotland Road

Edinboro, PA 16444



814-732-2813



ghoffman at edinboro.edu<mailto:ghoffman at edinboro.edu>


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20170725/d83f1c84/attachment-0001.html>


More information about the Molpro-user mailing list