Coulomb interactions

Forces between macroscopic objects result from a complex interplay of the interaction between molecules in the two objects and the medium separating 

1) Johannes Diderik van der Waals, 1837—1923. Dutch physicist, professor in Amsterdam. Nobel Prize in Physics, 1910. 

The potential energy between two electrical charges that are a distance D apart is 

This is 56 times higher than the thermal energy keT = 4.12 x 10 J at room temperature. 

---

For dilute solutions, the Debye-Huckel equation by calculations based on these Coulombic interactions is... 

The influence of electrical charges on surfaces is very important to their physical chemistry. The Coulombic interaction between charged colloids is responsible for a myriad of behaviors from the formation of opals to the stability of biological cells. Although this is a broad subject involving both practical application and fundamental physics and chemistry, we must limit our discussion to those areas having direct implications for surface science. 

There are tliree important varieties of long-range forces electrostatic, induction and dispersion. Electrostatic forces are due to classical Coulombic interactions between the static charge distributions of the two molecules. They are strictly pairwise additive, highly anisotropic, and can be either repulsive or attractive. 

Each electron in the system is assigned to either molecule A or B, and Hamiltonian operators and for each molecule defined in tenns of its assigned electrons. The unperturbed Hamiltonian for the system is then 0 = - A perturbation XH consists of tlie Coulomb interactions between the nuclei and... 

The theory of strong electrolytes due to Debye and Htickel derives the exact limiting laws for low valence electrolytes and introduces the idea that the Coulomb interactions between ions are screened at finite ion concentrations. 

The main difficulty in these simulations is the long-range nature of the Coulomb interactions, since both mirror-plane images and real charges must be included, and the finite nature of the simulated volume must also be mchided. A more detailed discussion is given by Benjamin [29], and the following conclusions have been reached. 

Diflfiision-controlled reactions between ions in solution are strongly influenced by the Coulomb interaction accelerating or retarding ion diffiision. In this case, die dififiision equation for p concerning motion of one reactant about the other stationary reactant, the Debye-Smoluchowski equation. 

Wlien the potential consists of electron-electron and electron-nucleus Coulombic interactions,... 

Many realistic simulations will involve the Coulomb interaction between charges, which decreases with... 

Flere we distinguish between nuclear coordinates R and electronic coordinates r is the single-particle kinetic energy operator, and Vp is the total pseudopotential operator for the interaction between the valence electrons and the combined nucleus + frozen core electrons. The electron-electron and micleus-micleus Coulomb interactions are easily recognized, and the remaining tenu electronic exchange and correlation... 

Atomistically detailed models account for all atoms. The force field contains additive contributions specified in tenns of bond lengtlis, bond angles, torsional angles and possible crosstenns. It also includes non-bonded contributions as tire sum of van der Waals interactions, often described by Lennard-Jones potentials, and Coulomb interactions. Atomistic simulations are successfully used to predict tire transport properties of small molecules in glassy polymers, to calculate elastic moduli and to study plastic defonnation and local motion in quasi-static simulations [fy7, ( ]. The atomistic models are also useful to interiDret scattering data [fyl] and NMR measurements [70] in tenns of local order. 

The f operators are the usual kinetic energy operators, and the potential energy V(r,R) includes all of the Coulomb interactions ... 

In the work of King, Dupuis, and Rys [15,16], the mabix elements of the Coulomb interaction term in Gaussian basis set were evaluated by solving the differential equations satisfied by these matrix elements. Thus, the Coulomb matrix elements are expressed in the form of the Rys polynomials. The potential problem of this method is that to obtain the mabix elements of the higher derivatives of Coulomb interactions, we need to solve more complicated differential equations numerically. Great effort has to be taken to ensure that the differential equation solver can solve such differential equations stably, and to... 

Nevertheless, the examination of the applicability of the crude BO approximation can start now because we have worked out basic methods to compute the matrix elements. With the advances in the capacity of computers, the test of these methods can be done in lower and lower cost. In this work, we have obtained the formulas and shown their applications for the simple cases, but workers interested in using these matrix elements in their work would find that it is not difficult to extend our results to higher order derivatives of Coulomb interaction, or the cases of more-than-two-atom 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. 

Separation of Short- and Long-Range Ewald and PPPM Methods If we split the total Coulomb interaction in a short- and a long-range contribution, chosen to be smooth functions of the distance, the two... 

Abstract. Molecular dynamics (MD) simulations of proteins provide descriptions of atomic motions, which allow to relate observable properties of proteins to microscopic processes. Unfortunately, such MD simulations require an enormous amount of computer time and, therefore, are limited to time scales of nanoseconds. We describe first a fast multiple time step structure adapted multipole method (FA-MUSAMM) to speed up the evaluation of the computationally most demanding Coulomb interactions in solvated protein models, secondly an application of this method aiming at a microscopic understanding of single molecule atomic force microscopy experiments, and, thirdly, a new method to predict slow conformational motions at microsecond time scales. 

Fig. 1. Structure adapted hierarchical description of Coulomb interactions in biological macromolecules. Filled circles (level 0) represent atoms, structural units (li vel 1) are surrounded by a single-line border, and clusters (level 2) are surrounded by a double-line border.

Task 4 Explicitly calculate the Coulomb interactions between atoms which are closer than about 10 A. 

N is the number of point charges within the molecule and Sq is the dielectric permittivity of the vacuum. This form is used especially in force fields like AMBER and CHARMM for proteins. As already mentioned, Coulombic 1,4-non-bonded interactions interfere with 1,4-torsional potentials and are therefore scaled (e.g., by 1 1.2 in AMBER). Please be aware that Coulombic interactions, unlike the bonded contributions to the PEF presented above, are not limited to a single molecule. If the system under consideration contains more than one molecule (like a peptide in a box of water), non-bonded interactions have to be calculated between the molecules, too. This principle also holds for the non-bonded van der Waals interactions, which are discussed in Section 7.2.3.6. 

Many problems in force field investigations arise from the calculation of Coulomb interactions with fixed charges, thereby neglecting possible mutual polarization. With that obvious drawback in mind, Ulrich Sternberg developed the COSMOS (Computer Simulation of Molecular Structures) force field [30], which extends a classical molecular mechanics force field by serai-empirical charge calculation based on bond polarization theory [31, 32]. This approach has the advantage that the atomic charges depend on the three-dimensional structure of the molecule. Parts of the functional form of COSMOS were taken from the PIMM force field of Lindner et al., which combines self-consistent field theory for r-orbitals ( nr-SCF) with molecular mechanics [33, 34]. 

Figure 6.25 reprinted from Chemical Physics Letters, 196, Ding H-Q, N Karasawa and W A Goddard III, T he Reduced Cell Multipole Method for Coulomb Interactions in Periodic Systems with Million-Atom Unit Cells, 6-10, 1992, with permission of Elsevier Science. 

---