Self-consistent field calculations, computer codes

If we except the Density Functional Theory and Coupled Clusters treatments (see, for example, reference [1] and references therein), the Configuration Interaction (Cl) and the Many-Body-Perturbation-Theory (MBPT) [2] approaches are the most widely-used methods to deal with the correlation problem in computational chemistry. The MBPT approach based on an HF-SCF (Hartree-Fock Self-Consistent Field) single reference taking RHF (Restricted Hartree-Fock) [3] or UHF (Unrestricted Hartree-Fock ) orbitals [4-6] has been particularly developed, at various order of perturbation n, leading to the widespread MPw or UMPw treatments when a Moller-Plesset (MP) partition of the electronic Hamiltonian is considered [7]. The implementation of such methods in various codes and the large distribution of some of them as black boxes make the MPn theories a common way for the non-specialist to tentatively include, with more or less relevancy, correlation effects in the calculations. 

The methodology that uses the dielectric model is instead the simpler and in principle the more suitable for the study of chemical reactions involving large molecular systems. In 1998, Amovilli et al [13] developed a computer code in which the solvent reaction field, including all the basic solute-solvent interactions, has been considered for Complete Active Space Self Consistent Field (CASSCF) calculations. 

Abstract. The direct self-consistent-field (SCF) method recalculates all two-electron integrals each time they are needed in an SCF calculation. This perspective article discusses how the original paper on direct SCF by Almlof et al. developed the principles by which this could be made efficient and thereby provided an example of the semantic approach to computational chemistry in which algorithm development and coding are not compartmentalized. 

With respect to standard molecular-cluster techniques, this approach has some attractive features explicit reference is made to the HF LCAO periodic solution for the unperturbed (or perfect) host crystal. In particular, the self-embedding-consistent condition is satisfied, that is, in the absence of defects, the electronic structure in the cluster region coincides with that of the perfect host crystal there is no need to saturate dangling bonds the geometric constraints and the Madelung field of the environment are automatically included. With respect to the supercell technique, this approach does not present the problem of interaction between defects in different supercells, allows a more flexible definition of the cluster subspace, and permits the study of charged defects. The perturbed-cluster approach is implemented in the computer code EMBEDOl [703] and applied in the calculations of the point defects both in the bulk crystal, [704] and on the surface [705]. The difficulties of this approach are connected with the lattice-relaxation calculations. 

---