Potentials Coulomb

The attractive part of the bonding force is the Coulomb potential, which results from the electrostatic force between two charged particles. 

If the ith and/th ions have the same charge, the potential is positive and the force repulsive if the charges are opposite, the force is attractive. Note the 4ttsq in the denominator. We will use rationalized meter, kilogram, second (MKS) units throughout. 

Another very common case is a potential that behaves as l/ r, known as the Coulomb potential, from the electrostatic interaction between particles with electrical charges at distance r. We will discuss this case for the simplest physical system where it applies, the hydrogen atom. For simplicity, we take the proton fixed at the origin of the coordinate system and the electron at r. The hamiltonian for this system takes the form 

Since in the above equation the left-hand side is exclusively a function of r while the right-hand side is exclusively a function oiO,(p, they each must be equal to a constant, which we denote by A, giving rise to the following two differential equations  

We consider the equation for Y(9, f) first. This equation is solved by the functions 

It is a straightforward exercise to show from these expressions that 

f) = i(i + i)h Yi ie,ci ), Ljue,(P) = mnYue,(P) (b.36) as might have been expected from our earlier identification of the quantity + 1) with the square of the angular momentum. This is another example of simultaneous eigenfunctions of two operators, which according to our earlier discussion must commute [L, LJ = 0. Thus, the spherical harmonics determine the angular momentum I of a state and its z component, which is equal to mh. 

---

The miderstanding of the quantum mechanics of atoms was pioneered by Bohr, in his theory of the hydrogen atom. This combined the classical ideas on planetary motion—applicable to the atom because of the fomial similarity of tlie gravitational potential to tlie Coulomb potential between an electron and nucleus—with the quantum ideas that had recently been introduced by Planck and Einstein. This led eventually to the fomial theory of quaiitum mechanics, first discovered by Heisenberg, and most conveniently expressed by Schrodinger in the wave equation that bears his name. 

Coulomb potential multiplied by -p. The graphical representation of the virial coefficients in temis of Mayer/ -bonds can now be replaced by an expansion in temis ofy bonds and Coulomb bonds ). 

This has the fomi of a second virial coefficient in which the Debye screened potential has replaced the Coulomb potential. Expressions for the other excess themiodynamic properties are easily derived. 

The solutions to this approximation are obtained numerically. Fast Fourier transfonn methods and a refomuilation of the FINC (and other integral equation approximations) in tenns of the screened Coulomb potential by Allnatt [M are especially useful in the numerical solution. Figure A2.3.12 compares the osmotic coefficient of a 1-1 RPM electrolyte at 25°C with each of the available Monte Carlo calculations of Card and Valleau [ ]. 

Where Coulomb potential of nucleus At atomic shell electrons... 

Many-body problems wnth RT potentials are notoriously difficult. It is well known that the Coulomb potential falls off so slowly with distance that mathematical difficulties can arise. The 4-k dependence of the integration volume element, combined with the RT dependence of the potential, produce ill-defined interaction integrals unless attractive and repulsive mteractions are properly combined. The classical or quantum treatment of ionic melts [17], many-body gravitational dynamics [18] and Madelung sums [19] for ionic crystals are all plagued by such difficulties. 

C) All mean-field models of electronic. structure require large corrections. Essentially all ab initio quantum chemistry approaches introduce a mean field potential F that embodies the average interactions among the electrons. The difference between the mean-field potential and the true Coulombic potential is temied [20] the "fluctuationpotentiar. The solutions Ef, to the true electronic... 

The magnitude and shape of such a mean-field potential is shown below [21] in figure B3.1.4 for the two 1 s electrons of a beryllium atom. The Be nucleus is at the origin, and one electron is held fixed 0.13 A from the nucleus, the maximum of the Is orbital s radial probability density. The Coulomb potential experienced by the second electron is then a function of the second electron s position along the v-axis (coimecting the Be nucleus and the first electron) and its distance perpendicular to the v-axis. For simplicity, this second electron... 

Coulomb potential felt by a 2p orbital s electron at a point r in the ls 2s 2p 2p configuration description of the carbon atom is ... 

In the connnonly used atomic sphere approximation (ASA) [79], the density and the potential of the crystal are approximated as spherically synnnetric within overlapping imifiBn-tin spheres. Additionally, all integrals, such as for the Coulomb potential, are perfonned only over the spheres. The limits on the accuracy of the method imposed by the ASA can be overcome with the fiill-potential version of the LMTO (FP-LMTO)... 

If V is the total Coulombic potential between all the nuclei and electrons in the system, then, in the absence of any spin-dependent terms, the electronic Hamiltonian is given by... 

In modem quantum chemistry packages, one can obtain moleculai basis set at the optimized geometry, in which the wave functions of the molecular basis are expanded in terms of a set of orthogonal Gaussian basis set. Therefore, we need to derive efficient fomiulas for calculating the above-mentioned matrix elements, between Gaussian functions of the first and second derivatives of the Coulomb potential ternis, especially the second derivative term that is not available in quantum chemistry packages. Section TV is devoted to the evaluation of these matrix elements. 

Next, we shall consider four kinds of integrals. The first is the expectation value of the Coulomb potential by one nucleus for one of the primitive basis function centered at that nucleus. The second is the expectation value of the Coulomb potential by one nucleus for one of the primitive basis function centered at a different point (usually another nucleus). Then, we will consider the matrix element of a Coulomb term between two primitive basis functions at different centers. The third case is when one basis function is centered at the nucleus considered. The fourth case is when both basis functions are not centered at that nucleus. By that we mean, for two Gaussian basis functions defined in Eqs. (73) and (74), we are calculating... 

In fact, the Coulomb integrals discussed in Section IV.C are available in contemporary quantum chemistry packages. We do not really need to develop our own method to calculate them. However, it is necessary to master the algebra so that we can calculate the matrix elements of the derivatives of the Coulomb potential. In the following, we shall demonstrate the evaluation of these matrix elements. 

Efficient Integrators for Systems with Coulomb Potentials... 

This represents an attractive (Coulombic) potential coupled with a repulsive soft wall, relative to a plane situated just below the rigid body. The rigid body is repeated drawn toward the plane, then repelled sharply from the wall. 

Coulombic potential energy is calculated by modification and fitting of some form of Coulomb s equation... 

There is a very convenient way of writing the Hamiltonian operator for atomic and molecular systems. One simply writes a kinetic energy part — for each election and a Coulombic potential Z/r for each interparticle electrostatic interaction. In the Coulombic potential Z is the charge and r is the interparticle distance. The temi Z/r is also an operator signifying multiply by Z r . The sign is - - for repulsion and — for atPaction. 

The expression for the Coulombic potential energy e/4neo can be canied through the entire derivation in Exercise 6-3 to anive at Eq. (6-17). Show that this is so. 

---