Anharmonic term

Harmonic analysis (normal modes) at given temperature and curvature gives complete time behavior of the system in the harmonic limit [1, 2, 3]. Although the harmonic model may be incomplete because of the contribution of anharmonic terms to the potential energy, it is nevertheless of considerable importance because it serves as a first approximation for which the theory is highly developed. This model is also useful in SISM which uses harmonic analysis. 

It is clear that nonconfigurational factors are of great importance in the formation of solid and liquid metal solutions. Leaving aside the problem of magnetic contributions, the vibrational contributions are not understood in such a way that they may be embodied in a statistical treatment of metallic solutions. It would be helpful to have measurements both of ACP and A a. (where a is the thermal expansion coefficient) for the solution process as a function of temperature in order to have an idea of the relative importance of changes in the harmonic and the anharmonic terms in the potential energy of the lattice. 

The interatomic potential function for the diatomic molecule was described in Section 6 5. In the Taylpr-series development of this function (6-72)3 cubic and higher terms were neglected in the harmonic approximation. It is now of interest to evaluate the importance of these so-called anharmonic terms with the aid of the perturbation theory outlined above. If cubic and quartic... 

The operator (2.153) is the lowest-order approximation to the monopole term. Anharmonic terms can be introduced by considering the operator... 

However, most of the examples quoted in these earlier papers do not include the higher melting-point elements such as W, where a detailed treatment shows that the total entropy (at least of the solid phases) must include many other components such as the electronic specific heat, anharmonicity terms and the temperature dependence of 9d (Grimwall et al. 1987, Moroni et al. 1996). An estimate for the Debye temperature of the high-temperature 0 phase was included in the seminal... 

This result indicates that, according to the OPP model, the higher-order anharmonic terms have a stronger temperature dependence than the leading harmonic terms. The quartic terms, not specifically included in Eq. (2.43), have an even larger temperature dependence proportional to T3. Thus, the effect of... 

The vibrational displacements corresponding to the anharmonic terms in the potential are most pronounced in the directions away from the stronger bonding interactions, in which restoring forces are weaker. Thus, for the tetrahedral site symmetry of the diamond structure, the anharmonicity causes a larger mean-square displacement in directions opposite to the covalent bonds. At lower... 

More quantitatively, the effect of the thermal motion follows from the anharmonic thermal motion formalisms discussed in chapter 2. In the bcc structure, the relevant nonzero anharmonic term in the one-particle potential is the anisotropic, cubic site-symmetry allowed, part of uJuku um in expression (2.39). The modified potential for the cubic sites is given by (Willis 1969, Willis and Pryor 1975)... 

Since the anharmonic term is small relative to the leading harmonic term, the corresponding temperature factor can be written as... 

When the anharmonic terms are included in the full Hamiltonian, the resulting eigenstates, being mixtures of the ZO states, carry some brightness indicated by... 

Within the harmonic approximation the choice of a system of internal coordinates is irrelevant provided they are independent and that a complete potential function is considered ). For example, the vibrations of HjO can be analysed in terms of valence coordinates (r, >2, or interatomic coordinates (r, r, 3) and any difference in the accuracy to which observed energy levels are fitted (considering all the isotopic species H2O, HDO and D2O) will be due to the neglect of anharmonic terms. If one makes the approximation of a diagonal force field so that one is comparing the two potentials... 

To obtain the anharmonic terms in the potential, on the other hand, the choice of coordinates is important 130,131). The reason is that the anharmonic terms can only be obtained from a perturbation expansion on the harmonic results, and the convergence of this expansion differs considerably from one set of coordinates to another. In addition it is usually necessary to assume that some of the anharmonic interaction terms are zero and this is true only for certain classes of internal coordinates. For example, one can define an angle bend in HjO either by a rectilinear displacement of the hydrogen atoms or by a curvilinear displacement. At the harmonic level there is no difference between the two, but one can see that a rectilinear displacement introduces some stretching of the OH bonds whereas the curvilinear displacement does not. The curvilinear coordinate follows more closely the bottom of the potential well (Fig. 12) than the linear displacement and this manifests itself in rather small cubic stretch-bend interaction constants whereas these constants are larger for rectilinear coordinates. A final and important point about the choice of curvilinear coordinates is that they are geometrically defined (i.e. independent of nuclear masses) so that the resulting force constants do not depend on isotopic species. At the anharmonic level this is not true for rectilinear coordinates as it has been shown that the imposition of the Eckart conditions, that the internal coordinates shall introduce no overall translation or rotation of the body, forces them to have a small isotopic dependence 132). 

Force fields up to quartic anharmonic terms are now known with reasonably high accuracy for several triatomic molecules and the results shown in Table 3 for H2O are typical. However, even for these there has had to be an assumption that some of the quartic interaction terms are zero in order that the equations from which the constants are derived shall have unique solutions. It can be seen moreover that some of the cubic and quartic terms have uncertainties which are larger than the values of the constants themselves. 

Figure 4 clearly illustrates that polarizability is a function of the frequency of the applied field. Changing the restoring force constant, k (equation (2)) is another way to modify the linear polarizability. Another alternative is to add anharmonic terms to the potential to obtain a surface such as that shown in Figure 13. The restoring force on the electron is no longer linearly proportional to its displacement during the polarization by the light wave, it is now nonlinear (Figure 14). As a first approximation (in one dimension) the restoring force could be written as ...

Figure 13. Plot of potential energy vs. distortion coordinate for a material with a harmonic potential and a material with an additional cubic anharmonic term.

The structure/property relationships that govern third-order NLO polarization are not well understood. Like second-order effects, third-order effects seem to scale with the linear polarizability. As a result, most research to date has been on highly polarizable molecules and materials such as polyacetylene, polythiophene and various semiconductors. To optimize third- order NLO response, a quartic, anharmonic term must be introduced into the electronic potential of the material. However, an understanding of the relationship between chemical structure and quartic anharmonicity must also be developed. Tutorials by P. Prasad and D. Eaton discuss some of the issues relating to third-order NLO materials. 

Since the ab initio calculation on coronene monoanion indicates that the stable configuration has C2h symmetry, we must take higher-order anharmonic terms up to sixth order into consideration which is not considered in the derivation of the epikernel principle to obtain a JT surface with the adequate structure. We will discuss this point later. 

This result implies that the minimum associated with the 3z2 — r2 electronic state is the lowest one if the anharmonic term dominates and A3 > 0 or if the quadratic coupling is the leading term and Vq is positive. The curvature at the minimum of E(3z2 r2 Qe) and E(x2 — y2 Qe) curves also depends on both the anharmonic and quadratic coupling terms... 

---