Perturbation formulation

A careful analysis of the SCF energy decomposition for an inter-molecular interaction of A and B of various perturbation energy expressions indicate that for certain perturbation formulations there is a one-to-one correspondence between certain SCF energy decomposition terms and certain terms in the perturbation expressions. Thus one can calculate the values for the terms from energy decomposition of ab-initio or ab-initio MODPOT/VRDDO SCF wave functions and compare these to the values for the same type term resulting from the perturbation theory expressions. Care must be taken to correct for possible basis set incompleteness. 

K s> mentioned earlier, these formulations are applicable to structureless particles (bare ions). The only one of them that may be easily extended to complex ions is the perturbative formulation, either in the form of Bethe s theory for atomic targets or Lindhard s theory (dielectric formalism, DF) for the electron gas model. In addition, a comprehensive semiclassical approach, which extends the Bohr model to complex ions, has been developed more recently by Sigmund et al. [26]. 

In the following we will consider the energy loss of ions using two formulations the linear formulation based on the DF (i.e., a perturbative approach), and the non-linear (or non-perturbative) formulation, which will be extensively applied in the rest of the work. 

We shall refer to as the prompt-mode reactivity. The two expressions presented in Eq. (28) for the prompt-mode reactivity can be used—one for its calculation in the perturbation formulation, and the other for its experimental determination. 

The distribution-perturbation formulation of Eq. (156) also holds for inhomogeneous systems when the flux and flux perturbation are the solutions of, respectively, Eqs. (29) and (44). The corresponding adjoint function and its perturbation are the solutions of the equations... 

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

By employing the angle t defined by f68), the perturbative Hamiltonina H can be formulated in the form completely analogous to the Pople and Longuet-Higgins ansatz [69] ... 

Many physical properties of a molecule can be calculated as expectation values of a corresponding quantum mechanical operator. The evaluation of other properties can be formulated in terms of the "response" (i.e., derivative) of the electronic energy with respect to the application of an external field perturbation. 

Essentially all experimentally measured properties can be thought of as arising through the response of the system to some externally applied perturbation or disturbance. In turn, the calculation of such properties can be formulated in terms of the response of the energy E or wavefunction P to a perturbation. For example, molecular dipole moments p are measured, via electric-field deflection, in terms of the change in energy... 

Correlation can be added as a perturbation from the Hartree-Fock wave function. This is called Moller-Plesset perturbation theory. In mapping the HF wave function onto a perturbation theory formulation, HF becomes a hrst-order perturbation. Thus, a minimal amount of correlation is added by using the second-order MP2 method. Third-order (MP3) and fourth-order (MP4) calculations are also common. The accuracy of an MP4 calculation is roughly equivalent to the accuracy of a CISD calculation. MP5 and higher calculations are seldom done due to the high computational cost (A time complexity or worse). 

A second method is to use a perturbation theory expansion. This is formulated as a sum-over-states algorithm (SOS). This can be done for correlated wave functions and has only a modest CPU time requirement. The random-phase approximation is a time-dependent extension of this method. 

RPA, and CPHF. Time-dependent Hartree-Fock (TDFIF) is the Flartree-Fock approximation for the time-dependent Schrodinger equation. CPFIF stands for coupled perturbed Flartree-Fock. The random-phase approximation (RPA) is also an equivalent formulation. There have also been time-dependent MCSCF formulations using the time-dependent gauge invariant approach (TDGI) that is equivalent to multiconfiguration RPA. All of the time-dependent methods go to the static calculation results in the v = 0 limit. 

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. 

The derivative formulation is perhaps the easiest to understand. In this case the energy is expanded in a Taylor series in the perturbation strength A. 

Just as the variational condition for an HF wave function can be formulated either as a matrix equation or in terms of orbital rotations (Sections 3.5 and 3.6), the CPFIF may also be viewed as a rotation of the molecular orbitals. In the absence of a perturbation the molecular orbitals make the energy stationary, i.e. the derivatives of the energy with respect to a change in the MOs are zero. This is equivalent to the statement that the off-diagonal elements of the Fock matrix between the occupied and virtual MOs are zero. 

A lower max response at resonance was noted for poly butadiene-acrylic acid-containing pro-pints compared with polyurethane-containing opaque proplnts. Comparison of the measured response functions with predictions of theoretical models, which were modified to consider radiant-heat flux effects for translucent proplnts rather than pressure perturbations, suggest general agreement between theory and expt. The technique is suggested for study of the effects of proplnt-formulation variations on solid-proplnt combustion dynamics... 

Relations (2.46) and (2.47) are equivalent formulations of the fact that, in a dense medium, increase in frequency of collisions retards molecular reorientation. As this fact was established by Hubbard within Langevin phenomenology [30] it is compatible with any sort of molecule-neighbourhood interaction (binary or collective) that results in diffusion of angular momentum. In the gas phase it is related to weak collisions only. On the other hand, the perturbation theory derivation of the Hubbard relation shows that it is valid for dense media but only for collisions of arbitrary strength. Hence the Hubbard relation has a more general and universal character than that originally accredited to it. 

If, now, Si lies within the limits 1/t i < S2< /< i, which seem quite reasonable, the over-all effect of the perturbation is a transference of negative electricity from the various carbon atoms to the nitrogen atom at position 1, with the consequence, in accordance with the rule formulated above, that the molecule is deactivated so that substitution is more difficult than in benzene and furthermore the a and y positions (2 and 6, and 4, respectively) are most affected by this transference of electricity, so that substitution will take place at the 0 positions (3 and 5), which have the smallest deficiency of electrons. Both of these conclusions are borne out by experiment. 

Let us consider lithium as an example. In the usual treatment of this metal a set of molecular orbitals is formulated, each of which is a Bloch function built from the 2s orbitals of the atoms, or, in the more refined cell treatment, from 2s orbitals that are slightly perturbed to satisfy the boundary conditions for the cells. These molecular orbitals correspond to electron energies that constitute a Brillouin zone, and the normal state of the metal is that in which half of the orbitals, the more stable ones, are occupied by two electrons apiece, with opposed spins. 

---