Physical approximations

In MOT the goal is compute the electronic wavefnnction by solving the Schrodinger equation (i.e., non-relativistic quantum mechanics). Many textbooks are available on the topic. The classic text by Hehre et al. [61] emphasizes understanding and planning actual calculations. This is complemented by the theoretical approach taken by Szabo and Ostlund [62]. Unlike SEMOT, no empirical parameters are used. However, approximations are necessary to make the problem computationally tractable. The approximations can be classihed as either physical or numerical. The most serious physical approximations involve the treatment of the electrostatic repulsion among the electrons. The principal numerical approximation is the choice of mathematical functions for describing the molecular orbitals. Popular approximations are described in this section. 

Nearly all calculations ignore relativistic effects, as mentioned earlier, and accept the Born-Oppenheimer approximation, which is that the nuclei are so much heavier than the electrons that they may be regarded as inhnitely more massive. This separates the problems of determining the electronic and nuclear wavefunctions. Eor the electronic problem, the simplest physical approximation is embodied in Hartree-Eock (HE) theory. This is a mean-held approximation, in which each electron moves independently of the other electrons. Electrons inhuence each other only 

One solution to this problem is to use spin-unrestricted (UHF) theory, in which electrons are no longer paired within orbitals. Each electron has its own spatial orbital, leading to a better description of bond dissociation (Fig. 2, dash-double-dotted curve labeled UHF ). The UHF approach is sufficiently popular that it is the default for open-shell molecules (e.g., free radicals) in many software packages. However, 

Although the HF approximation is successful in many applications, it is usually inadequate for quantitative thermochemistry. Since it is a mean-held theory, it neglects the instantaneous repulsion between electrons, known as the electron correlation. Electron correlation is important in electron pairing and weak intermolecular forces. Thus, HF theory is inappropriate when electron pairs (such as bonding pairs) are broken apart, or when van der Waals interactions are important. More sophisticated theories include electron correlation in some manner. Most of these correlated theories are based on HF theory, so they are sometimes called post-HF theories. There are many. Only those most popular for thermochemical applications are described here. 

A special Cl calculation, called full configuration interaction (FCI), is the most accurate, and most expensive, for a given basis set. All possible electron configurations are included. This is equivalent to complete-order coupled-cluster theory. For an A-electron molecule, the FCI result is attained when Ath excitations are included in the Cl or CC calculation. Unlike truncated Cl, such as CISD, FCI is size-consistent. Continuing the 

---

It is important to note that the only physical approximation here is the LDA. All other approximations are of a numerical nature and their convergence can be monitored and improved in a systematic way. Thus, this method allows a probe of the LDA limit for molecules and clusters without any further approximations such... 

We note that the calculation of At/ will depend primarily on local information about solute-solvent interactions i.c., the magnitude of A U is of molecular order. An accurate determination of this partition function is therefore possible based on the molecular details of the solution in the vicinity of the solute. The success of the test-particle method can be attributed to this property. A second feature of these relations, apparent in Eq. (4), is the evaluation of solute conformational stability in solution by separately calculating the equilibrium distribution of solute conformations for an isolated molecule and the solvent response to this distribution. This evaluation will likewise depend on primarily local interactions between the solute and solvent. For macromolecular solutes, simple physical approximations involving only partially hydrated solutes might be sufficient. 

This expression separates solute-solvent interactions into an outer shell contribution, which we submit can be described using simple physical approximations, and an inner shell contribution, which we treat using quasi-chemical theory. 

Because QRRK theory was developed long before computing became readily available, it had to employ significant physical approximations to obtain a tractable result. The most significant assumption was that the molecule is composed of s vibrational modes with identical frequency i and that other molecular degrees of freedom are completely ignored. RRKM theory relies on neither approximation and thus has a much sounder physical basis. In the limit of infinite pressure, RRKM theory matches the transition state theory discussed in Section 10.3. 

Equation (7.8) offers a clear separation of inner-shell and outer-shell contributions so that different physical approximations might be used in these different regions, and then matched. The description of inner-shell interactions will depend on access to the equilibrium constants K. These are well defined, observationally and computationally (see Eq. (7.10)), and so might be the subject of either experiments or statistical thermodynamic computations. Eor simple solutes, such as the Li ion, ab initio calculations can be carried out to obtain approximately the Kn (Pratt and Rempe, 1999 Rempe et al, 2000 Rempe and Pratt, 2001), on the basis of Eq. (2.8), p. 25. With definite quantitative values for these coefficients, the inner-shell contribution in Eq. (7.8) appears just as a function involving the composition of the defined inner shell. We note that the net result of dividing the excess chemical potential in Eq. (7.8) into inner-shell and outer-shell contributions should not depend on the specifics of that division. This requirement can provide a variational check that the accumulated approximations are well matched. 

