Numerical methods basis function

Two-dimensional semiclassical studies described in section 4 and applied to some concrete problems in section 6 show that, when no additional assumptions (such as moving along a certain predetermined path) are made, and when the fluctuations around the extremal path are taken into account, the two-dimensional instanton theory is as accurate as the one-dimensional one, and for the tunneling problem in most cases its answer is very close to the exact numerical solution. Once the main difficulty of going from one dimension to two is circumvented, there seems to be no serious difficulty in extending the algorithm to more dimensions that becomes necessary when the usual basis-set methods fail because of the exponentially increasing number of basis functions with the dimension. 

Here, n corresponds to the principal quantum number, the orbital exponent is termed and Ylm are the usual spherical harmonics that describe the angular part of the function. In fact as a rule of thumb one usually needs about three times as many GTO than STO functions to achieve a certain accuracy. Unfortunately, many-center integrals such as described in equations (7-16) and (7-18) are notoriously difficult to compute with STO basis sets since no analytical techniques are available and one has to resort to numerical methods. This explains why these functions, which were used in the early days of computational quantum chemistry, do not play any role in modem wave function based quantum chemical programs. Rather, in an attempt to have the cake and eat it too, one usually employs the so-called contracted GTO basis sets, in which several primitive Gaussian functions (typically between three and six and only seldom more than ten) as in equation (7-19) are combined in a fixed linear combination to give one contracted Gaussian function (CGF),... 

When the FMS method was first introduced, a series of test calculations were performed using analytical PESs. These calculations tested the numerical convergence with respect to the parameters that define the nuclear basis set (number of basis functions and their width) and the spawning algorithm (e.g., Xo and MULTISPAWN). These studies were used to validate the method, and therefore we refrained from making any approximations beyond the use of a... 

Since the exact solution of the Hartree-Fock equation for molecules also proved to be impossible, numerical methods approximating the solution of the Schrodinger s equation at the HF limit have been developed. For example, in the Roothan-Hall SCF method, each SCF orbital is expressed in terms of a linear combination of fixed orbitals or basis sets ((Pi). These orbitals are fixed in the sense that they are not allowed to vary as the SCF calculation proceeds. From n basis functions, new SCF orbitals are generated by... 

BSSE also affects the shape of the potential energy surface and the energy derivatives. There have been numerous attempts to find a general scheme to eliminate this error, and both a posteriori [2] and a priori [3] schemes are available. The counterpoise approach (CP) by Boys and Bemardi [5] and related methods are the most common a posteriori procedures. Within this method, the monomer electrons are described by the same basis functions as those used in the complex by means of the so... 

This representation permits analytic calculations, as opposed to fiiUy numerical solutions [47,48] of the Hartree-Fock equation. Variational SCF methods using finite expansions [Eq. (2.14)] yield optimal analytic Hartree-Fock-Roothaan orbitals, and their corresponding eigenvalues, within the subspace spanned by the finite set of basis functions. 

Variational one-center restoration. In the variational technique of one-center restoration (VOCR) [79, 80], the proper behavior of the four-component molecular spinors in the core regions of heavy atoms can be restored as an expansion in spherical harmonics inside the sphere with a restoration radius, Rvoa, that should not be smaller than the matching radius, Rc, used at the RECP generation. The outer parts of spinors are treated as frozen after the RECP calculation of a considered molecule. This method enables one to combine the advantages of two well-developed approaches, molecular RECP calculation in a gaussian basis set and atomic-type one-center calculation in numerical basis functions, in the most optimal way. This technique is considered theoretically in [80] and some results concerning the efficiency of the one-center reexpansion of orbitals on another atom can be found in [75]. 

Slater-type orbitals were introduced in Section 5.2 (Eq. (5.2)) as the basis functions used in extended Hiickel theory. As noted in that discussion, STOs have a number of attractive features primarily associated with the degree to which they closely resemble hydrogenic atomic orbitals. In ab initio HF theory, however, they suffer from a fairly significant limitation. There is no analytical solution available for the general four-index integral (Eq. (4.56)) when the basis functions are STOs. The requirement that such integrals be solved by numerical methods severely limits their utility in molecular systems of any significant size. 

This is the basis for least-squares methods of functional approximation, for example, and is certainly useful as a some measure of scatter in data points. However, as we shall see, the errors we must guard against in quantum chemistry tend to be systematic errors the failure of some approximation we have made. In this sense, it is equally important to know the worst-case error in our quantities. An appropriate analogy is with a numerical analyst writing a routine to evaluate some function. For a user, the important issue is how laxge is the maximum error when I use this routine . We should therefore consider also the maximum error... 

A recent numerical development is the introduction of the slow or smooth variable discretization (SVD) technique [101-103]. In the diabatic-by-sector method, the basis functions to expand the total wavefunction are fixed within each sector. In the SVD method, the hyperangular basis functions are constructed using the discrete variable representation (DVR) [104], The requirement is only that the total wavefunction be smooth in the adiabatic parameter p. By expanding the hyperradial wavefunctions using DVR basis functions, a new set of hyperangular basis functions are determined and they... 

---