The Analytic Derivative Method

Let us consider again a system in the presence of a general perturbation Pa that could be a component of the electric field a, field gradient a, magnetic induction Ba or nuclear moment with the first- and second-order perturbation Hamiltonians as given in Eq. (2.101) and Eqs. (2.108)-(2.110). 

The first derivative of the energy of such a perturbed system with respect to Pa ean then be written for most methods as 

However, for non-variational wavefunctions, as for example in the case of Moller-Plesset perturbation theory and the coupled cluster methods, the density matrices here are not consistent with the definition in Eq. (2.20) and thus Eqs. (9.107) and (9.108). The density matrices defined in Eq. (12.1) are for those methods therefore 

The relaxed density matrix can in general be decomposed in an SCF and correlation 

The SCF density is given in Eq. (9.112) and the correlation part consists of two parts 

---

For k = 0, the Fock matrix and its derivatives with respect to the displacements of the nuclei are always block diagonal. Then one can directly apply the analytical derivative methods developed for finite systems to extended systems [69,86,87,88]. But when the displacements break the translational symmetry, the Fock matrix and its derivatives are no longer block diagonal. To solve the CPHF equations, one needs to use the symmetrized (normal mode) coordinates instead of the Cartesian coordinates of the nuclei. Efficient analytical methods have been developed to calculate the energy derivatives for k / 0 with both plane wave [89-90] and general basis functions [85]. The latter can be functions of nuclear coordinates and have linear dependence. These methods reduce the computational cost required to calculate the phonon spectrum with k 7 0 to the same as that needed for the spectrum at k = 0. 

This is obtained directly in ab initio calculations, particularly those based upon the analytic derivative methods discussed later. Another advantage of working with the polar tensors is that the relationship between the Cartesian and the normal coordinate derivatives is straightforward... 

The various response tensors are identified as terms in these series and are calculated using numerical derivatives of the energy. This method is easily implemented at any level of theory. Analytic derivative methods have been implemented using self-consistent-field (SCF) methods for a, ft and y, using multiconfiguration SCF (MCSCF) methods for ft and using second-order perturbation theory (MP2) for y". The response properties can also be determined in terms of sum-over-states formulation, which is derived from a perturbation theory treatment of the field operator — [iE, which in the static limit is equivalent to the results obtained by SCF finite field or analytic derivative methods. 

Density Functional Theory does not require specific modifications, in relation to the solvation terms [9], with respect to the Hartree-Fock formalism presented in the previous section. DFT also absorbs all the properties of the HF approach concerning the analytical derivatives of the free energy functional (see also the contribution by Cossi and Rega), and as a matter of fact continuum solvation methods coupled to DFT are becoming the routine approach for studies of solvated systems. 

In order to correctly design analytical procedures used for the detection of food allergens, it is necessary to have basic knowledge of food product chemistry to know how to collect, prepare, and store food samples to be able to fragment, mix, disintegrate, and extract samples to know (or be able to find quickly) relevant food quality standards and admissible contents of particular food ingredients and finally to understand precision of determinations, their sensitivity, and detection threshold levels, reproducibility, and errors of determination methods. In addition, it is essential to be able to gather the results of assays, process them with the aid of a computer and statistical methods, and to present the analytically derived data. 

Probably even more important to computational quantum chemistry is the development in the analytical evaluation of first, second and higher derivatives of the potential energy with respect to nuclear coordinates . These analytical derivative methods are indispensable to the location and characterization of the stationary points (minima or transition states) on the potential energy surface and have greatly advanced the scope of applicability of ab initio calculations. Ab initio calculations are in a position to predict many new types of the heavier group 15 compounds and provide valuable information for the interpretation of complex experimental data. 

The second-derivative method is an extension of the first-derivative method. The second-derivative of the data changes sign at the point of inflection in the titration curve. This change is often used as the analytical signal in automatic titrators. 

