Electronic Hamiltonian function

In a diabatic representation, the electronic wave functions are no longer eigenfunctions of the electronic Hamiltonian. The aim is instead that the functions are so chosen that the (nonlocal) non-adiabatic coupling operator matrix, A in Eq. (52), vanishes, and the couplings are represented by (local) potential operators. The nuclear Schrddinger equation is then written... 

Yarkoni [108] developed a computational method based on a perturbative approach [109,110], He showed that in the near vicinity of a conical intersection, the Hamiltonian operator may be written as the sum a nonperturbed Hamiltonian Hq and a linear perturbative temr. The expansion is made around a nuclear configuration Q, at which an intersection between two electronic wave functions takes place. The task is to find out under what conditions there can be a crossing at a neighboring nuclear configuration Qy. The diagonal Hamiltonian matrix elements at Qy may be written as... 

Until now we have implicitly assumed that our problem is formulated in a space-fixed coordinate system. However, electronic wave functions are naturally expressed in the system bound to the molecule otherwise they generally also depend on the rotational coordinate 4>. (This is not the case for E electronic states, for which the wave functions are invariant with respect to (j> ) The eigenfunctions of the electronic Hamiltonian, v / and v , computed in the framework of the BO approximation ( adiabatic electronic wave functions) for two electronic states into which a spatially degenerate state of linear molecule splits upon bending. 

Since the form of the electronic wave functions depends also on the coordinate p (in the usual, parametric way), the matrix elements (21) are functions of it too. Thus it looks at first sight as if a lot of cumbersome computations of derivatives of the electronic wave functions have to be carried out. In this case, however, nature was merciful the matrix elements in (21) enter the Hamiltonian matrix weighted with the rotational constant A, which tends to infinity when the molecule reaches linear geometry. This means that only the form of the wave functions, that is, of the matrix elements in (21), in the p 0 limit are really needed. In the above mentioned one-elecbon approximation... 

We write them as i / (9) to shess that now we use the space-fixed coordinate frame. We shall call this basis diabatic, because the functions (26) are not the eigenfunction of the electronic Hamiltonian. The matrix elements of are... 

A convenience of electronic basis functions (53) is that they reduce at infinitesimal-amplitude bending to (28) with the same meaning of the angle 9 we may employ these asymptotic forms in the computation of the matrix elements of the kinetic energy operator and in this way avoid the necessity of carrying out calculations of the derivatives of the electronic wave functions with respect to the nuclear coordinates. The electronic part of the Hamiltonian is represented in the basis (53) by... 

Let us examine a special but more practical case where the total molecular Hamiltonian, H, can be separated to an electronic part, W,.(r,s Ro), as is the case in the usual BO approximation. Consequendy, the total molecular wave function fl(R, i,r,s) is given by the product of a nuclear wave function X uc(R, i) and an electronic wave function v / (r, s Ro). We may then talk separately about the permutational properties of the subsystem consisting of electrons, and the subsystemfs) formed of identical nuclei. Thus, the following commutative laws Pe,Hg =0 and =0 must be satisfied X =... 

In an ab initio method, all the integrals over atom ic orbital basis function s are com puted and the bock in atrix of th e SCK com puta-tiori is formed (equation (6 1) on page 225) from the in tegrals. Th e Kock matrix divides inui two parts the one-electron Hamiltonian matrix, H, and the two-electron matrix, G, with the matrix ele-m en ts... 

HyperChem s semi-empirical calculations solve (approximately) the Schrodinger equation for this electronic Hamiltonian leading to an electronic wave function I eiecW for the electrons ... 

For this simple Hamiltonian, which involves the sum of one-electron Hamiltonians, we can use a wave function of the form... 

The Hamiltonian function for an electron in a constant magnetic field of strength H parallel to the s axis is1... 

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. 

It is also of interest to study the "inverse" problem. If something is known about the symmetry properties of the density or the (first order) density matrix, what can be said about the symmetry properties of the corresponding wave functions In a one electron problem the effective Hamiltonian is constructed either from the density [in density functional theories] or from the full first order density matrix [in Hartree-Fock type theories]. If the density or density matrix is invariant under all the operations of a space CToup, the effective one electron Hamiltonian commutes with all those elements. Consequently the eigenfunctions of the Hamiltonian transform under these operations according to the irreducible representations of the space group. We have a scheme which is selfconsistent with respect to symmetty. 

The electronically adiabatic wave functions v(/f ad(r q ) are defined as eigenfunctions of the electronic Hamiltonian Hel with electronically adiabatic potential energies ad(qjJ as their eigenvalues ... 

In practice, each CSF is a Slater determinant of molecular orbitals, which are divided into three types inactive (doubly occupied), virtual (unoccupied), and active (variable occupancy). The active orbitals are used to build up the various CSFs, and so introduce flexibility into the wave function by including configurations that can describe different situations. Approximate electronic-state wave functions are then provided by the eigenfunctions of the electronic Hamiltonian in the CSF basis. This contrasts to standard HF theory in which only a single determinant is used, without active orbitals. The use of CSFs, gives the MCSCF wave function a structure that can be interpreted using chemical pictures of electronic configurations [229]. An interpretation in terms of valence bond structures has also been developed, which is very useful for description of a chemical process (see the appendix in [230] and references cited therein). 

Further, the time-independent electronic basis functions are taken to be the eigenfunctions of the electronic Hamiltonian,... 

---