Theoretical Approach to Potential Energy Functions for Linear AB2 and ABC and Bent AB2 Triatomic Molecules

Abstract

Let the coordinate system for a linear AB2 triatomic molecule have its origin on the A nucleus with R1 the distance to one B nucleus, R2 the distance to the other, and θ the apex angle. Then a working formula for the Born-Oppenheimer potential energy near equilibrium, W(R1, R2, θ), is W(R1, R2, θ) =WD(R1)+WD(R2)+A/|R3|N−B/(R1+R2)N,where W D(R1) and W D(R2) are potential functions for the ground state diatomic molecule AB, R3 is the vector sum R1+R2, A and B are constants and N is an integer. The potential energy for linear ABC tri-atomic molecules is given by Eq. (i) with R2 [or R1] scaled: R2→ η R2. Equation (i) is tested for CO2, CS2, OCS, HCN, and N2O by predicting all force constants up to fourth order except for the harmonic bending constant which is used in the parameterization. For bent AB2 molecules the working formula for the potential energy is W(R1, R2, θ) = WD(R1)+WD(R2)+A/|R3|N−B/(R1+R2)N+C/|R3|M−Aθ cos(θ − θe), (M≠ N),where θe is the equilibrium angle. Equation (ii) is tested for H2O, SO2, and O3 by predicting all force constants up to fourth order except for the harmonic bending and stretch-bend interaction constants which are used in the parameterization. A theory is presented for the presence of the W D(R) components in Eqs. (i) and (ii). It is suggested that the inverse |R3| terms partially represent averaged multipole interaction energies while the inverse (R1+R2) terms and angular term in Eq. (ii) partially represent valence orbital effects. The formulas for the parameters in Eqs. (i) and (ii) are given. Suggestions for writing down potential functions for larger molecules are included.