19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

12 EFFECTIVE CORE POTENTIALS 81<br />

12.3 Example for explicit ECP input<br />

***,CU<br />

! SCF d10s1 -> d9s2 excitation energy of the Cu atom<br />

! using the relativistic Ne-core pseudopotential<br />

! and basis of the Stuttgart/Koeln group.<br />

gprint,basis,orbitals<br />

geometry={cu}<br />

basis<br />

ECP,1,10,3; ! ECP input<br />

1; ! NO LOCAL POTENTIAL<br />

2,1.,0.;<br />

2; ! S POTENTIAL<br />

2,30.22,355.770158;2,13.19,70.865357;<br />

2; ! P POTENTIAL<br />

2,33.13,233.891976;2,13.22,53.947299;<br />

2; ! D POTENTIAL<br />

2,38.42,-31.272165;2,13.26,-2.741104;<br />

! (8s7p6d)/[6s5p3d] BASIS SET<br />

s,1,27.69632,13.50535,8.815355,2.380805,.952616,.112662,.040486,.01;<br />

c,1.3,.231132,-.656811,-.545875;<br />

p,1,93.504327,16.285464,5.994236,2.536875,.897934,.131729,.030878;<br />

c,1.2,.022829,-1.009513;C,3.4,.24645,.792024;<br />

d,1,41.225006,12.34325,4.20192,1.379825,.383453,.1;<br />

c,1.4,.044694,.212106,.453423,.533465;<br />

end<br />

rhf;<br />

e1=energy<br />

{rhf;occ,4,1,1,1,1,1,1;closed,4,1,1,1,1,1;wf,19,7,1;}<br />

e2=energy<br />

de=(e2-e1)*toev<br />

! Delta E = -0.075 eV<br />

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

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

Saved successfully!

Ooh no, something went wrong!