Computational chemistry Hamiltonians

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. 

Formulation of the methods of computational chemistry is reasonably straightforward. The hard work comes in its implementation. We begin with the electronic Hamiltonian for a molecule, a generalization of that for a many-electron atom given in Eq (9.2) ... 

In recent years, dramatic advances in computational power combined with the marketing of packaged computational chemistry codes have allowed quantum chemical calculations to become fairly routine in both prediction and verification of experimental observations. The 1998 Nobel Prize in Chemistry reflected this impact by awarding John A. Pople a shared prize for his development of computational methods in quantum chemistry. The Hartree-Fock approximation is a basic approach to the quantum chemical problem described by the Schrodinger equation, equation (3.10), where the Hamiltonian (//) operating on the wavefunction OP) yields the energy (E) multiplied by the wavefunction. 

Inclusion of Pauli s exclusion principle leads to the standard methods of ab initio computational chemistry. Within these methods, molecular systems containing the same nuclei and the same number of electrons, but having a different total electronic spin, can roughly speaking be said to be different systems. Thus, matrix elements of the Hamiltonian between Slater determinants corresponding to different spin states will all be zero, and they will not interact or mix at all. The wavefunctions obtained will be pure spin states. ... 

Operators that result from a DK transformation are directly given in the momentum representation. Hess et al. [29,31] developed a very efficient strategy to evaluate the corresponding matrix elements in a basis set representation it employs the eigenvectors of the operator as approximate momentum representation [29,31]. In practice, the two-component DK Hamiltonian is built of matrix representations of the three operators p, V, pVp + id(pV x p). This Douglas-Kroll-Hess (DKH) approach became one of the most successful two-component tools of relativistic computational chemistry [16,74]. In particular, many applications showed that the second-order operator 2 Is variationally stable [10,13,14,31,75,76,87]. 

The methods today more in use in computational chemistry belong to the ASC, MPE, GB and FD families. For every type of method there are now QM versions, many thus far limited to infinite isotropic distributions. Several of those methods may introduce, via the effective Hamiltonian or with less formal procedures, solute-solvent interaction effects of non electrostatic origin. 

We will start with a description of FDE and its ability to generate diabats and to compute Hamiltonian matrix elements—the EDE-ET method (ET stands for Electron Transfer). In the subsequent section, we will present specific examples of FDE-ET computations to provide the reader with a comprehensive view of the performance and applicability of FDE-ET. After FDE has been treated, four additional methods to generate diabatic states are presented in order of accuracy CDFT, EODFT, AOM, and Pathways. In order to output a comprehensive presentation, we also describe those methods in which wavefunctions methods can be used, in particular GMH and other adiabatic-to-diabatic diabatization methods. Finally, we provide the reader with a protocol for running FDE-ET calculations with the only available implementation of the method in the Amsterdam Density Functional software [51]. In closing, we outline our concluding remarks and our vision of what the future holds for the field of computational chemistry applyed to electron transfer. 

K. Hirao. Relativistic Multireference Perturbation Theory Complete Active-Space Second-Order Perturbation Theory (CASPT2) With The Four-Component Dirac Hamiltonian. In Radiation Induced Molecular Phenomena in Nucleic Adds, Volume 4 of Challenges and Advances in Computational Chemistry and Physics, p. 157-177. Springer, 2008. 

Moller-Plesset second-order perturbation theory (MP2) is a common method used in computational chemistry to include electron correlation as an extension to Hartree-Fock (HF) theory which neglects Coulomb correlation and thus also misses all dispersion effects. The perturbation is the difference between the Fock-operator and the exact electronic Hamiltonian. 

The determination of the ground state energy and the ground state electron density distribution of a many-electron system in a fixed external potential is a problem of major importance in chemistry and physics. For a given Hamiltonian and for specified boundary conditions, it is possible in principle to obtain directly numerical solutions of the Schrodinger equation. Even with current generations of computers, this is not feasible in practice for systems of large total number of electrons. Of course, a variety of alternative methods, such as self-consistent mean field theories, also exist. However, these are approximate. 

Most chemists are quite comfortable thinking of chemical structure and reactivity in terms of valence bond notions - die resonance structures so often invoked in organic chemistry are one example of this phenomenon - so this approach has conceptual appeal. From a computational standpoint, the issue is how to derive a Hamiltonian operator tliat will act on VB wave functions so as to deliver useful energies. 

---