Effective bond Hamiltonians

The matrix elements of the effective bond Hamiltonians are defined as (with the MINDO/3 parameterization for the Hamiltonian taken for the sake of definiteness)  

Pseudospin representation and the perturbative estimates of the bond-geminal ESVs. To provide the required explanation, we notice that the effective Hamiltonians for the bond geminals can be represented as a sum of the unperturbed part which, when diagonalized yields invariant, i.e. exactly transferable, values of the ESVs, and of a perturbation responsible for the specificity of electronic structure for different chemical compositions and environments of the bond. 

Pseudospin operator of the bond geminal. Let us introduce a pseudospin operator Tm corresponding to the pseudospin value rTO = 1. The matrices of its components in the basis of the configurations defining the geminal are given by  

The configurations corresponding to (fzm) = 1 are the ionic ones with both electrons located on the same end of the chemical bond (right or left, respectively). In terms of the pseudospin operator, the elements of the density matrices in eqs. (2.78), (2.81) can be presented as follows  

The effective bond Hamiltonians also can be rewritten in terms of the pseudospin operators. Indeed, the effective Hamiltonian eq. (3.1) for each of the bond geminals can be presented in the form  

---

Using the separation of the effective Hamiltonian into the unperturbed part and the perturbation, the total ionic contribution to the geminal is calculated exactly (vari-ationally). Only the bond polarity needs to be estimated perturbatively in the linear response approximation, but now the correlated ground state of the symmetric effective bond Hamiltonian is taken for evaluating the response function. In this context, it is convenient to use a dimensionless bond asymmetry parameter ... 

As mentioned previously, the density ESVs must be obtained from the effective bond Hamiltonian eq. (4.1). In terms of the geminal amplitudes, the ESVs are given by eq. (2.78). To get the required direct estimates of the ESVs, we use again the projection operator technique. In terms of the geminal amplitudes (subject to the normalization condition) the projection operator upon die ground state of a geminal has the form ... 

The last relation for the product xy is not an independent equation but it must be inserted into that for y and the system becomes one for x and y. Solving this system will be equivalent to solving the original 3x3 eigenvalue problem for the effective bond Hamiltonian. In a perturbative manner we get for the first order approximation ... 

The chemical bonds of a molecule are not isolated, but they are in the electrostatic field of each other—an effect which can be accounted for even in strictly localized models (Malrieu 1977). This can be done by defining the effective bond Hamiltonians analogously to the group function theory (McWeeny 1959, Mehler 1977, 1981). The effective Hamiltonians will be used to write the local Schrodinger equations for each bond, which determine the expansion coefficients of the SLGs ... 

In this approximate treatment, one simply freezes the nonreactive or spectator bond. The 6D Hamiltonian in this case becomes the effective 5D Hamiltonian by eliminating the kinetic energy operator for the r2 coordinate and fixed the r2 distance in the potential energy surface. The 5D potential is simply given by... 

In analogy with the molecular mechanics method, whose force fields derive from intramolecular and intermolecular potentials, the effective electronic hamiltonian can be comprised of intra and intermolecular matrix elements. The intramolecular matrix elements of the effective tight-binding hamiltonian must describe principally the energies of the chemical bondings. In the EHT method this is done by the matrix elements // = + H )Sij/2, where the diagonal elements // . corre-... 

Figure C3.2.18.(a) Model a-helix, (b) hydrogen bonding contacts in tire helix, and (c) schematic representation of tire effective Hamiltonian interactions between atoms in tire protein backbone. From [23].

For this reason, there has been much work on empirical potentials suitable for use on a wide range of systems. These take a sensible functional form with parameters fitted to reproduce available data. Many different potentials, known as molecular mechanics (MM) potentials, have been developed for ground-state organic and biochemical systems [58-60], They have the advantages of simplicity, and are transferable between systems, but do suffer firom inaccuracies and rigidity—no reactions are possible. Schemes have been developed to correct for these deficiencies. The empirical valence bond (EVB) method of Warshel [61,62], and the molecular mechanics-valence bond (MMVB) of Bemardi et al. [63,64] try to extend MM to include excited-state effects and reactions. The MMVB Hamiltonian is parameterized against CASSCF calculations, and is thus particularly suited to photochemistry. 

The first Hamiltonian was used in the early simulations on two-dimensional glass-forming lattice polymers [42] the second one is now most frequently used in two and three dimensions [4]. Just to illustrate the effect of such an energy function, which is given by the bond length, Fig. 10 shows two different states of a two-dimensional polymer melt and, in part. 

Hydrogen abstraction reactions potential surfaces for, 25-26,26,41 resonance structures for, 24 Hydrogen atom, 2 Hydrogen bonds, 169,184 Hydrogen fluoride, 19-20, 20,22-23 Hydrogen molecules, 15-18 energy of, 11,16,17 Hamiltonian for, 4,15-16 induced dipoles, 75,125 lithium ion effect on, 12... 

The QM/MM interactions (Eqm/mm) are taken to include bonded and non-bonded interactions. For the non-bonded interactions, the subsystems interact with each other through Lennard-Jones and point charge interaction potentials. When the electronic structure is determined for the QM subsystem, the charges in the MM subsystem are included as a collection of fixed point charges in an effective Hamiltonian, which describes the QM subsystem. That is, in the calculation of the QM subsystem we determine the contributions from the QM subsystem (Eqm) and the electrostatic contributions from the interaction between the QM and MM subsystems as explained by Zhang et al. [13],... 

Figure 4-2. Computed potential energy surface from (A) ab initio valence-bond self-consistent field (VB-SCF) and (B) the effective Hamiltonian molecular-orbital and valence-bond (EH-MOVB) methods for the S 2 reaction between HS- and CH3CI...

---