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.

16 THE SCF PROGRAM 99<br />

MAXIT=maxit Maximum number of iterations (default 60)<br />

SHIFTA|SHIFTC=shifta<br />

SHIFTB|SHIFTO=shiftb<br />

Level shift for closed-shell orbitals in RHF (default −0.3) and<br />

α-spin orbitals in UHF (default 0).<br />

Level shift for open-shell orbitals in RHF and β-spin orbitals<br />

in UHF (default 0)<br />

NITORD|NITORDER]=nitord In open-shell calculations, the orbitals are reordered after each<br />

iteration to obtain maximum overlap with the orbitals from the<br />

previous iteration. This takes only effect after nitord iterations.<br />

The default is depends on the quality of the starting guess.<br />

NITSH|NITSHIFT=nitsh<br />

NITCL|NITCLOSED=nitcl<br />

NITOCC=nitocc<br />

NITORT|NITORTH=nitort<br />

If the iteration count is smaller than nitsh, the shifts are set to<br />

zero. The default depends on the quality of the starting guess.<br />

TORT—<br />

If the iteration count is smaller than nitcl, only the closed-shell<br />

part of the Fock matrix is used (default nitcl= 0).<br />

Starting with iteration nitocc the occupation pattern is kept fixed.<br />

The default depends on the quality of the starting guess.<br />

The orbitals are reorthonormalized after every nitort iterations.<br />

The default is nitort= 10.<br />

POTFAC=potfac Scale factor for potential energy in first iteration (default 1.0).<br />

Note that in case of a restart the iteration count starts with 3.<br />

16.1.2 Options for the diagonalization method<br />

In calculations with very large basis sets, the diagonalization time becomes a significant fraction<br />

of the total CPU time. This can be reduced using the orbital rotation method as described in R.<br />

Polly, H.-J. Werner, F. R. Manby, and Peter J. Knowles, Mol. Phys. 102, 2311 (2004))<br />

MINROT=minrot<br />

NEXPR=nexpr<br />

DEROT=nexpr<br />

JACOBI=jacobi<br />

If minrot≥ 0, the orbital rotation method is employed. Explicit<br />

diagonalization of the full Fock matrix is performed in the first<br />

minrot iterations and in the last iteration. If minrot=0, a default<br />

is used which depends on the starting guess.<br />

Number of terms used in the exponential expansion of the unitary<br />

orbital transformation matrix (default 4).<br />

Energy gap used in the orbital rotation method. For orbitals<br />

within ±derot hartree of the HOMO orbital energy the Fock<br />

matrix is constructed and diagonalized (default 1.0)<br />

If nonzero, use Jacobi diagonalization.<br />

16.1.3 Options for convergence acceleration methods (DIIS)<br />

For more details, see IPOL directive.<br />

IPTYP=iptyp<br />

IPNIT|DIIS START=ipnit<br />

Interpolation type (default DIIS, see IPOL directive).<br />

First iteration for DIIS interpolation.<br />

IPSTEP|DIIS STEP=ipstep Iteration increment for DIIS interpolation.<br />

MAXDIS|MAXDIIS=maxdis<br />

Max number of Fock matrices used in DIIS interpolation (default<br />

10).

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

Saved successfully!

Ooh no, something went wrong!