Perturbation theory solution

A FIGURE 4.8 Radial distribution functions for the perturbation theory solution to helium. The average radial distribution functions for each of the two electrons in the zero-order and first-order wavefunctions are shown. The zero-order case is the same distribution as occupied by the single electron in Is He, because we have turned off the electron-electron repulsion. Once we turn on the repulsion (by adding the first-order correction to the energy), each electron pushes the other away. As a result, the electron density expands to a larger average distance from the nucleus. 

The determination of the good actions describing vibration-rotation motion requires the solution of the molecular Hamilton-Jacobi equation, which is a nonlinear partial differential equation in 3Na"5 variables (including rotation), where is the number of atoms. Even for = 3 (a triatomic molecule) an exact solution to this equation is extremely complex computationally, and it is not practical for collisional applications. Several approximations can be used to simplify this treatment, however, including (i) the separation of vibration from rotation (valid in the limit of an adequate vibration-rotation time scale separation), and (ii) the use of classical perturbation theory (in 2nd and 3rd order) to solve the three-dimensional vibrational Hamilton-Jacobi equation which remains after the separation of rotation. Details of both the separation procedures and the perturbation-theory solution are discussed elsewhere. For the present application, the validity of the first... 

Another subject with important potential application is discussed in Section XIV. There we suggested employing the curl equations (which any Bohr-Oppenheimer-Huang system has to obey for the for the relevant sub-Hilbert space), instead of ab initio calculations, to derive the non-adiabatic coupling terms [113,114]. Whereas these equations yield an analytic solution for any two-state system (the abelian case) they become much more elaborate due to the nonlinear terms that are unavoidable for any realistic system that contains more than two states (the non-abelian case). The solution of these equations is subject to boundary conditions that can be supplied either by ab initio calculations or perturbation theory. 

In applying quantum mechanics to real chemical problems, one is usually faced with a Schrodinger differential equation for which, to date, no one has found an analytical solution. This is equally true for electronic and nuclear-motion problems. It has therefore proven essential to develop and efficiently implement mathematical methods which can provide approximate solutions to such eigenvalue equations. Two methods are widely used in this context- the variational method and perturbation theory. These tools, whose use permeates virtually all areas of theoretical chemistry, are briefly outlined here, and the details of perturbation theory are amplified in Appendix D. 

Perturbation theory is then used to express the eorreetions to these zeroth order solutions as indieated in Appendix D. 

This Introductory Section was intended to provide the reader with an overview of the structure of quantum mechanics and to illustrate its application to several exactly solvable model problems. The model problems analyzed play especially important roles in chemistry because they form the basis upon which more sophisticated descriptions of the electronic structure and rotational-vibrational motions of molecules are built. The variational method and perturbation theory constitute the tools needed to make use of solutions of... 

There is another manner in whieh perturbation theory is used in quantum ehemistry that does not involve an externally applied perturbation. Quite often one is faeed with solving a Sehrodinger equation to whieh no exaet solution has been (yet) or ean be found. In sueh eases, one often develops a model Sehrodinger equation whieh in some sense is designed to represent the system whose full Sehrodinger equation ean not be solved. The... 

The fugacity coefficient of thesolid solute dissolved in the fluid phase (0 ) has been obtained using cubic equations of state (52) and statistical mechanical perturbation theory (53). The enhancement factor, E, shown as the quantity ia brackets ia equation 2, is defined as the real solubiUty divided by the solubihty ia an ideal gas. The solubiUty ia an ideal gas is simply the vapor pressure of the sohd over the pressure. Enhancement factors of 10 are common for supercritical systems. Notable exceptions such as the squalane—carbon dioxide system may have enhancement factors greater than 10. Solubihty data can be reduced to a simple form by plotting the logarithm of the enhancement factor vs density, resulting ia a fairly linear relationship (52). 

Free energy calculations rely on the following thermodynamic perturbation theory [6-8]. Consider a system A described by the energy function = 17 + T. 17 = 17 (r ) is the potential energy, which depends on the coordinates = (Fi, r, , r ), and T is the kinetic energy, which (in a Cartesian coordinate system) depends on the velocities v. For concreteness, the system could be made up of a biomolecule in solution. We limit ourselves (mostly) to a classical mechanical description for simplicity and reasons of space. In the canonical thermodynamic ensemble (constant N, volume V, temperature T), the classical partition function Z is proportional to the configurational integral Q, which in a Cartesian coordinate system is... 

