Perturbative methods

Let us assume that stress gradient in axial direction is present but smooth. Then we can use a perturbation method and expand the solution of equation (30) in a series. The first term of this expansion will be a solution of the plane strain problem and potential N will be equal to zero. The next terms of the stress components will contain potential N also. 

In integrated photoelasticity it is impossible to achieve a complete reconstruction of stresses in samples by only illuminating a system of parallel planes and using equilibrium equations of the elasticity theory. Theory of the fictitious temperature field allows one to formulate a boundary-value problem which permits to determine all components of the stress tensor field in some cases. If the stress gradient in the axial direction is smooth enough, then perturbation method can be used for the solution of the inverse problem. As an example, distribution of stresses in a bow tie type fiber preforms is shown in Fig. 2 [2]. 

Claverie P 1971 Theory of intermolecular forces. I. On the inadequacy of the usual Rayleigh-Schrddinger perturbation method for the treatment of intermolecular forces Int. J. Quantum Chem. 5 273... 

Zwanzig R 1954 High temperature equation of state by a perturbation method I. Nonpolar Gases J. Chem. Phys. 22 1420... 

Ramsden J J 1998 Towards zero-perturbation methods for investigating biomolecular interactions Coiioids Surfaces A 141 287-94... 

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. 

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. 

Bash, P.A., Field, M.J.,Karplus, M. Free energy perturbation method for chemical reactions in the condensed phase A dynamical approach baaed on a combined quantum and molecular dynamics potential. J. Am. Chem. Soc. 109 (1987) 8092-8094. 

This principle has been applied in a contribution by Mark, Schafer, Liu and van Gunsteren to this volume, and in section 6 of this article. For a review of free energy perturbation methods see [8]. 

The problems that occur when one tries to estimate affinity in terms of component terms do not arise when perturbation methods are used with simulations in order to compute potentials of mean force or free energies for molecular transformations simulations use a simple physical force field and thereby implicitly include all component terms discussed earlier. We have used the molecular transformation approach to compute binding affinities from these first principles [14]. The basic approach had been introduced in early work, in which we studied the affinity of xenon for myoglobin [11]. The procedure was to gradually decrease the interactions between xenon atom and protein, and compute the free energy change by standard perturbation methods, cf. (10). An (issential component is to impose a restraint on the... 

T.P. Lybrand, Computer simulations of biomolecular systems using molecular dynamics and free energy perturbation methods, in Reviews in Computational Chemistry, Vol. 1, K.B. Lipkowitz, D.B. Boyd (Eds.), VCH, New York, 1990, pp. 295-320. 

This perturbation method is claimed to be more efficient than the fluctuating dipole method, at least for certain water models [Alper and Levy 1989], but it is important to ensure that the polarisation (P) is linear in the electric field strength to avoid problems with dielectric saturation. 

An early application of the free energy perturbation method was the determination of t] tree energy required to create a cavity in a solvent. Postma, Berendsen and Haak determin the free energy to create a cavity (A = 1) in pure water (A = 0) using isothermal-isobai... 

Calculations of relative partition coefficients have been reported using the free energy perturbation method with the molecular dynamics and Monte Carlo simulation methods. For example, Essex, Reynolds and Richards calculated the difference in partition coefficients of methanol and ethanol partitioned between water and carbon tetrachloride with molecular dynamics sampling [Essex et al. 1989]. The results agreed remarkably well with experiment... 

Energy Perturbation Methods. In Lipkowitz K B and D B Boyd (Editors) Reviews in Compiitalio Oieniistry Volume 1. New York, VCH Publishers, pp. 295-320. 

Zwanzig R W 1954. High-temperature Equation of State by a Perturbation Method. 1. Nonpolar Gases. Journal of Chemical Physics 22 1420-1426. 

The MoIIer-PIesset perturbation method (MPPT) uses the single-eonfiguration SCF proeess (usually the UHF implementation) to first determine a set of LCAO-MO eoeffieients and, henee, a set of orbitals that obey F( )i = 8i (jii. Then, using an unperturbed Hamiltonian equal to the sum of these Foek operators for eaeh of the N eleetrons =... 

The various studies attempting to increase our understanding of turbulent flows comprise five classes moment methods disregarding probabiUty density functions, approximation of probabiUty density functions using moments, calculation of evolution of probabiUty density functions, perturbation methods beginning with known stmctures, and methods identifying coherent stmctures. For a thorough review of turbulent diffusion flames see References 41—48. 

Perturbation Methods If the ordinary differential equation has a parameter that is small and is not multiplying the highest derivative, perturbation methods can give solutions for small values of the parameter. 

Various techniques exist that make possible a normal mode analysis of all but the largest molecules. These techniques include methods that are based on perturbation methods, reduced basis representations, and the application of group theory for symmetrical oligomeric molecular assemblies. Approximate methods that can reduce the computational load by an order of magnitude also hold the promise of producing reasonable approximations to the methods using conventional force fields. 

To overcome the limitations of the database search methods, conformational search methods were developed [95,96,109]. There are many such methods, exploiting different protein representations, objective function tenns, and optimization or enumeration algorithms. The search algorithms include the minimum perturbation method [97], molecular dynamics simulations [92,110,111], genetic algorithms [112], Monte Carlo and simulated annealing [113,114], multiple copy simultaneous search [115-117], self-consistent field optimization [118], and an enumeration based on the graph theory [119]. 

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). 

The main limitation of perturbation methods is the assumption that the zero-order wave function is a reasonable approximation to the real wave function, i.e. the perturbation operator is sufficiently small . The poorer the HF wave function describes... 

Just as single reference Cl can be extended to MRCI, it is also possible to use perturbation methods with a multi-detenninant reference wave function. Formulating MR-MBPT methods, however, is not straightforward. The main problem here is similar to that of ROMP methods, the choice of the unperturbed Hamilton operator. Several different choices are possible, which will give different answers when the tlieory is carried out only to low order. Nevertheless, there are now several different implementations of MP2 type expansions based on a CASSCF reference, denoted CASMP2 or CASPT2. Experience of their performance is still somewhat limited. 

Perturbation methods add all types of corrections (S, D, T, Q etc.) to the reference wave function to a given order (2, 3, 4 etc.). The idea in Coupled Cluster (CC) methods is to include all corrections of a given type to infinite order. The (intermediate normalized) coupled cluster wave function is written as... 

---