There are many ways one can try to reduce the computational burden. Ideally, one would find numerical methods which are guaranteed to retain accuracy while speeding the calculations, and it would be best if the procedure were completely automatic i.e. it did not rely on the user to provide any special information to the numerical routine. Unfortunately, often one is driven to make physical approximations in order to make it feasible to reach a solution. Common approximations of this type are the quasi-steady-state approximation (QSSA), the use of reduced chemical kinetic models, and interpolation between tabulated solutions of the differential equations (Chen, 1988 Peters and Rogg, 1993 Pope, 1997 Tonse et al., 1999). All of these methods were used effectively in the 20th century for particular cases, but all of these approximated-chemistry methods share a serious problem it is hard to know how much error is... 

Proton transfer at the surface of a protein or biomembrane is a cardinal reaction in the biosphere, yet its mechanism is far from clarification. The reaction, in principle, should be considered as a quantum chemistry event, and the reaction space as a narrow layer, 3-5 water molecules deep. What is more, local forces are intensive and vary rapidly with the precise molecular features of the domain. For this reason, approximate models that are based on pure chemical models or on continuum physical approximations are somewhat short of being satisfactory models with quantitative prediction power. 

In order to define the basis states, the interaction potential V, - is, as usual in atomic physics, approximated by a sum over individual potentials Vi, where Vi must be Hermitian, but can be chosen arbitrarily. We can thus define an unperturbed or independent particle Hamiltonian... 

The first class of approximations are referred to as 44 mathematical approximations because they are based upon purely mathematical considerations. The second class of approximations are called physical approximations because they are based upon the physical characteristics of a particular system. 

Before getting more involved in the finer points of particle adhesion studies by molecular mechanics it must be pointed out that such an approach suffers from considerable drawbacks. Molecular mechanics and dynamics by definition involve interactions between clearly defined types of atoms, with clearly defined atomic characteristics, placed in clearly defined molecular structures. Thus, a generalized, fictitious surface of only, let s say nitrogens, or even worse of generalized spheres, is a rather extreme physical approximation. It is by definition incorrect in a molecular mechanics and dynamics investigation. Such a drawback needs to be pointed out to put in perspective and understand... 

All other thermodynamic and dynamic properties determined using the XP3P model in Table 2.3 are in reasonable accord with experiments and are of similar accuracy in comparison with other empirical models. We note that in contrast to the large number of PMMFFs in the literature that are based on parameterization using different physical approximations, the electronic polarization from the present XP3P model is explicitly described based on a quantum chemical formalism. 

Since the exact solution of the SchrBdinger equation for most such processes of interest remains far beyond our current capabilities, both in terms of the algorithms for exact solution and the computational resources, both software and hardware, to carry them out, most of this Workshop has focused on the adequacy of various physical approximations... 

A theoretical model corresponds to a choice of physical approximations. Performing an actual calculation also requires numerical approximations. In particular, a basis set must be selected. This is the set of functions used to describe the molecular orbitals. A large basis set contains many basis functions and describes orbital shapes better than a small basis set. However, the improved accuracy of a larger basis set is counterbalanced by greater computational cost. The nomenclature and notation for basis functions may appear mysterious, but do have logical structure some are described below. 

If the configuration in question is a viable physical approximation (i.e. i), then the support of the wavelet transform must be consistent with the quantization scale. That is a, 6 (og, oo), or uq < a -... 

As noted previously, a similcir relation was used in the recent work of HMBB, which developed a TPQ analysis within the MRF representation, the physical solutions were identified by the extent to which the condition physQj- (aQ,cxi) was satisfied, where oquantization scale, flg, then the solution was considered as a viable physical approximation. 

The BLA in conjugated polymers can be very well described by an improved Hiickel-type method, the sem-iempirical Longuet-Higgins-Salem (LHS) model [7, 8]. The corresponding solid-state physics approximation is the Su-Schrieffer-Heeger model [9, 10]. The relationship between the LHS and the SSH models was pointed out in Ref. [11]. 

---