Oscillation function functions

An alternative to using a superposition of Gaussian functions is to extend the basis set by using Hermite polynomials, that is, hamonic oscillator functions [24]. This provides an orthonormal, in principle complete, basis set along the bajectoiy, and the idea has been taken up by Billing [151,152]. The basic problem with this approach is the slow convergence of the basis set. 

Xvf (Ra - Ra,e) Xvi> will be non-zero and probably quite substantial (because, for harmonic oscillator functions these "fundamental" transition integrals are dominant- see earlier) ... 

Let us consider an example, that of the derivative operator in the orthonormal basis of Harmonic Oscillator functions. The fact that the solutions of the quantum Harmonic... 

Figure 3.14 General representations of stress and strain out of phase by amount 5 (a) represented by oscillating functions and (b) represented by vectors.

Owing to the effects of mechanical anharmonicity - to which we shall refer in future simply as anharmonicity since we encounter electrical anharmonicity much less frequently -the vibrational wave functions are also modified compared wifh fhose of a harmonic oscillator. Figure 6.6 shows some wave functions and probabilify densify functions (iA A ) for an anharmonic oscillator. The asymmefry in and (iA A ) 5 compared wifh fhe harmonic oscillator wave functions in Figure f.i3, increases fheir magnitude on the shallow side of the potential curve compared with the steep side. 

We proceed now to the calculation of B, following [Benderskii et al. 1992a]. The denominator in (4.11) (apart from normalization) is equal to the harmonic-oscillator partition function [2 sinh(ico+ )] The numerator is the product of the s satisfying an equation of the Shrodinger type... 

The potential energy is often described in terms of an oscillating function like the one shown in Figure 10.9(a) where the minima correspond to the relative orientations in which the interactions are most favorable, and the maxima correspond to unfavorable orientations. In ethane, the minima would occur at the staggered conformation and the maxima at the eclipsed conformation. In symmetrical molecules like ethane, the potential function reflects the symmetry and has a number of equivalent maxima and minima. In less symmetric molecules, the function may be more complex and show a number of minima of various depths and maxima of various heights. For our purposes, we will consider only molecules with symmetric potential functions and designate the number of minima in a complete rotation as r. For molecules like ethane and H3C-CCI3, r = 3. 

It Is now well established experimentally that the solvation force, fg, of confined fiuld Is an oscillating function of pore wall separation. In Figure 4 we compare the theoretical and MD results for fg as a function of h. Given that pressure predictions are very demanding of a molecular theory, the observed agreement between our simple theory and the MD simulations must be viewed as quite good. The local maxima and minima In fg coincide with those In n y and therefore also refiect porewldths favorable and unfavorable to an Integral number of fiuld layers. 

At x = 0, B(0) is equal to the uniform density of electrons. The first term of the right hand side makes a bulk peak around x = 0. It sharply damps outside, because the k-integration over the occupied states is similar in structure to the following damping oscillation function ... 

Some of the Hermite polynomials and the corresponding harmonic-oscillator wave functions are presented in Thble 1. The importance of the parity of these functions under the inversion operation, cannot be overemphasized. 

To solve this equation, an appropriate basis set ( >.,( / ) is required for the nuclear functions. These could be a set of harmonic oscillator functions if the motion to be described takes place in a potential well. For general problems, a discrete variable representation (DVR) [100,101] is more suited. These functions have mathematical properties that allow both the kinetic and potential energy... 

The eigenvalues of this Hamitonian can calculated by numerical diagonalization of the truncated matrix of the quantum system in the basis of the harmonic oscillator wave functions. The matrix elements of Hq and V are... 

In this approximation the nuclear wavefunctions are a product of N harmonic oscillator functions, one for each normal mode ... 

In the harmonic approximation the functions Xi and Xf are products of harmonic oscillator functions. We therefore specify the initial state by a set of quantum numbers n — (ni, ri2,..., n/v), and those for the final state by m = (mi,m2,..., tun)- So the nuclear wavefunctions are henceforth denoted by Xi,n and Xf,m- Equation (19.21) tells us how to calculate the rate of transition from one particular initial quantum mode n to a final quantum state m. This is more than we want to know. All we are interested in is the total rate from any initial state to any final state. The ensemble of reactants is in thermal equilibrium therefore... 

Now it becomes apparent why it was useful to replace the delta function by its Fourier transform. The wavefUnctions Xin are products of harmonic oscillator functions, the Hamiltonians Hi and H/ are sums of harmonic oscillator terms. Therefore the terms in the brackets factorize in the form ... 

Fig. 4.1 The zero point energy or low temperature approximation As temperature drops and u increases above u 4 the harmonic oscillator partition function Q (Harm. Osc.) is better and better approximated by the zero point energy term, exp(—u/2). For a typical CH stretching frequency, v = 3000 cm-1, u 4 at 1050 K and it is reasonable to use the ZPE approximation for that frequency at temperatures below 1000 k...

Using harmonic oscillator partition functions to describe both internal and external modes, the logarithmic Q ratios introduced above, ln(Qg/Qc/QgQc/) = ln(Qc/QcO + ln(Qg7Qg), become... 

Ascribe to q the form of a harmonic-oscillator partition function so that the preexponential factor becomes ... 

For the special case of the self-correlation function (n=m) B n,n,t) reveals the mean-square displacement of a polymer segment. For large p the cos in Eq. 3.16 is a rapidly oscillating function which may be replaced by the mean-value 1/2. With this approximation we can convert the sum into an integral and obtain ... 

In another example of differences in complexity, the bondstretching energy in CHARMM is calculated with a harmonic oscillator function. MM3 solves the problem described by French, Tran and Perez in this book for MM2 s cubic stretching function by using a quartic function for bond stretching. Additional complexity in MM3 is described in Ref. 12. 

---