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.

34 PROPERTIES AND EXPECTATION VALUES 263<br />

***,H2O finite fields<br />

memory,4,m<br />

R = 0.96488518 ANG<br />

THETA = 101.90140469<br />

geometry={H1<br />

O,H1,R;<br />

H2,O,R,H1,THETA}<br />

{hf;wf,10,1}<br />

!scf without field<br />

f=0.05<br />

dip,,,f<br />

hf<br />

field,dmz,f<br />

hf<br />

quad,,,f<br />

hf<br />

!add dipole (z) field to h0<br />

!do scf with modified h0<br />

!add dipole (z) field to H0<br />

!same result as previous example<br />

!do scf with modified h0<br />

!add quadrupole (qmzz) field to h0<br />

!do scf with modified h0<br />

field,qmzz,f !add quadrupole (qmzz) field to h0;<br />

!same result as previous example<br />

hf<br />

!do scf with modified h0<br />

field,zz,f,xx,-0.5*f,yy,-0.5*f<br />

!add general field; same result as quad above<br />

hf<br />

!do scf with modified h0<br />

field,zz,f<br />

field+,xx,-0.5*f<br />

field+,yy,-0.5*f<br />

hf<br />

field<br />

hf<br />

!same as before with separate field commands<br />

!do scf with modified h0<br />

!remove field<br />

!scf without field<br />

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

The second example shows how to compute dipole moments and polarizabilities using finite<br />

fields.

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

Saved successfully!

Ooh no, something went wrong!