Coulomb interactions diatomic molecules

In Chapter IX, Liang et al. present an approach, termed as the crude Bom-Oppenheimer approximation, which is based on the Born-Oppen-heimer approximation but employs the straightforward perturbation method. Within their chapter they develop this approximation to become a practical method for computing potential energy surfaces. They show that to carry out different orders of perturbation, the ability to calculate the matrix elements of the derivatives of the Coulomb interaction with respect to nuclear coordinates is essential. For this purpose, they study a diatomic molecule, and by doing that demonstrate the basic skill to compute the relevant matrix elements for the Gaussian basis sets. Finally, they apply this approach to the H2 molecule and show that the calculated equilibrium position and foree constant fit reasonable well those obtained by other approaches. 

Like atomic orbitals (AOs), molecular orbitals (MOs) are conveniently described by quantum mechanics theory. Nevertheless, the approach is more complex, because the interaction involves not simply one proton and one electron, as in the case of AOs, but several protons and electrons. For instance, in the simple case of two hydrogen atoms combined in a diatomic molecule, the bulk coulombic energy generated by the various interactions is given by four attractive effects (proton-electron) and two repulsive effects (proton-proton and electron-electron cf figure 1.20) ... 

Figure 1.20 Coulombic interactions in a diatomic molecule composed of hydrogen ions Pa and are the proton and electron of atom A, and pg and Cb are those of atom B.

ICRU Report 55 [36], and references therein. Fast electrons predominantly interact via the coulomb force with the bound electrons of the medium resulting in ionization, i.e., leading to formation of free electrons and residual positive ions. The residual ions can be left in any of a wide variety of final states ranging from their ground ionic state, to states resulting from simultaneous ionization and excitation, and/or dissociation of molecular constituents of the target material. The interaction of a primary electron gp with the diatomic molecule AB can precede via many channels, e.g., ionization can occur by any of the following reaction pathways... 

Recently, detailed molecular pictures of the interfacial structure on the time and distance scales of the ion-crossing event, as well as of ion transfer dynamics, have been provided by Benjamin s molecular dynamics computer simulations [71, 75, 128, 136]. The system studied [71, 75, 136] included 343 water molecules and 108 1,2-dichloroethane molecules, which were separately equilibrated in two liquid slabs, and then brought into contact to form a box about 4 nm long and of cross-section 2.17 nmx2.17 nm. In a previous study [128], the dynamics of ion transfer were studied in a system including 256 polar and 256 nonpolar diatomic molecules. Solvent-solvent and ion-solvent interactions were described with standard potential functions, comprising coulombic and Lennard-Jones 6-12 pairwise potentials for electrostatic and nonbonded interactions, respectively. While in the first study [128] the intramolecular bond vibration of both polar and nonpolar solvent molecules was modeled as a harmonic oscillator, the next studies [71,75,136] used a more advanced model [137] for water and a four-atom model, with a united atom for each of two... 

Equation (6.13) was derived on a model of covalent bonds between nearest neighbors. It is not strictly applicable to ionic solids. The repulsion part of the potential energy must be similar for ionic and covalent cases, but the attraction part for ionic solids must also include the sum of the coulombic interactions with the remainder of the lattice. In effect, the number of bonds is increased. To see the magnitude of this effort, compare Equations (6.7) and (6.8) with their counterparts for a diatomic molecule, or ion-pair. 

The interaction energy represents the binding energy of a diatomic molecule. It is usually known as a function of inter-nuclear distance from an experimentally determined equation, as, for example, a Morse curve. If we have some notion as to the relative proportions of and a/3, we can estimate all the coulomb and exchange integrals for the pairs and substitute into an expression such as Eq. (98) to obtain E for more complex systems. By judicious variation of the so-called coulomb fraction one can get a fair agreement of E with experimental determinations thereof. The coulomb fractions ordinarily need not be varied more than from about 1/10 to 3/10, which is close to the theoretical calculation for hydrogen atoms. 

In sohds, then, the forces which tend to impose order on the individual chemical units may be van der WaaJs forces, as with the inert gases, with elements which form diatomic molecules, and with organic substances they may be Coulomb forces between ions, as in many simple salts and with the positive ion-electron interaction in metals and, finally, they may be the covalencies binding atom to atom throughout extended regions of space. 

A diatomic homonuclear molecule, origin of the BFCS in the centre of the molecule, potential energy of the Coulombic interactions equals V. The total non-relativistic Hamiltonian is equal to ... 

Harmonium represents the two-electron Hooke atom. A Hooke diatomic molecule means two heavy particles (nuclei) interacting by Coulomb forces. The same is true with electrons, but the heavy particle-light particle interactions are harmonic. 

Furthermore their incidence are very dependent upon the nature of the properties [52,53], Due to eomputer limitations, basis sets cannot be extended indefinitely and direct numerical evaluations seem the ultimate solution for molecules [54], In position space this is a viable alternative for diatomic molecules [55,56], but it cannot be extended easily to polyatomic systems. Formulated in momentum space, the HF equations have not explicit solutions and the difficulties to express them in terms of basis functions are analogous to those encountered in the r-space. However the momentum space HF equations give way to numerical approaches in which Coulombic interactions become tractable even for polyatomic molecules [7] among other advantages, these equations, Eqs. 13 and 20, do not require coordinate systems adapted to the geometry of the molecules to remove Coulombic singularities. In both equations the only singular contribution comes from the q 2 factor. 

PM3, developed by James J.P. Stewart, is a reparameterization of AMI, which is based on the neglect of diatomic differential overlap (NDDO) approximation. NDDO retains all one-center differential overlap terms when Coulomb and exchange integrals are computed. PM3 differs from AMI only in the values of the parameters. The parameters for PM3 were derived by comparing a much larger number and wider variety of experimental versus computed molecular properties. Typically, non-bonded interactions are less repulsive in PM3 than in AMI. PM3 is primarily used for organic molecules, but is also parameterized for many main group elements. 

---