Fourth-order, generally approximations

A method similar to the iterative, is the partial closure method [37], It was formulated originally as an approximated extrapolation of the iterative method at infinite number of iterations. A subsequent more general formulation has shown that it is equivalent to use a truncated Taylor expansion with respect to the nondiagonal part of T instead of T-1 in the inversion method. An interpolation of two sets of charges obtained at two consecutive levels of truncations (e.g. to the third and fourth order) accelerates the convergence rate of the power series [38], This method is no longer in use, because it has shown serious numerical problems with CPCM and IEFPCM. 

All calculations in Ref. [22] were performed utilizing the Gaussian-98 code [30]. The potential energy scan was performed by means of the Mqller-Plesset perturbation theory up to the fourth order (MP4) in the frozen core approximation. The electronic density distribution was studied within the population analysis scheme based on the natural bond orbitals [31,32], A population analysis was performed for the SCF density and MP4(SDQ) generalized density determined applying the Z-vector concept [33]. 

Published kinetic data were generally obtained in batch reactors (1-5). The data obtained were observed to fit first order kinetics notwithstanding the complexity of the feeds studied and the constant change in the nature of the product-forming intermediates. It has been shown (3,5) that batch coking has a definite induction period and that the usually observed first order coking behavior of complex feeds is only apparent. Also it was determined that the rates observed could fit third order kinetics for decomposition and fourth order for polymerization/ condensation (5). Clearly, kinetics of batch coking are only approximations. 

Table 2 explores the aforementioned notion that MP2 serves as a good approximation to full MP4 calculations of H-bond energies. The basis set notation is taken directly from ref. 42, where M refers to a minimal basis set, D to double-C, and Q to quadruple- . The numbers of polarization functions are displayed in the following parentheses. We first focus on the data on the left side of the table, which has been corrected for BSSE by the counterpoise procedure. The full fourth order correction A (4) (AT-(- A -l- A is generally on the order of -0.5 kcal/mol. Note that the minimal type of basis set, even with... 

We see that each additional teim in the Taylor series requires storage of another set of N delayed variables. Our experience suggests that a fourth-order expansion generally gives satisfactory results and that increasing the number of steps per interval N is more efficacious with stiff DDE systems than is increasing the order of the Taylor series approximation. 

In the above discussion we assumed that the pde is in the form of (3.26) and that the spatial grid is uniform. For more general pdes involving first spatial derivatives as well, Bieniasz [15] developed a three-point, fourth order accurate compact boundary gradient approximation. In the case of nonuniform spatial grids the derivation of compact gradient approximations becomes more complicated (see, for example, Bieniasz [16]) and relevant boundary gradient formulae are not yet available. 

General experience, however, is that fourth-order MBPT overestimates the triples correction. A somewhat more reliable approach is obtained when the second-order MBPT amplitudes in equation (44) are replaced by the converged CCSD amplitudes. This leads to the so-called CCSD-I-T(CCSD) approximation. Its performance is reasonable, but it still tends to exaggerate the triples contribution in more difficult cases. ... 

If /I > 4, there is no formula which gives the roots of the general equation. For fourth and higher order (even third order), the roots can be found numerically (see Numerical Analysis and Approximate Methods ). However, there are some general theorems that may prove useful. 

---