The simplest esterification reagents are the corresponding alcohols ROH themselves. Different esters have been used as the analytical derivatives of carboxylic adds Me, Et, Pr, iso-Pr, isomeric Bu (excluding tert-B i esters, owing to their poorer synthetic yields), and so forth. This method requires the use of excess of dry alcohol and acid catalysis by BCI3, BF3, CH3COCI, SOCI2,... 

In indirect methods, heat content measurements are performed over a large temperature range, for instance by drop calorimetry, and Cp is derived by the analytic derivation of heat content plots versus temperature. 

Also the algebraic diagrammatic construction (ADC) method that was discussed in Section VI has been applied to the polarization propagator (Schirmer, 1982). Diagrammatic rules, rather than the analytic derivation used in SOPPA, are applied to formulate the second-order ADC(2) and basically the same approximation is obtained. 

Since their introduction in the late 1960s, gradient methods, or more properly analytical derivative methods, have become one of the most vigorously developing topics on modern quantum chemistry. They have also acquired considerable significance for the solution of practical chemical problems. The first review on this subject was written in 1974-75, although it was published much later (Pulay, 1977) a short chapter, limited to first derivatives,... 

The non-relativistic PolMe (9) and quasirelativistic NpPolMe (10) basis sets were used in calculations reported in this paper. The size of the [uncontractd/contracted] sets for B, Cu, Ag, and Au is [10.6.4./5.3.2], [16.12.6.4/9.7.3.2], [19.15.9.4/11.9.5.2], and [21.17.11.9/13.11.7.4], respectively. The PolMe basis sets were systematically generated for use in non-relativistic SCF and correlated calculations of electric properties (10, 21). They also proved to be successful in calculations of IP s and EA s (8, 22). Nonrelativistic PolMe basis sets can be used in quasirelativistic calculations in which the Mass-Velocity and Darwin (MVD) terms are considered (23). This follows from the fact that in the MVD approximation one uses the approximate relativistic hamiltonian as an external perturbation with the nonrelativistic wave function as a reference. At the SCF and CASSCF levels one can obtain the MVD quasi-relativistic correction as an expectation value of the MVD operator. In perturbative CASPT2 and CC methods one needs to use the MVD operator as an external perturbation either within the finite field approach or by the analytical derivative schems. The first approach leads to certain numerical accuracy problems. 

NACMEs based on state-averaged MCSCF wave functions and analytic derivative methods. This should provide NACME s of similar overall accuracy to that obtained for the adiabatic potentials. 

The importance of analytic derivative methods in quantum chemistry cannot be overstated. Analytic methods have been demonstrated to be more efficient than are corresponding finite difference techniques. Calculation of the first derivatives of the energy with respect to the nuclear coordinates is perhaps the most common these provide the forces on the nuclei and facilitate the location of stationary points on the potential energy hypersurface. Differentiating the electronic energy with respect to a parameter x (which may be, but is not required to be, a nuclear coordinate), leads to the well-known expression... 

B.R. Brown The organic chemistry of aliphatic nitrogen compounds Y. Yamaguchi, Y. Osamura, J.D. Goddard, and H.F. Schaeffer A new dimension to quantum chemistry analytic derivative methods in ab initio molecular electronic structure theory... 

The farsighted development of methods to calculate the analytical derivatives of the ab initio energy with respect to the Cartesian coordinates of the nuclei by Pulay [36] played a key role in the practical implementation of TST. Further developments by Handy and Schaefer [37a] Schlegel et al. [37b], and Johnson and Frisch [38], among others greatly expanded the number of ab initio methods for which these derivatives were available. The tools to evaluate the partition functions were at hand, but first we must locate the TS. 

Some of the above reasons for preferring the energy derivative over the expectation value will only hold for variational wavefunctions. However, Diercksen et have argued that this technique is more suitable with MBPT methods as well. The techniques developed in analytic derivative methods can also be applied to the calculation of MBPT properties. In Moller-Plesset theory (the simplest form of MBPT), the zeroth-order wavefunction is SCF and the Hamiltonian is partitioned so that... 

---