Back now to the simpler case of a single perturbation. Perturbation theory aims to write solutions for... 

A Perturbation Theory is developed for treating a system of n electrons in which the Hartree-Fock solution appears as the zero-order approximation. It is shown by this development that the first order correction for the energy and the charge density of the system is zero. The expression for the second order correction for the energy greatly simplifies because of the special property of the zero order solution. It is pointed out that the development of the higher order approximation involves only calculations based on a definite one-body problem. 

If we add a perturbation A then the self-consistency is destroyed and we need to re-do the iterative HF-LCAO calculation. The idea of self-consistent perturbation theory is to seek solutions of the perturbed HF-LCAO equations... 

The idea in perturbation methods is that the problem at hand only differs slightly from a problem which has already been solved (exactly or approximately). The solution to the given problem should therefore in some sense be close to the solution of the already known system. This is described mathematically by defining a Hamilton operator which consists of two part, a reference (Hq) and a perturbation (H )- The premise of perturbation methods is that the H operator in some sense is small compared to Hq. In quantum mechanics, perturbational methods can be used for adding corrections to solutions which employ an independent particle approximation, and the theoretical framework is then called Many-Body Perturbation Theory (MBPT). 

Up to this point we are still dealing with undetermined quantities, energy and wave funetion corrections at each order. The first-order equation is one equation with two unknowns. Since the solutions to the unperturbed Schrddinger equation generates a complete set of functions, the unknown first-order correction to the wave function can be expanded in these functions. This is known as Rayleigh-Schrddinger perturbation theory, and the equation in (4.32) becomes... 

In developing perturbation theory it was assumed that the solutions to the unpermrbed problem formed a complete set. This is general means that there must be an infinite number of functions, which is impossible in actual calculations. The lowest energy solution to the unperturbed problem is the HF wave function, additional higher energy solutions are excited Slater determinants, analogously to the Cl method. When a finite basis set is employed it is only possible to generate a finite number of excited determinants. The expansion of the many-electron wave function is therefore truncated. 

This is an eigenvalue problem of the form of Eq. III.45 referring to the truncated basis only, and the influence of the remainder set is seen by the additional term in the energy matrix. The relation III.48 corresponds to a solution of the secular equation by means of a modified perturbation theory,19 and the problem is complicated by the fact that the extra term in Eq. III.48 contains the energy parameter E, which leads to an iteration procedure. So far no one has investigated the remainder problem in detail, but Eq. III.48 certainly provides a good starting point. 

For one-dimensional rotation (r = 1), orientational correlation functions were rigorously calculated in the impact theory for both strong and weak collisions [98, 99]. It turns out in the case of weak collisions that the exact solution, which holds for any happens to coincide with what is obtained in Eq. (2.50). Consequently, the accuracy of the perturbation theory is characterized by the difference between Eq. (2.49) and Eq. (2.50), at least in this particular case. The degree of agreement between approximate and exact solutions is readily determined by representing them as a time expansion... 

A simple repetition of the iteration procedure (2.20)-(2.22) results in divergence of higher order solutions. However, a perturbation theory series may be summed up so that all unbound diagrams are taken into account, just as is usually done for derivation of the Dyson equation [120]. As a result P satisfies the integral-differential equation... 

In the DC-biased structures considered here, the dynamics are dominated by electronic states in the conduction band [1]. A simplified version of the theory assumes that the excitation occurs only at zone center. This reduces the problem to an n-level system (where n is approximately equal to the number of wells in the structure), which can be solved using conventional first-order perturbation theory and wave-packet methods. A more advanced version of the theory includes all of the hole states and electron states subsumed by the bandwidth of the excitation laser, as well as the perpendicular k states. In this case, a density-matrix picture must be used, which requires a solution of the time-dependent Liouville equation. Substituting the Hamiltonian into the Liouville equation leads to a modified version of the optical Bloch equations [13,15]. These equations can be solved readily, if the k states are not coupled (i.e., in the absence of Coulomb interactions). 

Since the spin-orbit interaction energy is small, the solution of equations (7.43) to obtain E is most easily accomplished by means of perturbation theory, a technique which is presented in Chapter 9. The evaluation of E is left as a problem at the end of Chapter 9. 

Perturbation theory provides a procedure for finding approximate solutions to the Schrodinger equation for a system which differs only slightly from a system for which the solutions are known. The Hamiltonian operator H for the system of interest is given by... 

---