19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

28 SMILES 203<br />

Atom type index is used to establish the correspondence between the basis functions and the<br />

centers (atoms). All records having the same value of I will define the basis set for all the atoms<br />

of the corresponding type.<br />

Example:<br />

For a calculation on the CO 2 molecule with the following geometry data:<br />

Geometry={<br />

3<br />

C1, 0.000000000000E+00, 0.000000000000E+00, 0.000000000000E+00<br />

O2, 0.000000000000E+00, 0.000000000000E+00,-0.117963799946E+01<br />

O3, 0.000000000000E+00, 0.000000000000E+00, 0.117963799946E+01<br />

}<br />

Clementi and Roetti’s Single Zeta basis set could be supplied in an external file like:<br />

1 1 0 5.67263 12<br />

1 2 0 1.60833 12<br />

1 2 1 1.56788 12<br />

2 1 0 7.65781 12<br />

2 2 0 2.24588 12<br />

2 2 1 2.22662 12<br />

The first three records define the basis set for atoms of the first type (carbon in this example),<br />

and the following three, the basis set for atoms of the second type (oxygen).<br />

NG is mandatory even in case of diatomics, though gaussian expansions will not be used, the<br />

value of NG being irrelevant in this case.<br />

28.3 Example<br />

Example using internal basis set for H 2 O.<br />

memory,20,m<br />

if(.not.modul_slater) then<br />

exit<br />

end if<br />

geomtyp=zmat<br />

geometry={<br />

o;h1,o,r;h2,o,r,h1,theta<br />

}<br />

r=0.96 ang<br />

theta=102<br />

intyp=’SLATER’<br />

basis=VB1<br />

hf<br />

ccsd(t)<br />

multi<br />

mrci<br />

http://www.molpro.net/info/current/examples/slater.com

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!