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.

55 ORBITAL MERGING 378<br />

55.13 Examples<br />

55.13.1 H 2 F<br />

This example merges the orbitals of H 2 and F<br />

***,example for merge<br />

print,orbitals,basis<br />

rh2=1.4<br />

rhf=300.<br />

basis=vdz<br />

symmetry,x,y<br />

geometry={F}<br />

!use C2v symmetry<br />

text,F<br />

{rhf;wf,9,1,1;occ,3,1,1;orbital,2130.2} !rhf for f-atom<br />

text,H2<br />

symmetry,x,y<br />

geometry={<br />

H1,<br />

H2,H1,rh2}<br />

{hf;orbital,2100.2}<br />

{multi;occ,2;orbital,2101.2}<br />

text,FH2<br />

geometry={F;<br />

H1,F,rhf<br />

H2,H1,rh2,F,180}<br />

!use C2v symmetry<br />

!scf for h2<br />

!mcscf for h2<br />

!linear geometry for F+H2<br />

{merge<br />

orbital,2130.2<br />

!rhf orbitals for F-atom<br />

move,1.1,2.1,1.1 !move orbitals 1.1, 2.1<br />

move,3.1,0.4,4.1; !move all remaining, starting at 4.1<br />

orbital,2100.2<br />

!hf orbitals for H2<br />

move,1.1,0.4<br />

!move these to free positions<br />

save,2131.2}<br />

!save merged orbitals<br />

{rhf;occ,4,1,1;start,2131.2<br />

orbital,2132.2}<br />

!rhf for F+H2<br />

{merge<br />

orbital,2130.2<br />

!rhf orbitals for F-atom<br />

move,1.1,2.1,1.1 !move orbitals 1.1, 2.1<br />

move,3.1,3.1,4.1; !move orbital 3.1 to 4.1<br />

move,4.1,0.4,6.1 !move all remaining, starting at 6.1<br />

orbital,2101.2<br />

!mcscf orbitals for H2<br />

move,1.1,0.4<br />

!move these to free positions<br />

save,2141.2}<br />

!save merged orbitals<br />

{multi;occ,5,1,1;start,2141.2}<br />

!casscf for F+H2 using valence space<br />

http://www.molpro.net/info/current/examples/h2f_merge.com<br />

55.13.2 NO<br />

This example merges the SCF orbitals of N and O to get a full valence space for NO. In the<br />

simplest case the atomic calculations are performed in the individual separate basis sets, but<br />

using the same symmetry (C 2v ) as the molecular calculation.

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

Saved successfully!

Ooh no, something went wrong!