SCF with zero iterations

[jvi...@gmail.com]

Dear Dirac Experts,

is it possible in Dirac just to load orbitals (from DFCOEF, DFPCMO or DFACMO), re-compute SCF energy, but do no SCF iterations and let these orbitals to be accepted in further post-HF calculations as converged?

Maybe, I will set ?

.ERGCNV 

1D+10 1D+10

I am interested in this in particular as I have some D17 computed orbitals and would like to compare SCF energy with D19 version, but keep the orbitals the same (no iterations, but SCF energy to be recomputed, not just read-in from DF* file).

Best regards,

Jakub

[Trond SAUE]

HI Jakub,

in each SCF iteration, DIRAC will build the Fock-matrix and then calculate the energy. It then checks for convergence. If convergence, the cycle will exit. If not, DIRAC proceeds to diagonalization of the Fock matrix in orthonormal basis.

This means that if you read in a set of coefficients, DIRAC will calculate the energy associated with these coefficients. You can check orbital energies and the coefficients themselves through the analysis module.

If DIRAC finds that the calculation has converged, it will stop and print the usual information about total and orbital energies and you have your original coefficients.

If DIRAC finds that the calculation has not converged, it will proceed to diagonalization, so you get new coefficients.

For your purpose you can either play with the convergence thresholds, as you suggest, or you can play with the MAXITR keyword.

All the best,

   Trond

--

Trond Saue
Laboratoire de Chimie et Physique Quantiques
UMR 5626 CNRS --- Université Toulouse III-Paul Sabatier
118 route de Narbonne, F-31062 Toulouse, France

Phone : +33/561556361 Fax: +33/561556065
Mail : trond...@irsamc.ups-tlse.fr
Web : https://dirac.ups-tlse.fr/saue
DIRAC : http://www.diracprogram.org/
ESQC : http://www.esqc.org/

Book: Principles and Practices of Molecular Properties: Theory, Modeling, and Simulations