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 127<br />

2.) Within each group of orbitals found for a given primary centre, those containing only<br />

one centre (lone pairs) are included first. The remaining ones are ordered so that the<br />

secondary atoms are in the order of centrelist. Orbitals with secondary centres which are<br />

not in centrelist come last.<br />

3.) If REVERT is given, the order in each localization group is reverted.<br />

4.) If GROUP is given, only the orbitals in the given group are reordered. igrp is 2 for closed<br />

shells and inactive orbitals, 1 for open-shells in single reference methods, and 3 for active<br />

orbitals in CASSCF calculations.<br />

5.) If THREIG is given, only orbitals with energies larger than the given value are reordered.<br />

eps must be negative. The remaining orbitals come last (first if REVERT is given).<br />

Note that core orbitals are neither localized nor reordered.<br />

18.7.3 Defining reference orbitals (REFORB)<br />

REFORB,record.file,specifications<br />

The localized orbitals are reordered such that the overlap with the reference orbitals read from<br />

record.file is maximized. This is useful for local correlation treatments for keeping the order of<br />

the localized constant for different geometries. A state specific orbital set can be selected using<br />

specifications, as explained in section 4.11.<br />

18.7.4 Selecting the fock matrix (FOCK)<br />

FOCK,record.file<br />

This specifies a record holding a Fock operator to be used for ordering the orbitals. Note that<br />

only SCF dump records hold fock operators. Default is the Fock operator from the energy<br />

calculation which produced the input orbitals.<br />

18.7.5 Selecting a density matrix (DENSITY)<br />

DENSITY,record.file,specifications<br />

This specifies a record holding a density matrix for construction of a fock operator used for<br />

ordering the orbitals. This can be used if no fock operator is available, and has only an effect for<br />

MCSCF localizations. By default, the (state averaged) MCSCF density is used. A state specific<br />

density matrix can be selected using specifications as described in section 4.11.<br />

18.8 Localization thresholds (THRESH)<br />

THRESH,thresh,eorder<br />

thresh is a threshold for localization (default 1.d-12). If eorder is nonzero (default 1.d-4), the<br />

orbitals whose energy difference is smaller then eorder are considered to be degenerate and<br />

reordered according to the position of their charge centres (see section 18.7).

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

Saved successfully!

Ooh no, something went wrong!