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.

18 ORBITAL LOCALIZATION 124<br />

geometry={c;n,c,r}<br />

r=1.1 angstrom<br />

df=[b,lyp]<br />

rhf;method(1)=program<br />

dft;edf(1)=dftfun<br />

uhf;method(2)=program<br />

dft;edf(2)=dftfun<br />

uks;method(3)=program,edf(3)=dftfun<br />

dft;method(4)=program,edf(4)=dftfun<br />

table,dftname,dftfuns<br />

table,method,edf<br />

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

18 ORBITAL LOCALIZATION<br />

Localized orbitals are calculated according to the Boys, Pipek-Mezey or NLMO criteria. Localization<br />

takes place within each symmetry species separately. If complete localization is desired,<br />

no symmetry should be used. All subcommands can be abbreviated by three characters.<br />

The localization program is invoked by the LOCALI command<br />

LOCALI [,method],[LOCMETHOD=locmethod],[REFORB=record]<br />

The keyword method can be either BOYS, PIPEK or NATURAL. By default, the valence orbitals<br />

from the last energy calculation are localized using the Boys criterion. Only orbital subsets<br />

which leave the energy invariant are transformed. These defaults can be modified using the<br />

optional commands described in the following sections.<br />

The option LOCMETHOD only applies to Pipek-Mezey localization. The value locmethod can<br />

take the following values:<br />

locmethod=1:<br />

Standard iterative localization method<br />

locmethod=2: Use second-order localization method. Redundant rotations<br />

will be eliminated.<br />

locmethod=3:<br />

First do a first iterations using the standard method, then invoke<br />

locmethod=2. This option is recommended in cases with<br />

redundant rotations, e.g., benzene.<br />

The option REFORB has the same effect as the directive REFORB described further below.<br />

18.1 Defining the input orbitals (ORBITAL)<br />

ORBITAL,record.file,specifications<br />

The orbitals to be localized are read from dump record record.file. A state specific orbital set can<br />

be selected using specifications, as explained in section 4.11. Default are the orbitals calculated<br />

last.<br />

18.2 Saving the localized orbitals (SAVE)<br />

SAVE,record.file<br />

This specifies the dump record where the localized orbitals are stored. If the dump record<br />

already exists, the localized orbitals are added to it. Default is the input record (cf. ORBITAL).

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

Saved successfully!

Ooh no, something went wrong!