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.

17 THE DENSITY FUNCTIONAL PROGRAM 109<br />

TOLORB=value<br />

MATRIX=matrix<br />

Threshold for orbital screening (current default 1.d-15).<br />

Option to select integrator. matrix=0: use old (slow) integrator;<br />

matrix=1: Use new matrix-driven integrator (default).<br />

In addition, all options valid for HF (see section 16.1) can be given.<br />

17.2 Directives<br />

The following options may be used to control the operation of the DFT modules. In the Kohn-<br />

Sham case, these may come in any order before or after directives for the SCF program as<br />

described in Section 16.<br />

17.2.1 Density source (DENSITY, ODENSITY)<br />

DENSITY,orbc.filec,... ODENSITY,orbo.fileo,...<br />

For non-self-consistent DFT calculations, specifies the source of the density matrix. The total<br />

density is read from orbc.filec, with further options specifying density sets in the standard way<br />

as described in Section 4.11. ODENSITY can be used to specify the spin density. The defaults<br />

are the densities last written by an SCF or MCSCF program.<br />

17.2.2 Thresholds (DFTTHRESH)<br />

DFTTHRESH,key1=value1,key2=value2...<br />

Sets various truncation thresholds. key can be one of the following.<br />

TOTAL<br />

ORBITAL<br />

DENSITY<br />

FOCK<br />

Overall target accuracy (per atom) of density functional. Defaults<br />

to the value of the global threshold GRID or the value<br />

specified by option GRID. For proper use of this threshold,<br />

other thresholds should be left at their default value of zero.<br />

Orbital truncation threshold.<br />

Density truncation threshold.<br />

Fock matrix truncation threshold.<br />

17.2.3 Exact exchange computation (EXCHANGE)<br />

EXCHANGE,factor<br />

For Kohn-Sham calculations, compute exchange energy according to Hartree-Fock formalism<br />

and add the contribution scaled by factor to the fock matrix and the energy functional. Otherwise,<br />

the default is factor=0, i.e., the exchange is assumed to be contained in the functional, and<br />

only the Coulomb interaction is calculated explicitly.<br />

DFTFACTOR,fac1, fac2, ...<br />

Provide a factor for each functional specified. The functionals will be combined accordingly.<br />

By default, all factors are one.

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

Saved successfully!

Ooh no, something went wrong!