Potential energy computation

How is electronic potential energy computed Electrons, which are more than three orders of magnitude lighter than nuclei, are too small for classical mechanics calculations. Electronic energy must... 

Let us briefly state where we stand, valence denotes occupied valence orbitals only and is simply represented by the straight sum of their pertinent kinetic and potential energies computed over the entire coordinate space. in contrast, described by I (T + 2V ), represents all occupied orbitals but integrated only over specifled (core and valence) regions of space. The relationship between the two, Eq. (4.36), is surprisingly simple considering the basic differences between the two models. 

Fig. 6.1. Potential energy computed in the RHF, UHF, HL and MC-SCF HF-HL approximations and experimental data. For H2 we include also the potential for a six configuration MC-SCF (see text).

Recent CNDO/2 computations give also similar results (Warner, D. and Steward, E. C. /. Mol. Struct. 25,403 (1975)). These authors also show that the magnitude of the stabilization of the folded form is in agreement with a classical coulomb potential energy computation. 

Another statistical mechanical approach makes use of the radial distribution function g(r), which gives the probability of finding a molecule at a distance r from a given one. This function may be obtained experimentally from x-ray or neutron scattering on a liquid or from computer simulation or statistical mechanical theories for model potential energies [56]. Kirkwood and Buff [38] showed that for a given potential function, U(r)... 

The theory coimecting transport coefficients with the intemiolecular potential is much more complicated for polyatomic molecules because the internal states of the molecules must be accounted for. Both quantum mechanical and semi-classical theories have been developed. McCourt and his coworkers [113. 114] have brought these theories to computational fruition and transport properties now constitute a valuable test of proposed potential energy surfaces that... 

Statistical mechanical theory and computer simulations provide a link between the equation of state and the interatomic potential energy functions. A fluid-solid transition at high density has been inferred from computer simulations of hard spheres. A vapour-liquid phase transition also appears when an attractive component is present hr the interatomic potential (e.g. atoms interacting tlirough a Leimard-Jones potential) provided the temperature lies below T, the critical temperature for this transition. This is illustrated in figure A2.3.2 where the critical point is a point of inflexion of tire critical isothemr in the P - Vplane. 

In this section we present several numerical teclmiques that are conmronly used to solve the Sclirodinger equation for scattering processes. Because the potential energy fiinctions used in many chemical physics problems are complicated (but known to reasonable precision), new numerical methods have played an important role in extending the domain of application of scattering theory. Indeed, although much of the fomial development of the previous sections was known 30 years ago, the numerical methods (and computers) needed to put this fomialism to work have only been developed since then. 

Tellinghuisen J A 1974 A fast quadrature method for computing diatomic RKR potential energy curves Comput. Phys. Commun. 6 221-8... 

Jensen F 1994 Transition structure modeling by intersecting potential energy surfaces J. Comput. Chem. 15 1199... 

Static properties of some molecules ([193,277-280]). More recently, pairs of ci s have been studied [281,282] in greater detail. These studies arose originally in connection with a ci between the l A and 2 A states found earlier in computed potential energy surfaces for C2H in symmetry [278]. Similar ci s appear between the potential surfaces of the two lowest excited states A2 and B2 iit H2S or of 82 and A in Al—H2 within C2v symmetry [283]. A further, closely spaced pair of ci s has also been found between the 3 A and 4 A states of the molecule C2H. Here the separation between the twins varies with the assumed C—C separation, and they can be brought into coincidence at some separation [282]. 

The full dynamical treatment of electrons and nuclei together in a laboratory system of coordinates is computationally intensive and difficult. However, the availability of multiprocessor computers and detailed attention to the development of efficient software, such as ENDyne, which can be maintained and debugged continually when new features are added, make END a viable alternative among methods for the study of molecular processes. Eurthemiore, when the application of END is compared to the total effort of accurate determination of relevant potential energy surfaces and nonadiabatic coupling terms, faithful analytical fitting and interpolation of the common pointwise representation of surfaces and coupling terms, and the solution of the coupled dynamical equations in a suitable internal coordinates, the computational effort of END is competitive. 

The combination is in this case an out-of-phase one (Section I). This biradical was calculated to be at an energy of 39.6 kcal/mol above CHDN (Table ni), and to lie in a real local minimum on the So potential energy surface. A normal mode analysis showed that all frequencies were real. (Compare with the prebenzvalene intermediate, discussed above. The computational finding that these species are bound moieties is difficult to confimi experimentally, as they are highly reactive.)... 

An alternative approximation scheme, also proposed by Bom and Oppenheimer [5-7], employed the straightforward perturbation method. To tell the difference between these two different BO approximation, we call the latter the crude BOA (CBOA). A main purpose of this chapter is to study the original BO approximation, which is often referred to as the crude BO approximation and to develop this approximation into a practical method for computing potential energy suifaces of molecules. 

Figure 3. Low-energy vibronic spectrum in a. 11 electronic state of a linear triatomic molecule, computed for various values of the Renner parameter e and spin-orbit constant Aso (in cm ). The spectrum shown in the center of figure (e = —0.17, A o = —37cm ) corresponds to the A TT state of NCN [28,29]. The zero on the energy scale represents the minimum of the potential energy surface. Solid lines A = 0 vibronic levels dashed lines K = levels dash-dotted lines K = 1 levels dotted lines = 3 levels. Spin-vibronic levels are denoted by the value of the corresponding quantum number P P = Af - - E note that E is in this case spin quantum number),...

HCR and co-workers carried out a number of studies by employing 3D potential energy surfaces calculated by means of highly sophisticated ab initio approaches [88,91-101]. The results of these computations are in impressive agreement with the corresponding experimental findings. The discrepancies in the order of 100 wavenumbers, as in early ab initio studies [16,17], have been reduced in the HCR studies to only a few wavenumbers. In conclusion of their paper on the ( H ) system of NH2, Gabriel et al. state We believe... 

In Chapter VI, Ohm and Deumens present their electron nuclear dynamics (END) time-dependent, nonadiabatic, theoretical, and computational approach to the study of molecular processes. This approach stresses the analysis of such processes in terms of dynamical, time-evolving states rather than stationary molecular states. Thus, rovibrational and scattering states are reduced to less prominent roles as is the case in most modem wavepacket treatments of molecular reaction dynamics. Unlike most theoretical methods, END also relegates electronic stationary states, potential energy surfaces, adiabatic and diabatic descriptions, and nonadiabatic coupling terms to the background in favor of a dynamic, time-evolving description of all electrons. 

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. 

In the q = l limit, the effective temperature equals the standard temperature. Otherwise, adding a constant shift to the potential energy is equivalent to rescaling the temperature at which the canonical probability distribution is computed. 

For each pair of interacting atoms (/r is their reduced mass), three parameters are needed D, (depth of the potential energy minimum, k (force constant of the par-tictilar bond), and l(, (reference bond length). The Morse ftinction will correctly allow the bond to dissociate, but has the disadvantage that it is computationally very expensive. Moreover, force fields arc normally not parameterized to handle bond dissociation. To circumvent these disadvantages, the Morse function is replaced by a simple harmonic potential, which describes bond stretching by Hooke s law (Eq. (20)). 

---