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.

34 PROPERTIES AND EXPECTATION VALUES 262<br />

34.4.2 Saving the NLMO orbitals (SAVE)<br />

SAVE,record.file;<br />

The NLMO orbitals are saved in the specified record, together with the NPA charges.<br />

34.5 Finite field calculations<br />

Dipole moments, quadrupole moments etc. and the corresponding polarizabilities can be obtained<br />

as energy derivatives by the finite difference approximation. This is most easily done<br />

with the DIP, QUAD, or FIELD commands. An error will result if the added perturbation is not<br />

totally symmetric (symmetry 1). Note that the orbitals must be recomputed before performing<br />

a correlation calculation.<br />

34.5.1 Dipole fields (DIP)<br />

DIP,xfield,yfield,zfield;<br />

DIP+,xfield,yfield,zfield;<br />

Add a finite dipole field to the one electron Hamiltonian and the core energy. The field strength<br />

is given by xfield,yfield,zfield. DIP+ adds to any existing field, otherwise any previous field is<br />

removed.<br />

34.5.2 Quadrupole fields (QUAD)<br />

QUAD,xxfield,yyfield,zzfield,xyfield,xzfield,yzfield;<br />

QUAD+,xxfield,yyfield,zzfield,xyfield,xzfield,yzfield;<br />

Exactly as the DIP command, but adds a quadrupole field.<br />

34.5.3 General fields (FIELD)<br />

FIELD,oper1,fac1, oper2,fac2, . . . ;<br />

FIELD+,oper1,fac1, oper2,fac2, . . . ;<br />

Adds one-electron operators oper1, oper2, . . . with the corresponding factors fac1, fac2, . . . to<br />

the one-electron hamiltonian. The available operators are given in section 6.13. An error will<br />

result if the added perturbation is not totally symmetric (symmetry 1).<br />

FIELD+ adds to any existing field, otherwise any previous field is removed.<br />

Note that FIELD does currently not modify core polarization potentials (CPP). If CPPs are<br />

present, only DIP and QUAD should be used.<br />

34.5.4 Examples<br />

The first examples shows various possibilities to add perturbations to the one-electron hamiltonian.

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

Saved successfully!

Ooh no, something went wrong!